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Measurements of emission wave-lengths in the Kf group of Ca and S in CaS and CaSO,, of 
Mg and S in Mg§, of S in SrS and BaS, of Ca and Cl in CaCl, have been made from the 
secondary radiation from pressed blocks of the salts. An attempt is made to determine the 
origin of as many x-ray lines as possible by means of electronic energy band systems in the 
cases of CaS, MgS, and ZnS. The K@; line of Zn in ZnS and the K@; line of Mg in MgS are 
found to be due to cross transitions from the outermost filled band of the negative ion. 


T has been pointed out that information in 
regard to electronic energy bands in solids 
may be obtained by a study of x-ray emission! 
and absorption? spectra. In the second paper of 
reference 1, it is shown that the principal lines in 
the x-ray emission spectra of alkali halides fit the 
energy level diagrams for the solid state with 
fair accuracy, leaving no doubt as to the levels 
and transitions which account for some of these 
lines. However, there are still several prominent 
lines which cannot be accounted for with cer- 
tainty. The present paper summarizes recent 
measurements of wave-lengths in the K£-radia- 
tion from alkaline earth sulphides and a few 
other compounds. The discussion reveals anoma- 
lies which were absent in the spectra of the alkali 
halides. 

The radiation was excited by placing pressed 
blocks of the compound near the focal spot on a 
tungsten-surfaced anticathode as described in a 
previous article.* An attempt was made at this 


1J. Valasek, Phys. Rev. 51, 832 (1937); Phys. Rev. 53, 
274 (1938). 

*R. Landshoff, Rev. Sci. Inst. 9, 143 (1938); Phys. Rev. 
55, 631 (1939). 

*J. Valasek, Phys. Rev. 47, 896 (1935). 


time to record as many of the satellites of K8, as 
possible with this indirect method of excitation. 
A Siegbahn vacuum spectrograph equipped with 
a quartz crystal was employed. The power input 
to the tube was kept low enough so that the 
material did not show any appreciable alteration. 
However, the material was changed every 24 
hours during the exposures, which ranged from 
100 to 500 hours. The power input was 1 to 3 ma 
at 12 kv depending on the stability of the 
substance. 

For the Mg Kf§,; line a focusing spectrograph 
with a curved mica crystal was employed. The 
radiation was from a pressed block of MgS 
which received radiation from an adjacent 
tungsten plate which was bombarded by elec- 
trons. The wave-length was measured relative to 
that of copper Ka; in the eighth and tenth 
orders as well as Al Ka in the second order. 

In the other cases, the wave-lengths were 
measured relative to those previously given for 
the principal Kf, lines' from the same com- 
pounds. The line intervals were measured with 
a Gaertner comparator. Each spectrogram con- 
tained a number of lines which were too weak 
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to measure directly since they disappeared into 
the granular background even under low mag- 
nification. When microphotometer traces were 
made of these lines, it was found difficult to 
locate accurately the exact maximum because of 
the line breadth and the accidental irregularities 
in the trace. It was found just as accurate to 
mark carefully the location of the line on the 
plate with a sharp pencil and to measure the 
mark. Repetition of this procedure on the same 
and on different plates gave results which were 
usually consistent to about 1 X.U. and the 
average should be good to about 0.5 X.U. as a 
rule. These lines are designated by b in Table I. 
The lines designated by a were always measured 
directly and are probably correct to 0.1 X.U. in 
general. 

The auxiliary data needed to assist in the 
interpretation of the lines are most complete in 
the case of CaS. In Fig. 1, the emission lines of 
Ca and S are plotted on a scale of electron volts 
below the K-absorption edge for the corre- 
sponding element. The wave-length of the K 
edge of S in CaS, 5006.6 X.U., was obtained 
from the measurements of Lindh.‘ That for the 
K edge of Ca in CaS, 3059.5 X.U., was obtained 
by applying a correction to the measurements of 
Aoyama, Kimura, and Nishina.’ These authors 
as well as Lindsay and Van Dyke® have measured 


TABLE I. Wave-lengths in X.U. 


LinE TyPE WaAVE-LENGTHS 
ion, S~~ (and Cl-) 
MgS aS SrS BaS CaSO, CaCl: 
5124.3 5083.5 5084.4 5093.4 5087. 
5067.3 5065.3 5065.3 5067.1 5067 
5046.2 5046 5046 4419.6 
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Positive ion, Ca++ (and Mg**) 


CaS CaSO, CaCl: MgS 
6 3122. 3122. 3125 
3111. 3110. 3114 
Bi std 43 3083. 43 pia 9545 
a@ 3079.2 3079.3 3079.3 
Bs @ 3067.3 3068.1 3068.1 
b 3048.3 3048.1 3048.1 
3041.3 3042.7 3041.2 


‘ A. Lindh, Handb. d. Exp. Physik, Vol. 24, No. 2, p. 296. 
5S. Aoyama, K. Kimura, and Y. Nishina, Zeits. f. 
Physik 44, 810 (1927). 
as A Lindsay, and G. Van Dyke, Phys. Rev. 28, 613 


the K-absorption wave-lengths of Ca in various 
compounds. There is an average difference of 
3.1 X.U. between the two sets of measurements 
on the same compound, although Lindsay and 
Van Dyke give no values for CaS. The difference 
between the two groups of observers is largely 
due to the fact that the former measured the 
middle of the first white line in the absorption 
spectrum, while Lindsay and Van Dyke meas- 
ured the beginning of absorption. In view of both 
sets of data, it is considered that the most 
probable value of the K-absorption edge of Ca 
in CaS is determined by the most rapid rise in 
absorption which is estimated to be 1.6 X.U. 
above the value given by Aoyama, Kimura, and 
Nishina, namely at 3059.5 X.U. 

In order to fit the energy levels of both ions 
to the same energy scale, it will be assumed at 
first that the K electron, at the beginning of the 
continuous K absorption, is raised to the lower 
edge of the same band of empty states of the 
crystal lattice regardless of whether it comes 
from the K shell of the positive or the negative 
ion. Accordingly, when they are both known, 
the absorption edges of both ions are placed at 
the same point in the energy level diagram. 
The locations of the emission lines should then 
show the locations of the initial occupied level 
from which the emitting electron falls into the K 
shell. This assumes a single excitation and a 
single electron transition, and hence does not 
apply to the spark lines. 

In the cases of NaCl and KCl reported in a 
previous article, the origin of four of the emission 
lines was plainly indicated when a diagram pre- 
pared as above was compared with an energy 
level diagram prepared from the atomic levels by 
the method of Slater and Shockley.” Two of 
these lines were, as expected, due to transitions 
of a 3p electron into a K shell vacancy in the 
same ion, while two were due to transitions of a 
3p electron into a K shell vacancy in a neighbor- 
ing ion of opposite sign. The latter are usually 
much weaker than the former except when the 
3p, or sometimes the 2p, shell is not occupied in 
one of the ions because of its low atomic number. 
Then the cross transitions give rise to principal 
lines. Further evidence on this point is presented 
below. 
~ TJ, Slater and W. Shockley, Phys. Rev. 50, 705 (1936). 
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ENERGY BANDS 


Because the alkaline earth sulphides also have 
a rocksalt type of structure, one can readily 
apply the same calculations to find the locations 
of their energy bands for the purpose of com- 
parison with similar results obtained from the 
x-ray data. Because of the double charge on each 
ion, the electrostatic potential energy term must 
be multiplied by two. Accordingly 


3.495 e* 300 


eo = electron volts, 


r 
where r is the distance between closest neighbors 
of opposite sign. The ionization potential, or 
electron affinity, of the S~~ ion is calculated by 
extrapolating the ionization potentials of the 
isoelectronic sequence Catt+, Kt, A, Cl- by 
means of a quadratic formula.’ The potentials 
used are taken from Bacher and Goudsmit’s 
Atomic Energy States. The work required to 
remove one electron from S~~ is thus found to 
be —4.18 electron volts. Seitz® has recently used 
a value of —5 volts which he obtained by sub- 
tracting an estimated binding energy for the 
first electron from the electron affinity found by 
the use of a Born cycle. The cycle gives the work 
done to remove both electrons. The negative 
sign applied to the electron affinity indicates that 
the free doubly charged ion is unstable. In a 
crystal lattice, however, more than enough bind- 
ing energy is supplied through the external 
screening by surrounding ions to result in a 
positive binding energy of some 10 to 15 electron 
volts, depending on the compound. 

The ionization potentials for the doubly 
charged positive ions of Mg and Ca are also 
from Bacher and Goudsmit’s tables. However 
those for Sr++ and Bat+ are not listed, and 
since the K-absorption edges of these ions in 
the corresponding sulphides are also unknown, 
they are omitted from the following discussion. 

Calculated positions of the energy bands for 
CaS are shown in Fig. 1. It is difficult to deter- 
mine the extent of the continuum of empty 
states, but the experimental data can be used 
to supply this information. In order to do this, 
it is assumed that the intervals between the 
occupied energy bands of the positive ion are 


*G. Glockler, Phys. Rev. 46, 111 (1934). 
*F. Seitz, J. Chem. Phys. 6, 454 (1938). 


OF COMPOUNDS 


Colcium in 


Calevieted Energy Bands in CaS 


3432 3226 B20 E420 
Energy— Electron volte below K ebsorption limit 


Fic. 1. Top: K-series emission lines of calcium in CaS on 
energy scale relative to the K-absorption edge. Center: 
~ 7. aman in CaS. Bottom: electron energy bands 


not altered by chemical combination. This is 
indicated by the exceedingly small effects of such 
combination on the emission wave-lengths of the 
positive ion. Hence, if the 3p level of Ca is 
placed at the same point on the energy scale as 
the Kf; line, the location of the K-absorption 
edge will give the location of the lower edge of 
the empty states in the compound. The other K 
lines should then indicate the upper filled states 
from which they originate. The Ca K&; line is 
thus seen to be due to a cross transition of an 
S 3p electron into the vacancy in the K shell 
of Ca. 

However, it is evident on Fig. 1 that the 
placing of the K-absorption edge of the negative 
ion in coincidence with that of the positive ion, 
results in a large discrepancy in the location of 
the S K®; line with respect to the S 3p band. It 
should be noted that this displacement did not 
occur as a result of the application of the same 
considerations to the spectra of NaCl and KCl 
in an earlier article.! Obviously, one can remove 
the discrepancy in this case by an appropriate 
shift in the position of the empty states belonging 
to the sulphur ion. This would require a dis- 
placement of 5 to 7 electron volts with selection 
rules to reserve the shifted states for the negative 
ion. The existence of two such empty bands 
“belonging” to the respective ions does not seem 
improbable, for an examination of the wave- 
mechanical problem leads to a result of this 
general nature. When K excitation of the positive 
ion takes place, the K electron goes into a band 
derived from the lowest vacant p state of the 
positive ion. On the other hand when the K 
excitation of S-~ takes place, the electron 
probably goes over into the vicinity of one of 
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10 L, Pauling, Phys. Rev. 34, 954 (1929). 


the neighboring Ca** ions. This is essentially the 
same state of Ca, but modified by the fact that 
the net charge on the ion core is one unit less. 
An application of Moseley’s law as a first approxi- 
mation gives much too small a displacement of 
the upper level to account for that observed. 
Thus it seems necessary to assume a great distor- 
tion of the outer levels of the negative ion either 
as a result of the nature of the chemical bond, 
or as a result of some secondary process which 
occurs between excitation and emission. 

The first interpretation is ruled out because 
the cross-transition line KB; of Ca shows that 
the initial level S~— 3 is in its normal state when 
it is the positive ion which is excited. This con- 
clusion is supported by similar lines in the spectra 
of MgS and ZnS as will be shown in the follow- 
ing. Thus the displacement of the S~— 3p level 
accompanies only the K excitation of the S 
ion itself and, by the Franck-Condon principle, 
should follow the excitation. Accordingly, the K- 
absorption frequency of S~~ will be at its ex- 
pected place in the compound but its K8-emission 
line will be displaced. 

Following Gurney, one may assume two prin- 
cipal sources of such secondary effects on the 
energy levels, namely the displacement of sur- 
rounding ions and polarization. The net effect is 
to elevate the energy levels of the electrons in 
the K-excited ion whether it is positive or 
negative, that is to say, to decrease the binding 
energy. Aoyama, Kimura, and Nishina® have 
calculated the energy wuz of rearrangement of 
surrounding ions when the central ion loses an 
electron. They used this as one of the terms in 
their explanation of the effect of chemical com- 
bination on absorption frequencies. Pauling’® 
points out that this is incorrect because the 
rearrangement cannot take place during absorp- 
tion, but must follow after it, and he substitutes 
an electron affinity of the crystal which he 
obtains from the diamagnetic susceptibility of 
the compound. 

Aoyama, Kimura, and Nishina find that the 
energy uz is only 1 electron volt in the case of 
chlorine in sodium and potassium chloride, but 
becomes 10 electron volts in the case of sulphur 
in CaS, and 11 electron volts for sulphur in 
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MgS and ZnS. Furthermore, it is very small in 
the cases of the positive ions in all these sub- 
stances. From the manner in which wz is calcu- 
lated, this energy includes energy of polarization 
as well as displacement, though the former is not 
mentioned in the original article. It is question- 
able whether the lifetime of the K-excited states 
is large enough to permit ion displacements, but 
polarization should be considered and it is 
estimated that this will lead to a displacement 
of the valence levels of the central ion by some 
5 to 6 electron volts. The discussion by Aoyama, 
Kimura, and Nishina shows that a similar effect 
escapes notice in the alkali halides and in the 
spectra of positive ions in general, because 1, is 
relatively small in these cases. 

A further indication that the energy levels of 
the negative ion can be shifted by several electron 
volts as a result of changes in the internal and 
external screening, is shown by the spectra of 
sulphur in the sulphates as compared with the 
corresponding sulphides.! The K-absorption edge 
of sulphur is found to be displaced toward 
higher energy values by 9.5 electron volts on 
the average. The Ka lines are displaced to higher 
frequencies by 1.7 volts, showing that the L 
levels are also displaced considerably, the shift 
being 7.8 electron volts. On the other hand, the 
KB lines are displaced by only 0.5 volt on the 
average. This shows an unusually large amount 
of distortion or relative displacement of the 
energy levels. 

The above considerations show that energy 
level displacements of the indicated value can 
occur when the region around a sulphur ion is 
made more negative. The principal displacement 
applies to all levels equally, like the crystal 
potential effect previously considered. However 
a secondary effect of opposite sign applies 
principally to the outermost levels and makes for 
a decrease in binding. This differential shift takes 
place between excitation of the negative ion and 
emission by the same ion and is negligible in the 
case of the positive ion. Any influence which acts 
equally on the M and K levels cannot explain 
the effects of chemical combination on the 
observed emission wave-lengths. It is just this 
point which makes it difficult to calculate the 
exact magnitude of the effect beyond the rough 
approximation discussed above. One might ex- 
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ENERGY BANDS OF COMPOUNDS 
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Magnesium in MgS 


Sulphur in Mgs 


| 7 


64 60 56 S2 48 444036 32 28242018 128 40 
Electron Volts 


Fic. 2. Top: K-series emission lines of Mg in MgS on 
energy scale. Center: ditto, sulphur in MgS showing 
relative location of K-absorption edge. Bottom: electron 
energy bands of MgS. 


pect that it would increase with the polarizability 
of the ion concerned and the valence of the sur- 
rounding ions and their proximity. This is in 
accord with the observations, for the effects are 
greatest in the doubly charged negative ion in a 
divalent lattice with small lattice constant. 

In the case of the spectrum of MgS, the wave- 
length of the K-absorption edge of Mg in the 
compound is not known, making it impossible to 
link together the spectra of S~~ and Mgt* by 
making the vacant bands coincide as was done 
in the case of CaS. However, the locations of 
the energy levels in the crystal can be estimated 
as before by applying the crystal potential term 
to the levels for the free ion. The lower limit of 
the continuum of vacant states can be estimated 
to lie at about 2 electron volts by analogy with 
CaS, or by calculating the average energy of a 
free electron in the crystal from a formula given 
by Pauling.'° Figure 2 shows the resulting 
diagram. The K edge of sulphur is made to 
coincide with the lower limit of the continuum 
and the magnesium Ka line is placed opposite 
the Mg++ 2p level. There is no 3p level in a free 
magnesium ion from which a Kf line might 
originate. Nevertheless such a line is observed. 
Its location in the energy level diagram shows 
that it must be interpreted as a cross transition 
from the S~~ 3p band. The same interpretation 
applies to the K; line of Na in NaCl. The dis- 
placement of the S Kf; line with respect to the 
parent S~— 3p energy band is about 6 electron 
volts which is in concordance with the result 
noted for CaS. 

Recently F. Seitz® has published an energy 
level diagram for ZnS which is reproduced in 


Zine in 


Sulphur in 
Be 


Caleviated Energy Bands in ZnS (Seitz) 


98 96 94 20 
Volts 


Fic. 3. Top: K-series emission lines of zinc in ZnS on 
energy scale. Center: ditto, sulphur in ZnS showing rela- 
tive location of K-absorption edge. Bottom: electron 
energy bands of ZnS according to Seitz. 


Fig. 3 with x-ray emission data for Zn" and S® 
in the same compound. The K-absorption edge 
of S in ZnS is taken from the tabulation by 
Lindh, but the wave-length of the K edge of Zn 
in this compound could not be found in the 
literature. When the K-absorption edge of S is 
lined up with the lower edge of the continuum 
given by Seitz it is seen that the S KA; line falls 
nearly seven electron volts above the S 3p band. 
This is the same as the effect observed in CaS 
and MgS and is explained in the same way. It 
should be noted that in the S KB group, the 
strongest line is taken to be the S K®; line and 
that it is accompanied by weaker satellites on 
the long wave-length side instead of the short 
wave-length side as in the alkaline earth sul- 
phides. This line structure and intensity relation 
was observed in zinc blende and wurtzite as well 
as in pellets pressed from the chemically pure 
reagent. 

To locate the Zn K&B; line correctly, the Zn 3p 
band was taken as the origin of the Zn K§; line. 
The Zn 3p level was located relative to the Zn 3d 
band by a combination of x-ray data," without 
the use of Zn K8;, the origin of which is doubtful. 
The equation used was: 


Mw,v— 


the La wave-length being obtained from the 
measurements of Tyren." The interval between 
Zn 3d and Zn 3p was found to be 76.0 electron 
volts. When the Zn Kf, line is placed in coinci- 
dence with the Zn 3p band, the K@; line falls at 


"J. C. McDonald, Phys. Rev. 50, 694 (1936). 
2 J. Valasek, Phys. Rev. 43, 612 (1933). 

Siegbahn, S yy (1931). 
(19s + Tyren, Archiv f. t., Astr. och Fysik 25A, 1 
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S-- 3p as shown in Fig. 3. This indicates that 
the Zn KB; line which is usually regarded as due 
to a “forbidden” or quadrupole transition (Al = 2) 
is due instead to a cross transition of a S~~ 3p 
electron into a K-shell vacancy in the zinc ion. 
One can thus understand the relatively low in- 
tensity of this line and its anomalous variation 
with atomic number. The above applies to zinc 
sulphide and probably to all salts of elements in 
the iron group. The Kf; line has also been 
observed in the x-ray emission spectrum of 
metallic zinc, in fact, accompanied by a weak 
and scarcely resolvable Kf, line.!* These lines 
show a breadth and lack of symmetry which is 
characteristic of the conduction band in metals. 


% J. A. Bearden and C. H. Shaw, Phys. Rev. 48, 18 
(1935). 


PERLOW 


Their origin is discussed by Beeman and Fried- 
man!® in connection with their study of the 
structure of the K-absorption edges of com- 
pounds in the neighborhood of the iron group. 
They also give excellent reasons against the usual 
interpretation of K®; as a forbidden or quadru- 
pole transition in those cases where there is 
conduction band formation in the solid. The 
present observations give support to a similar 
interpretation in regard to compounds. 

The writer takes this opportunity to thank the 


Graduate School of the University of Minnesota . 


for its support of this investigation, and to 
Mr. George Werner for his assistance in the 
laboratory. 


16 W. W. Beeman and H. Friedman, Phys. Rev. 56, 392 
(1939). 
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An Accurate Measurement of the Energy Released in the 
Disintegration of Li® by Protons 


GILBERT J. PERLOW 
Ryerson Physical Laboratory, University of Chicago, Chicago, Illinois 
(Received May 21, 1940) 


The energy release in the reaction Li*(p,~)He*® has been determined by a comparison of the 
He‘ range with that of the alpha-particles from Be*(p,a)Li*. The energy of the latter reaction 
is accurately known from electrostatic deflection measurements. The determination was carried 
out in a small, variable air pressure, absorption cell and ionization chamber, both being at 
the same pressure. The device had high dispersion, 1 mm of range corresponding to about 13 
cm Hg. Calibration of the cell was carried out by a method which utilized this high dispersion. 
The Q value for the reaction was found to be 3.945+0.06 Mev as compared to the previously 
accepted value due to Neuert of 3.72+0.08. With the new Q value, H? is found to be unstable 
against §-decay by either 0.19+0.09 or 0.10+0.1 Mev depending on whether Bainbridge’s or 
Aston’s value for the deuterium-helium bracket is used. The mass of He® is found to be 


3.01688 +0.00011 mass units. 


INTRODUCTION 


NLY three nuclear reactions are known at 
the present time which lead to the forma- 
tion of the He* nucleus. The one most studied 
has been D?+D*—He’+mn!, and this chiefly 
through measurement of the neutron energy ;} 
the alternative being to measure a 4-mm range. 


1T. W. Bonner and W. M. Brubaker, Phys. Rev. 49 


19 (1936); T. W. Bonner, Phys. Rev. 52, 686 (1937) a 
53, 711 (1938); E. Baldinger, H. Huber, and H. Staub, 
Helv. Phys. Acta 11, [3] 245 (1938), and others. 


Another reaction is the recently discovered? beta- 
radioactivity of H*. The third is the reaction 
Li®+H'!—He'’+ He’, which is the subject of this 
paper. These short range particles were first 
found by Oliphant, Kinsey, and Rutherford.’ The 
best measurement of the energy release hitherto 
published is that of Neuert,‘ who measured the 


2L. Alvarez and R. Cornog, Phys. Rev. 56, 613 (1939) 
and 57, 248A (1940). 

3M. L. E. Oliphant, B. B. Kinsey, and Lord Rutherford, 
Proc. Roy. Soc. 141, 722 (1933). 

*H. Neuert, Physik. Zeits. 36, 629 (1935). 
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range of the He® nuclei in a cloud chamber. This 
method is open to criticism chiefly because it 
depends directly on the range-energy curve in a 
region where it is not too well known. The method 
to be described here depends only on the slope 
of the range-energy curve over a narrow region. 
The percentage error in the Q value is 0.2 times 
the percentage error in the slope. 

In this work the He‘ range from Li*(p,a)He* 
was compared by an absorption cell and ioniza- 
tion chamber method with that of the He*’s 
emitted in Be®(p,a)Li®. The energy of the latter 
had been accurately determined in this labora- 
tory by electrostatic deflection,’ and since the 
expected difference in range was only of the order 
of a millimeter, it was felt that the measurements 
would be of high precision. The ranges were 
expected to be about 7 and 8 mm for the beryl- 
lium and lithium reaction, respectively. Thin 
targets were used in both cases. 


APPARATUS 


The high voltage set used was the Cockcroft- 
Walton type circuit described by Hatch.* A low 
voltage ion source designed by Allison furnished 
the protons. The beam was magnetically sepa- 
rated and passed through two circular apertures 
0.54 cm in diameter and 10 cm apart before 
striking the target. The disintegrations were ob- 
served at 90°. Figure 1 is a drawing of the target 
assembly and range measuring apparatus. The 
latter is seen to consist of two sections. The 
region from B to A acted as a variable air pres- 
sure absorption cell, that from A to the collecting 
electrode farther to the right as a parallel plate 
ionization chamber. The target spot and the 
opening A were the limiting apertures, each 0.54 
cm in diameter. The screen B was covered on the 
higher pressure side by a collodion window with 
an air equivalent of about 3.5 mm. This was the 
only window in the apparatus. There was a fine 
perforated screen over A also, but nothing else, 
so that the pressures in the ionization chamber 
and absorption cell were the same. The screen 
had 26 holes per sq. cm, each hole being about 


5S. K. Allison, L. S. Skaggs, and N. M. Smith, Jr., 
Phys. Rev. 54, 171 (1938); L. S. Skaggs, Phys. Rev. 56, 24 
(1939); S. K. Allison, L. S. Skaggs, and N. M. Smith, Jr., 
Phys. Rev. 57, 550 (1940). 

*G. T. Hatch, Phys. Rev. 54, 165 (1938). 
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Fic. 1. Target assembly and range measurement ap- 
paratus. There was a collodion window on the right-hand 
side of screen B. 


0.02 cm in diameter. The depth of the absorption 
cell was 0.54 cm and of the ionization chamber 
0.45 cm. The collecting voltage for the latter was 
applied to the outside brass case of the entire 
device. 


PERFORMANCE OF THE VARIABLE PRESSURE 
RANGE COMPARISON DEVICE 


The performance of the variable pressure 
ionization chamber will depend on the biasing of 
the counting circuit, that is, on the minimum 
number No of ion pairs which must be formed 
between the plates to record a count. Preliminary 
trials showed that acceptable number-pressure 
curves of the short range particles from beryllium 
could be produced with a bias of 9 volts on the 
first stage of the vacuum-tube scale-of-sixteen 
counting circuit. By a method described later it 
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CELL TEMPERATURE 
222°C. 
BIAS 9.0 VOLTS 


[ CMS. PRESSURE BELOW A 


CALIBRATION OF CELL 


TEMPERA 
23.5° C. 
BIAS 9.0 VOLTS 


35 


10 


Fic. 2. Calibration of the cell with polonium alpha-particles. Curves A and C are separated by 1.5-mm 
motion of micrometer at the prevailing barometric pressure and at the temperatures indicated in the 
drawing. Curves B and C represent the effect of changing the bias three volts for groups of the same 


was shown that this corresponded to the collec- 
tion of 1900 ion pairs. 

The range of a particle measured in the cell 
is a sum of three terms. The first is the air 
equivalent of the collodion window. This is as- 
sumed to be independent of the temperature and 
of the energy of the particle. The second term 
refers to the range expended in the first compart- 
ment of the device, i.e., the absorption cell 
proper. This will be the product of the cell length 
L with the air density p. In order to conform to 
the standard practice in range work, the unit of 
density is taken to be the density of dry air 
under 76 cm Hg pressure at a temperature of 
15°C. The third term is the product of the pene- 
tration 6 into the ionization chamber required 
to register a count with the corresponding air 
density. In the device used here, the latter is 
also p. Then for the difference in range of two 
separate groups of particles we have 


Ri — R2=L(p1— p2) + (5101 — 522). 


If now neither 6; nor 62 is greater than the 
ionization chamber depth D (=0.45 cm), the 
difference (51p:— 452p2) is zero since 6; and 42 are 
the linear path lengths required to produce Ny 
ion pairs before losing ionizing power and these 
lengths vary inversely as the density. The device 
now acts like the conventional type in which the 


residual range. The ordinate 10 represents 1000 counts in 70 seconds. 


ionization chamber is at constant pressure. Then 


288L 
76.0 LT, Ts 


where P is the pressure and 7 the absolute 
temperature. 

There is a certain critical pressure P., however, 
below which formula (1) becomes invalid. This is 
the pressure at which an alpha-particle will not 
produce enough ionization to record a count 
before reaching the end of the ionization cham- 
ber. It may be computed from: — 


Mo= f 


with D.= DP./76. The abscissa x = 0 corresponds 
to residual range zero in the specific ionization 
versus range curve I(x). If it is assumed that I(x) 
has the form of Fig. 9 in the paper of Holloway 
and Livingston,’ and that the maximum ordinate 
represents the creation of 60,000 ion pairs per cm, 
one gets P.=17 cm Hg for No=1900 ion pairs. 
Since none of the number-pressure curves started 
dropping off until somewhat above this pressure, 
formula (1) was applicable. 


7M. G. Holloway and M. S. Livingston, Phys. Rev. 54, 
18 (1938). 
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The quantity No was determined in the follow- 
ing way. A thin polonium source was placed so 
that its alpha-particles passed through the ab- 
sorption compartment and entered the ionization 
chamber near the end of their range. The pres- 
sure was then reduced in the two-compartment 
device and the diminishing pulse heights on the 
oscillograph observed. (In order to reduce the 
pressure without rupturing the collodion film, a 
vacuum was created in an auxiliary chamber, 
temporarily attached.) The pressure was lowered 
until this rough evidence indicated that the 
particles traversing the ionization compartment 
were producing approximately the critical num- 
ber of ion pairs corresponding to the selected bias. 
With the pressure held at this value, the distance 
of the polonium source was readjusted until the 
pulses showed a maximum height. Then an ac- 
curate measurement of the critical pressure was 
made by counting with the biased circuit and 
varying the pressure in this reduced region until 
the sharp low pressure cut-off of the number- 
pressure curve was found. Let the pressure at 
which the cut-off occurs be Pp. If the linear depth 
of the ionization chamber is D, then its equiva- 
lent depth at Po is PoD/76=Dp. The critical 
number of ion pairs No corresponding to the 
given bias may then be found from: 


a+Do 


where a is to be chosen to give the greatest value 
for the integral. For the 9-volt bias, Py) was 
measured to be 5.5 cm Hg. For D=0.45 cm Dy 
becomes 0.033 cm, and graphical integration 
gives Ny= 1900 ion pairs. 


CALIBRATION 


Since the difference in range of two groups of 
particles as measured with an absorption cell is 
proportional to the product of the cell length by 
the pressure difference between the two number- 
pressure curves, the dispersion of the device may 
be made extremely large by choosing a short 
enough cell. Such was the case in the present 
measurements, each mm of range in standard air 
corresponding to about 13 cm difference in 
pressure. Now measurement of the cell dimen- 
sions is complicated by such things as the bulge 
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of the screens, and the uncertainty of the bound- 
ary between absorption and ionization compart- 
ments. It is clear that a small absolute error in 
this measurement produces a large error in the 
range difference. Therefore the measured cell 
dimensions were given no weight in the calibra- 
tion. Instead a method suggested to the author 
by Dr. N. M. Smith of this laboratory was em- 
ployed. A thin polonium source was mounted on 
an accurate metric micrometer. The axis of the 
range measuring device was placed along the 
screw axis and the micrometer adjusted until 
particles ended their range in the cell at some 
particular pressure. The evacuated chamber 
mentioned previously was used to avoid ruptur- 
ing the cell window.A number pressure curve was 
taken for a given position of the source. The 
micrometer was then turned an integral number 
of times and another number-pressure curve 
taken. In Fig. 2 curves A and C are separated by 
1.5 mm motion of the micrometer at the prevail- 
ing barometric pressure. Table I gives the data 
of several such calibration runs suitably cor- 
rected for variations in temperature and baro- 
metric pressure during the run. The mean value 
is seen to be 12.95+0.2 cm Hg difference in 
pressure at 15°C for each mm difference in range 
in air at 76.0 cm Hg and 15°C. 

This method of calibration gives directly the 
range difference of two groups of alpha-particles 
measured with the cell. It automatically takes 
into account the irregularities of the screens, the 
variation of the air equivalent of window with the 
alpha-particle energy, any deviation from uni- 
formity of the collecting field in the ionization 
chamber, the protruding of the collodion window 
through the screen apertures and any variation 
of this with pressure, and finally it minimizes 
any correction to the range difference due to 
change in obliquity of the two groups,—for al- 
though the geometry for the calibration was not 


SIGNAL 


OPERATING 
REGION 


NOISE ( X~20) 


STAGES 1&2 


Fic. 3. Some operating characteristics of the amplifier and 
ionization chamber. 
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identical with that for the reactions investigated, 
it was not very different. 

The geometry for the experiment was ‘‘good”’ 
according to the criteria of Livingston and Bethe® 
whose notation is used below. The ratio of the 
limiting aperture widths to the separation be- 
tween them (one of them being the target spot) 
was a/b=0.12 while the quantity cos 0) was equal 
to 0.26 and 0.18 for the lithium and beryllium 
reactions, respectively. The obliquity correction 
was found to be less than 0.01 mm and was 
neglected. The angle straggling should be sym- 
metrical and therefore not affect the mean range. 

Pressure in the cell was measured with an 
open-tube mercury manometer. Atmospheric 
pressure was measured with an aneroid barom- 
eter frequently checked against the standard 
barograph at the U. S. weather bureau observa- 
tory a few hundred yards away. The temperature 
of the cell was read from a thermometer pressing 
against it. This temperature was usually about 
2°C higher than that of the room because of heat 
conduction from the target chamber. Diffusion 
of cooler air into the cell was impeded by making 
the pressure lead of small bore. 

Lithium targets were prepared by deposition 
of the fluoride from a flame upon a nickel button. 
For the beryllium targets, the fluoride was de- 
posited from an oxygen flame. Because of the 
high yields from the beryllium reaction,® it was 
found necessary to rub the targets clean before 
the intensity was reduced enough to enable ac- 
curate measurement. A high intensity of par- 
ticles produces an apparent shortening of the 
range as well as some straggling, due to particles 
crowding up on the “back-kick’”’ of a given 


amplifier pulse. 


Targets were heated after the fashion de- 
scribed by Hatch.* It is certain that the beryllium 
targets remained free of carbon. For the lithium 
ones, however, it was difficult to decide since 
alkali halides become discolored by bombard- 
ment with charged particles. It may have been 
the case that cohesive forces between the lithium 
fluoride and such organic materials as pump oil 
were stronger than those due to heat agitation. 
In such a case, a film of oil would be left, and 
would subsequently carbonize. This would reduce 


8M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 


245 (1937). 


the measured range, and because of the irregu- 
larity of the carbonization, give increased strag- 
gling. To guard against this, no lithium target 
in the data reported was used for more than one 
run. However, bombardment of badly discolored 
targets showed no unidirectional shift of the 
mean range point, nor any appreciable change 
in the slope of the number-pressure curves. 

The linear amplifier constructed for this re- 
search was of the usual resistance coupled type. 
A 954 acorn pentode was used in the first stage. 
There were two subsequent stages using 6]7 
pentodes followed by two employing 6CS triodes, 
Power for these was furnished by a power pack, 
except that the heater on the second stage was 
run in parallel with that on the first, a storage 


‘ battery supplying the current. All stages but the 


last were self-biased. The power supplies for the 
amplifier and ionization chamber were of a 
simple sort, stabilization being achieved by a 
voltage regulator of the saturated core type 
placed on the a.c. line. Figure 3 gives some of the 
characteristics of the amplifier and ionization 
chamber. The curves on the left show the de- 
pendence of the gain and noise upon heater 
voltage of the first two stages. It is seen that 
there was a broad plateau in the neighborhood of 
6 volts. This coincided with the discharge plateau 
of the storage battery, making stabilization auto- 
matic. The curve on the right shows that the 
ionization chamber was run at saturation. A 
scale-of-sixteen circuit similar to that described 
by Stevenson and Getting® was utilized. It was 
likewise stabilized by way of the a.c. line. The 
bias voltage was supplied by a battery and read 
with a good voltmeter. 


MEASUREMENTS 


All data were taken with the bombarding 
beam at 370 kv as measured with a high resist- 
ance voltmeter calibrated some time before. This 
was safely above the observed 330-kev resonance 
in F'*, The deflecting magnet was set on the 
maximum of the mass 1 spot. Considerations of 
the ripple in a Cockcroft-Walton type circuit 
show that most of the beam has the average 
voltage, so that the energy of the protons at the 
target could be assumed quite close to 370 kv. 


*E. C. Stevenson and I. A. Getting, Rev. Sci. Inst. 8, 
414 (1937). 
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Fic. 4. Data for the computation of the Q value in Li*(p,2)He*. Both the beryllium and lithium 
curves as plotted are corrected for barometric variation. Correction for the change in bias of the two 
curves has not been made. The extreme right of the beryllium curve shows the deuteron cut-off from 
Be*(p,d)Be*. The line marked 1 mm is from the calibration data of Table I. The numbers above the 
lithium curve indicate where the mean range point would have to occur to give various Q values if the 
alpha-particle energy from Be*(p,«)Li® is known exactly. 


Care was taken to keep the magnet current con- 
stant during a run, and both the magnet and the 
generator supplying it were allowed to warm up 
for more than an hour before data were taken. 
Beams of between two and fifteen microamperes 
of protons were used, the measurement being 
made with a current integrator. 

A considerable amount of adjustment of the 
absorption cell constants and the amplifier gain 
had to be made. Because of the high dispersion 
of the cell, and because it ran at pressures below 
atmospheric, only a small range interval could 
be measured. This interval had to be chosen to 
include the curves of interest and to exclude 
scattered protons (range~4 mm). It was also 
necessary that the deuterons from Be*(p,d)Be*® 
should have a threshold pressure sufficiently low 
that they might all be counted with the alphas 
from bery'lium. It was found (Fig. 4) that the 
deuterons had more range in the cell than the 
alpha-particles and contributed more ionization 
at all pressures. Pains had then to be taken not 
to bias out the alphas. A satisfactory solution 
was achieved by increasing the amplifier gain and 


decreasing the bias as much as could be done 
safely. The curve of Fig. 4 nrarked Be shows the 
alpha-particles and the deuterons. This is taken 
at such a bias that at least 1900 ion pairs must be 
produced to record a count. The yield from the 
beryllium targets was so great that incidental 
noise had a negligible effect on these curves. For 
the lithium reaction, however, the situation was 


quite different. The only isotopic Li® target avail- — 


able was found to be too thick. The yield from 
the thin fluoride targets was quite small and 
comparable to the number of small disturbances. 
It was therefore necessary to count at a higher 
bias. The bias could not be made too much 
higher or the low pressure cut-off for the long 
range alpha-particles from Li’ might fall in the 
pressure range being investigated and produce a 
rapid change in the background. A bias three 
volts higher (i.e., 12 volts) was found to be 
suitable. With this the He® particles and long 
range He*’s were counted at all pressures. Since 
the scaling circuit counted pulses of zero height 
for a 5.5-volt bias and 9.0 volts corresponded 
to 1900 ion pairs, 12 volts corresponds to 
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1900(12.0 —5.5)/(9.0 —5.5) or 3500 ion pairs. For 
this the critical pressure P. is computed to be 
25 cm Hg, which is safely below the region of 
interest. 

In Fig. 4 both the beryllium and lithium curves 
are shown. The mean range of the lithium had to 
be corrected for the three-volt change in bias. 
This was done with the calibration apparatus, 
curves being taken at both biases without moving 
the micrometer. B and C of Fig. 2 are such a pair. 
Two runs judged to be of equal weight gave 5.80 
and 6.00 cm Hg difference in pressure at 15°C for 
a three-volt bias change when suitable correc- 
tions were applied. The mean is then 5.90+0.04. 

An interesting point concerning the lithium 
reaction is that there were not as many He® as 
He! particles observed. This is explicable if the 
reaction does not have spherical symmetry ; for 
90° in the laboratory system corresponds to 
different angles in the center of mass system for 
the two kinds of particles. A curve taken by 
varying the bias showed that there were about 
as many He’’s from Li® as long range He*’s from 
Li’ and about as many short range He*’s as both 
together. A curve published by Neuert‘ shows 
the same inequality. This circumstance was quite 
helpful in that it reduced the probable error of 
the base-line on the number-pressure curves. 

The energy of the alpha-particles from beryl- 
lium was computed from the 1940 revision of the 
electrostatic analyzer data of Allison, Skaggs, 
and Smith.’ With a Q value of 2.115+0.04 
Mev, the Be alpha-particle energy becomes 
1.455+0.024 Mev at 370 kev proton energy. 
From Table II the pressure separation of the 
lithium and beryllium curves is 17.75+0.4 cm 
Hg. The bias correction is 5.90+0.04 cm Hg and 
addition to the uncorrected separation gives a 
true separation of 23.65+0.4 cm Hg. Table I 
gives the cell calibration as 12.95+0.4 cm 


TABLE I. Calibration of the absorption cell. 


A4P/AR 


MOTION OF cm Hg 
TRIAL Source (CM) PER MM RANGE WEIGHT 
i* 0.100 12.65 10 
2 0.150 12.45 10 
3 0.100 14.40 3 
4 0.150 13.38 8 


Weighted average result : 12.95+0.2 cm Hg at 15°C is 
equivalent to 1-mm range at 76 cm Hg and 15°C 


* Fig. 2, Curves A, C. 


Hg/mm range whence the true difference jp 
range of the two groups is 1.83+0.04 mm in 
standard air. From the 1938 Cornell range-energy 
curve based, presumably, upon the work of 
Holloway and Livingston’ the slope at 1.6 Mey 
is 0.190 Mev/mm. If we estimate a 2 percent un- 
certainty in this value we get the energy separa- 
tion to be 0.347+0.009 Mev, the alpha-particles 
from Li® being more energetic. This gives 1.802 
Mev for the lithium alpha-particles at 370-ky 
bombardment voltage. If a 10-kv uncertainty is 
assumed in both bombarding voltages, and exact 
masses are used in the computation, the Q value 
for the reaction becomes 3.945+0.06 Mev, prac- 
tically the entire error being due to the quoted 
uncertainty of the electrostatic deflection de- 
termination of the beryllium reaction, an uncer- 
tainty which is probably not as large as assumed. 


APPLICATION OF RESULTS 


From the present measurements and certain 
others some information of interest may be 
deduced. 


H*— He’ mass difference 
Consider the reactions: 
Li®+H!—He'+ He*+ 3.95 +0.06, (1) 
D?+ D?->H?+H!+3.98+0.02, (2) 
D*He'+ He*+22.20+0.04. (3) 
These three energies were determined with 
absorption cells, (1) and (3) were studied in this 
laboratory, (1) by the author and (3) by Smith." 
Reaction (2) is described by Oliphant, Kempton 


and Rutherford," and corrected by Livingston 
and Bethe.’ From these three we get 


(2D*— — (2H'— | 
= (H*—He*) +22.23 Mev. 


TABLE II. Pressure separation of alpha-particles from 
Li*(p,a)He?® and Be*(p,a)Li®. 


TRIAL AP InN cm Hg WEIGHT 
1 17.35 10 
2 17.85 8 


3 19.5 2 
Weighted average result : AP = 17.75+0.4 cm Hg at 15°C 


10 N. M. Smith, Jr., Phys. Rev. 56, 548 (1939). 
uM. L. E. Oliphant, A. E. Kempton, and Lord Ruther- 
ford, Proc. Roy. Soc. 149, 406 (1935). 
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The quantities on the left in parentheses are well- 
known mass-spectrograph brackets. If we use 
Bainbridge’s values we get H*— He*?=0.19+0.09 
while Aston’s predict H*— He*=0.10+0.1. That 
js we may expect H? to emit a low energy beta- 
ray. Neuert’s value of 3.72 for reaction (1) would 
make H® stable by 50 and 130 kev for the two 
cases above. In the second case H* would be 
stable within the experimental error. The radio- 
activity of H* was recently observed by Alvarez 
and Cornog? who estimated the energy of the 
beta-ray to be about 10 Kev. Still more recently 
O’Neal and Goldhaber” have given evidence for 
the reaction Be’+D*—Be*+H' and detected 
beta-rays of about the same energy. 

The present measurements give greater con- 
sistency than Neuert’s with another series of 
reactions, suggested by Rumbaugh, Roberts, and 
Hafstad."* These reactions are 


D*?+ D?-H'+H?+Q,, 
D?+ He®+Qz, 
Li®+n'—He'+ H*+Q;, 
H'—He'+ He? +Q,. 


From these we get 


(Q2+Qs) — (Qi+Qs) =f(Q) =0. 


Table III gives various values of the discrepancy 
f(Q) for different determinations of the Q’s. It is 
seen that the value Q,=3.95+0.06 given by the 
author is in better agreement with all the differ- 
ent combinations chosen than that of Neuert. 
The smallest discrepancy (0.03 Mev) is obtained 
with Bonner’s latest value of the energy of the 
neutrons from the D on D reaction, and with 
Livingston and Hoffman’s earliest value of the 
Li® plus neutron reaction. 


Mass of He* 


In the absence of a good measurement of the 
upper limit of the 8-ray spectrum from H®, the 


nen. D. O'Neal and M. Goldhaber, Phys. Rev. 57, 1086 
4L.H. Rumbaugh, R. B. Roberts, and L. R. Hafstad, 
Phys. Rev. 54, 657 (1938). 


TABLE III. Values of the discrepancy f(Q)=(Q:+Qs) 
—(Qit+Qs) for various combinations of the values 2 Q from 
the reactions D(d,p)H*, Qi; D(d,n)He*, Q2; Li®(n,a)H’, 
Qs; Li®(p,a)He’, Qs. 


ERROR 

REACTION Q QUOTED COMBINATION EV 
3.98 0.02 —0.08 
3.18 0.13 (Q10203'Q4’) +0.11 
3.29 0.08 (Q:0203"0,’) +0.22 
TOs 4.67 0.05 (Qi102'Q30,') +0.03 
3’ 4.86 0.04 ( +0.22 
4.97 +0.33 
3.72 0.08 (Q:1020304 +0.15 
3.95 0.06 (Q:0203'Q4) +0.34 
(Q10203"04) +0.45 

+0.26 

+0.45 

+0.56 


t by and Bethe.® 
Qi, M. L. E. Oliphant, A. E. Kempton, and Lord Rutherford, Proc. 
Roy. Soc. 149, 406 (1935). 
4% 4 2 Bonner and W. M. Brubaker, Phys. Rev. 49, 19 (1936). 
?, T. W. Bonner, Phys. Rev. 53, 711 (1938). 
3, M. S. Livingston and J. G. Hoffman, Phys. Rev. 50, 401 (1936). 
Sis. M. S. Livingston and J. G. Hoffman, Phys. Rev. 53, 227 (1938). 
3s’, L. H. Rumbaugh, R. B. Roberts and L. R. Hafstad, Phys. Rev. 
54, 657 (1938). 
8s H. Neuert, Physik. Zeits. 36, 629 (1935). 
«’, Present measurements. 


best method of computing the mass of He* ap- 
pears to be the following: 


Li‘+H'—He'+ He*+3.95+0.06 Mev, 
Li®+ D*—He‘+ He*+ 22.20+0.04 Mev, 
H'= 1.00812+0.000, 02, 
D?= 2.01472+0.000, 02, 
Het = 4.00388 +0.000, 07. 


From this we obtain He*=3.01688+0.000, 11 
mass units. This is to be compared with the value 
3.01701+0.000, 12 quoted in the paper by 
Barkas."* 

The 1.9-Mev excited state in He* reported by 
Bonner! would result in particles of too short 
range to be detected in this experiment. 

The author wishes to express his indebtedness 
to Professor S. K. Allison for his extremely 
helpful direction of this research; also to Dr. L. 
S. Skaggs for the considerable time he spent 
aiding in the performance of the experiment. The 
author is also grateful to Miss M. R. Jones of the 
Chemistry Department for her aid in preparing 
a strong polonium source. 


4 W. H. Barkas, Phys. Rev. 55, 691 (1939). 
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AUGUST 1, 1940 


The radiofrequency spectra of H2, D2, and HD molecules 
have been obtained for transitions corresponding to reori- 
entations of the rotational angular momenta of the mole- 
cules. From these spectra the rotational magnetic moments 
of Hz, HD, and Dz molecules in the first rotational states 
are found to be +0.8787+0.0070, +0.6601+0.0050, and 
+0.4406+0.0030 nuclear magneton, respectively, all of 
the moments being positive. The relative values are good 
to 0.2 percent and to this precision they are in the ratio 
4:3 : 2 as expected theoretically. The rotational magnetic 
moment of H; in the second rotational state was found to 
be just twice that in the first. From the rotational moment 
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values, the high frequency contribution to the diamagnetic 
susceptibility of Hz is deduced to be (0.093 +0.007) x 10- 
per mole. Certain molecular interaction constants, includ- 
ing the interaction of the quadrupole moment of the 
deuteron with the inhomogeneous electric field of the 
molecule, are also obtained. The results are in complete 
agreement with similar measurements from the nuclear 
reorientation spectrum. Furthermore, the dependence upon 
orientation of the diamagnetic susceptibility of the hydro- 
gen molecule in the first rotational state can be found. It 
is such that — (4.542.0) X10™ per molecule. 


INTRODUCTION 


HE molecular beam magnetic resonance 
method! was originally designed for the 
measurement of nuclear magnetic moments, but 
since rotational magnetic moments of molecules 
are of the same order of magnitude the method 
may be directly applied to the measurement of 
these moments as well. The chief difference is 
that a different ratio of frequency to magnetic 
field must be used so that resonance minima 
corresponding to reorientations of the rotational 
angular momenta instead of the nuclear spins 
may be observed. In addition to the rotational 
magnetic moments of the molecules, certain 
molecular interaction constants can be evaluated 
from the present experiments and compared with 
values of the same constants as found by Kellogg, 
Rabi, Ramsey, and Zacharias.'*? These constants 
measure the interaction energies of the nuclear 
magnetic moments with the rotational magnetic 
field of the molecule, the mutual interaction 
energy of the two nuclear magnetic moments 
with each other, and—most important of all— 
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the interaction energy of the deuteron quadru- 
pole moment?‘ with the inhomogeneous electric 
field of the molecule. A still different constant 
measuring the dependence of the diamagnetic 
susceptibility of a hydrogen molecule upon its 
orientation can also be found. 

The only previous measurements of rotational 
magnetic moments are those of Estermann and 
Stern,’ who found by measurements on para-H, 
molecules in rotational state J=2 that the 
rotational magnetic moment was between 0.8 
and 0.9 nuclear magneton per rotational quantum 
number. The present experiment not only pro- 


“a 
Fic. 1, The radiofrequency spectrum of H: corresponding 


TEMPERATURE 

to transitions in which the rotational angular momentum 


MAGNETIC FIELD iN GAUSS 
changes its orientation. The oscillator frequency is 2.4198 
aii and the oscillating field about 10 gauss (Jr =3.5 
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vides a much more accurate measurement of the 
rotational magnetic moment of H: by an entirely 
different method, but is also applied to D2 and 
HD molecules and can be used to determine the 
signs as well as the magnitudes of the moments. 
Furthermore, the rotational magnetic moments 
of H: both in the first and in the second rotational 


- states can be separately measured, whence the 


proportionality of the magnetic moment to the 
rotational angular momentum may be experi- 
mentally tested. This proportionality was as- 
sumed by Estermann and Stern® in their determi- 
nation of the magnetic moment of the proton 
and is predicted theoretically by Wick,*’ but up 
to now no experimental evidence has been 


obtainable. 


METHOD AND APPARATUS 


The method used in this experiment is that 
described by RMKZ* and KRRZ.'? As the 
latter authors have pointed out, the molecular 
beam magnetic resonance method may be re- 
garded as a means of obtaining the absorption 
and stimulated emission spectra of molecules at 
radiofrequencies. Thus, the minima of beam 
intensity in the accompanying curves occur when 
the oscillator frequency fo is equal to the Bohr 
frequency, »=AE/h, of the molecule, where AE 
is the energy difference corresponding to any 
two orientation states of the molecule between 
which the transitions are allowed, that is, states 


08 


Fic. 2. The radiofrequency spectrum of D, corresponding 
to transitions in which the rotational angular momentum 
changes its orientation. The oscillator frequency is 1.2410 
megacycles and the oscillating field about 14 gauss (Jr = 4.8 
amp.). The temperature is that of liquid nitrogen. 


for which the total magnetic quantum numbers 
differ by an amount AM= +1. Hence by obser- 
vations of the positions of these resonance 
minima the dependence of the energy of the 
molecule on its orientation may be deduced. 
When, as in the present experiment, the re- 
orientation transitions occur in the presence of a 
homogeneous magnetic field, the energy depends 


6G. C. Wick, Zeits. f. Physik 85, 25 (1933). 
7G. C. Wick, Nuovo Cimento 10, 118 (1933). 
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upon the values of the nuclear and rotational 
magnetic moments of the molecule as well as 
upon certain interaction constants within the 
molecule; so, from the observed radiofrequency 
spectrum, the magnetic moments and the inter- 
action constants may be deduced. 

A group of resonance minima may be obtained 
corresponding to reorientations either of the 
nuclear or of the rotational angular momentum, 
i.e., to a change either of m; or m,;. Experiments 


1 ' 
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Fic, 3. The radiofrequency spectrum of HD corre- 
sponding to transitions in which the rotational angular 
momentum changes its orientation. The oscillator fre- 
quency is 1.850 megacycles and the oscillating field about 
12 gauss (Ir=4 amp.). The temperature is that of liquid 


nitrogen. 


of the former type have been done by KRRZ*? 
in measuring nuclear moments, whereas the 
present experiments are limited to changes of my 
from which the rotational magnetic moments are 
obtainable. 

The molecular beam apparatus used for this 
work is that which has already been fully de- 
scribed by KRRZ.' Since the present observa- 
tions were taken with the same apparatus that 
was used in the proton and deuteron nuclear 
magnetic moment measurements the same mag- 
net calibration could be used. As a check that 
the magnet had not changed, measurements of 
the nuclear magnetic moments were made during 
the course of the rotational measurements and 
were found to agree with the earlier measure- 
ments. In this way values of the rotational 
magnetic moments relative to a value of the 
proton moment of 2.785 nuclear magnetons were 
obtained to a greater precision than the precision 
of the absolute calibration of the magnetic field. 


RESULTS 


The experimental results on the groups of 
resonance minima corresponding to Am,;=+1 
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TABLE 1. Magnetic fields at which Hz rotational resonances should occur for a fixed oscillator frequency fo lay a 


=(up/r*)w, and The numerical evaluations of the small s 


terms come from 


KRRZ’s values of H’ 

MAGNETIC FIELDS NUMERICAL 
ans First ORDER SECOND ORDER 


are shown in Figs. 1, 2, and 3 for Hz, De, and 
HD. Approximate values of the rotational mag- 
netic moments (good to about 1.5 percent) can 
be obtained by merely using the values of the 
oscillator frequencies and of the magnetic fields 
at the centers of the groups of minima in the 
formula for the Larmor precession frequency 
(v=ynH/Jh). However, because of the asymmetry 
of the curves, accurate determinations of the 
rotational magnetic moments as well as determi- 
nations of the interaction constants of the 
molecules depend on the detailed nature of the 
interaction functions as is discussed in subse- 
quent paragraphs. 

That the resonance minima shown in the 
accompanying diagrams correspond to changes of 
the rotational my, is indicated by the following 
results: (1) that unlike the nuclear cases these 
groups of resonance minima do not correspond to 
magnetic moments common to either Hz and HD 
or D2 and HD, (2) that a new minimum arises in 
the center of the group of H2 minima when the 
temperature is raised so that a higher rotational 
state is excited, (3) that the magnetic moments 
observed with Hz, HD, and D: are approximately 
inversely proportional to the moments of inertia 
of the molecules as expected theoretically for 
rotational magnetic moments, and (4) that this 
assumption makes possible a correct theoretical 
prediction of the separations of the various 
minima in a single group. 


It is easily possible to account for the number 
of observed resonance minima in Figs. 1, 2, 
and 3. All but one of the curves are taken with 
the source at liquid nitrogen temperature at 
which, because of the low moments of inertia of 
hydrogen molecules, most of the molecules are 
in the zeroth or first rotational state.' Since it is 
impossible to change the value of m, when the 
rotational quantum number, J, equals zero, the 
only molecules of interest are those for which J 
equals one. In the cases of Hz and Dz, where the 
two nuclei are identical, the quaantum-mechanical 
symmetry conditions require that the only 
varieties of molecules that can have J=1 are 
ortho-H; and para-Dz, for both of which the total 
nuclear spin, J, equals one. Therefore with H, 
and D; at liquid nitrogen temperature the only 
molecules concerned in the experiment are those 
with J=1 and J=1, whence there are nine energy 
levels corresponding to all allowed orientations of 
I and J and between these energy levels six 
different transitions corresponding to Am,;=+1 
can occur. This is the reason that six resonance 
minima are observed for Hz and Dz at liquid 
nitrogen temperature. The additional H_ mini- 
mum which appears when dry ice is placed in 
the source trap is due to molecules of para-H, 
with J=0 which are partially excited to rota- 
tional state two at the higher temperature. 
Since J=0, there is no spin-spin or spin-rotational 
interaction so only a single new peak is produced. 
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TaBLE II. Magnetic fields at which Dz rotational resonances should occur for a fixed oscillator frequency fo=ur—D»Ho/h. 
Sp=1/5(H" +H"), =(up/1)m, H’” = —Seq0/4up, a= and k= (E4:— fo) /2up-py. The are as in 
The numerical evaluations of the small second- and third-order terms come from KRRZ’'s values of H’ and Sp. 


MAGNETIC FIELDS 
NUMERICAL VALUES 
First SECOND THIRD or SECOND AND 
Amy mr ORDER ORDER ORDER THIRD ORDERS 
O<1 1 48,050 
a a 
6 2 8,050 14,140,000 
O<>1-1 Hot 441,300 , 7,070,000 
a a a a a H mm 
1 4 41,300 _ 7,070,000 
a a a a H i? 
-1<0 0 4.8,050 _ 14,140,000 
a a a H 
—-l<e0-1 8,050 
a a H 


For HD where the nuclei are not identical, there 
is no symmetry restriction so the number of 
states with J=1 is eighteen, corresponding to 
all possible orientations of the spins of the two 
nuclei and of the rotational angular momentum. 
Between these eighteen states twelve different 
transitions with Amy=+1 can occur, which 
correspond to the twelve observed resonance 
minima. 

By making the same assumptions as were 
made by KRRZ,'? that is, by assuming the only 
interaction energies within the molecule to be a 
spin-spin magnetic interaction energy, a spin- 
rotational magnetic interaction energy, and the 
interaction of a deuteron electrical quadrupole 
moment with the inhomogeneous electric field of 
the molecule, it is possible to account not only 
for their results but also quite well for the 
present ones. Such agreement should be expected 
since the present experiments concern merely 
different transitions between the same set of 
energy levels. However, if in addition to their 
assumption, one takes into account the dia- 
magnetic interaction energies of the molecules 
with the external fieid, the agreement between 
theory and experiment is even better, being 
perfect to well within the experimental error. 
This new assumption amounts to adding a term 
—}tm,H? to the Hamiltonian where tm, is the 
diamagnetic susceptibility of the molecule in the 


rotational state J when the rotational magnetic 
quantum number is m,. This term gives rise to a 
displacement of the rotational resonance minima 
because of the dependence of the diamagnetic 
susceptibility upon &m, clearly being 
least in magnitude when the line of the nuclei 
and hence the electron distributions are most 
nearly parallel to the field (m;=0) and greatest 
when they are most nearly perpendicular to the 
field (m;=+J). This added term of course does 
not affect the agreement already obtained by 
KRRZ,'? since my, and hence &m, is unaltered in 
these transitions whence its effect cancels out. 
With this assumption added to those of KRRZ 
the Hamiltonian of a diatomic molecule in a 
magnetic field H becomes (apart from terms 
independent of orientation) 


§ = J) 


1 3( ui: . 
2 


where wu: and w2 are the nuclear magnetic mo- 
ments, wr the rotational magnetic moment, 
J the rotational angular momentum in units of 
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TABLE III. M at which HD rotational resonances should occur for a fixed oscillator frequency fo=u qolie/h. 
ap =upR- umerical eval 


S= Sr*)av, 


/up, ap and k = (£41— Then 


tions of the small anal and third-order terms come from KRRZ's values of Hp’, Hp’, S, and 


MAGNETIC FIELDS 
0 +S n+ 32070 1220 _ 590,000 
ap 
—le0 3 0 ~Sz- 


h/2n, 1, the nuclear angular momentum in units 
of h/2m, r the radius vector between the two 
nuclei, Q; the quadrupole moment of the kth 
nucleus, g; the gradient of the electric field at 
the position of the kth nucleus divided by the 
charge of an electron, H,’ and H;’ the spin rota- 
tional interaction constants and essentially 
measures of the rotational magnetic fields at the 
positions of the nuclei, and &m, the diamagnetic 
susceptibility of the molecule in rotational state 
J » My. 

From their Hamiltonian and the ordinary 
perturbation theory, KRRZ*? have obtained 
expressions for the energies of the different 
molecular states when J=1 to the approximation 
of the third order perturbation theory. The 
energy expressions needed here are merely those 
of KRRZ with the above diamagnetic term 
added. These energy expressions may be differ- 


enced subject to the restriction Am,=—1 and 
from these differences and from the relation 
hfy=AE, expressions Yor the predicted fields at 
which resonance minima should occur may be 
obtained. These expressions for Hz, D2, and HD 
molecules in the first rotational state are given 
in Tables I, II, and III, respectively, the notation 
being that of KRRZ'? whenever possible. The 
numerical expressions for the small second- and 
third-order perturbations are based on the values 
of the interaction constants obtained by KRRZ. 
It should be noted that all the first-order correc- 
tions are symmetrical about Hp» so the second- 
order corrections are important in that they 
account for the observed asymmetry of the 
results. 

From the expressions of the tables and from 
the values of the interaction constants found by 
KRRZ"? in the nuclear resonance experiments, 
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the theoretical positions of the resonance minima 
may readily be predicted subject to the assump- 
tion of appropriate values for the ur’s and 
£o, the only remaining adjustable constants. 
When this is done, excellent agreement with the 
observations is obtained. This agreement makes 
possible an identification of the quantum num- 
bers associated with each resonance minimum 
and shows the complete consistency of these 
results with those of KRRZ. The quantum 
numbers to be identified with each minimum are 
given in Tables IV, V, and VI, the correlation 
between these tables and the curves being the 
designation letters Az, etc. The best determina- 
tions of the experimental positions of the reso- 
nances are also given in these tables. 

After this identification has once been made, 
the process may be performed in the reverse 
direction and the experimental results used to 
evaluate the various interaction constants. The 
values so found not only provide a check upon 
the values of the same quantities as obtained 
from the nuclear resonance minima, but in some 
cases, like Hz where the separations in gauss of 
the rotational minima are about seven times 
greater than the separations of the nuclear 
minima, the determinations can be made more 
accurately from the rotational than from the 
nuclear experiments. This evaluation has been 
performed by equating the theoretical resonance 
positions of Tables I, II, and III to the corre- 
sponding observed positions listed in Tables IV, 
V, and VI and by solving for the desired con- 
stants by the method of least squares. The values 
of the various quantities obtained in this way 
are as follows: 


fo=2.420 megacycles: gauss, H’ 
= 26.8+0.3 gauss, H’ =33.9+0.3 gauss, and 


—(3.942.5) X10-! 


per molecule ; 

(b) For Dg in the first rotational state when 
fo=1.2410 megacycles: Hy»=3696.2 gauss, H’ 
=13.7+0.8 gauss, Sp=19.55+0.20 gauss, and 


= —(5.142.5) K10-*! 


per molecule; 

(c) For HD in the first rotational state when 
fo=1.850 megacycles: H)»=3677.4 gauss, Hp’ 
=20.5+1.5 gauss, Hp’=20.23+40.5 gauss, S 
=11.70+0.5 gauss, E=17.30+0.3 gauss, and 


§+1—fo= —(4.643.5) X10-#! 


per molecule. 
The number of Po resonance minima 
exceeds the number of available parameters, but 


TABLE V. Observed positions of Dz rotational resonances 
and comparison with theory when fo= 1.2410 megacycles. The 
theoretical positions are based on the values of Ho, H’, Sp, 
and k evaluated from these data by the method of least squares. 


EXPERI- 

MENTAL 

PosITION 

DESIGNATION OF THEORETICAL 

OF RESONANCE POSITION 
Am, my] MINIMUM (Gauss) (Gauss) 
0<1 1 Ar 3787.3 3786.8 
Oeil 0 CL 3473.6 3473.9 
Br 3849.4 3849.5 
—-le0 1 Bi 3565.1 3565.2 
—le0 0 Cr 3922.6 3922.8 
AL 3610.5 3610.0 


TABLE VI. Observed positions of HD rotational resonances 
and comparison with theory when fo= 1.850 me, The 
theoretical positions are based on the values of H o, Hp’, Ho’, 


(a) For Hg in the first rotational state when 5, E, and k evaluated from these data by the method of least 
squares. 
TABLE IV. Observed positions of Hz rotational resonances 
and comparison with theory when fo=2.420 megacycles. The EXPERIMENTAL 
theoretical positions are based on the values of Ho, H', H", DESIGNA- 
and k evaluated from these data by the method of least squares. dm, tp mp i ia A rime 
1 Di 3556.5 3557.1 
Soren O<1 0 Ar 3737.7 3736.1 
PosITION Oeil -1 Ex 3499.6 3499.4 
DESIGNATION OF —-} 1 By, 3631.8 3629.2 
RESONAN 
1 AL 3573.0 3572.9 1 Ci 3583.2 3584.9 
0<1 0 Bi, 3354.3 3354.7 —l<0 ' 0 Fi 3469.6 3467.9 
-1 Cr 3920.3 3919.6 —le0 Br 3750.0 3750.9 
-l<e0 1 CL 3319.3 3319.7 -le0 —-} 1 Er 3857.1 3855.6 
-le0 0 Br 3871.3 3871.5 —-le0 —-} 0 AL 3631.8 3636.4 
-le0 -1 Ar 3653.6 3653.5 —-le0 -} -1 Dr 3806.5 3806.3 
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the consistency of the results is shown by a 
comparison of the experimental positions of the 
resonance minima with the theoretical positions 
predicted with the above values of the interaction 
constants. This comparison is given both on the 
graphs and in the last two columns of Tables IV, 
V, and VI. It will be seen that the agreement 
between theory and experiment is excellent— 
good to about one part in eight thousand for 
H, and D». Although the agreement in the case 
of HD is worse, it is still very good (about one 
part in two thousand). That the agreement 
should be worse in this case is not surprising 
since the experimental observations were by far 
the most difficult for HD. This was due to the 
fact that the rotational moment being reori- 
entated was less than one-fourth the magnitude 
of the proton moment which occurred in the 
same molecule. Hence a sufficiently strong in- 
homogeneous field to throw out reoriented mole- 
cules caused all the molecules to make such 
large excursions because of the proton moment 
that many nonreoriented molecules were lost, 
which made the observable intensity of the 
resonance minima very weak. In fact, the 
average depth of an HD rotational resonance 
minimum of Fig. 3 was less than one percent of 
the total beam intensity. As this refocusing 
difficulty occurred for all the rotational measure- 
ments in a greater or lesser degree, here, unlike 
the nuclear case, no fundamental significance can 
be attached to the measured intensities of the 
resonances. 

It will be noted that whereas most of the 
above quantities have been measured to a high 
precision, such is not true of §+;—&». This is due 
to the fact that at the magnetic fields of the 
experiments this term only slightly affects the 
positions of the resonance minima (1.8 gauss out 
of 3600 gauss at the most). At the time the experi- 
ments were performed it was not realized that 
such a quantity as §4:—& could be obtained 
from the results, so there was no special effort to 
make the observations under the most favorable 
conditions. It is hoped that at some future date 
the experiments will be carried out at a much 
higher frequency and field where this diamag- 
netic effect, which varies as H?, will be relatively 
much stronger. 

From the frequencies, from the least-squares 
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values of Ho listed above, and from the relation 
br=Jfoh/Ho (2) 


the values of the rotational magnetic moments 
of the molecules are obtained. These values are 
= 0.8787+0.0070, ur_pv, =0.4406+0.0030, 
and nuclear magneton. 
The moment values here quoted are relative to 
a value of the nuclear moment of the proton of 
2.785 nuclear magnetons.' The relative values 
should be accurate to 0.2 percent although the 
absolute values of the moments, like the absolute 
value of the proton moment, are subject to an 
error of possibly 0.7 percent due to possible errors 
in the absolute calibration of the magnetic 
field.! Except for the evaluation of the effect of 
the second-order perturbations which slightly 
shift the centers of gravity of the groups of 
resonance minima, the determinations of the 
rotational magnetic moments are completely 
independent of the above detailed theory of the 
molecular interactions. So even if this theory 
were completely wrong, which is very highly 
improbable in view of the extensiveness of its 
experimental verification, the rotational mag- 
netic moments should be correctly evaluated to 
an accuracy of at least 1.5 percent. Although the 
most accurate evaluation of the rotational mag- 
netic moments comes as above from the observed 
value of Ho/fo, an independent check of the 
values can be obtained from the separations of 
the different minima in the group. From Tables 
I, II, and III it will be noted that the separations 
of the minima depend on the ratios a=ypr/yp, 
etc. Therefore the agreement of the values of 
interaction constants obtained in this experiment 
with those obtained in the nuclear experiments 
indicates that the ywr’s have been correctly 


chosen. 
It should be noted that, in the theoretical 


prediction of the resonance minima positions 


given above, it was implicitly assumed that the 
ur’s were positive. Had the yr’s been negative 
there would have been no change in the predicted 
positions as far as the first-order perturbations 
are concerned since by Tables I, II, and III the 
first-order energy perturbations are symmetrical 
about H». However, the direction of the shift due 
to the second-order perturbation and hence the 
nature of the asymmetry would have been 
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TaBLe VII. Summary of results. The notation is in agreement with Table XII of KRRZ* The quantity in parenthesis in 
the wp/r* column ts the value obtained by multiplying up/r* by up/up. It is used in obtaining Se*qQ/4up from the D: observations. 
The first quantity for H’ in HD is its value at the position of the proton and the second 1s at the position of the deuteron. 


(up/P)m Constant =H’ SeqQ/4up XHF (10-8 PER 

MOLECULE (N. M.) (Gauss) (Gauss) (Gauss) (Gauss) (107-* PER MOLE) MOLECULE) 

H. when J=1 0.8787+0.0070 33.9 +0.3 26.8 +0.3 0.094+0.007 —3.942.5 
Hz when J=2 1.757+0.014 or 

2(0.8786+0.0070) 
when J=1 0.4406 +0.0030 (10.42) 13.7 40.8 87.36+41.2 0.09340.007 —5.1+2.5 
HD when J=1 0.6601+0.0050 34.2141.5 10.26+0.40 86.50+1.5 0.09340.007 —4.643.5 
5 +1. 


Average of Hz, De, and HD 


0.093+0.007 —4.5+2.0 


reversed. Therefore, the existence of the slight 
second-order asymmetry of the experimental 
curves makes possible a determination of the 
sign of the rotational magnetic moments. In all 
three cases the rotational magnetic moments are 
found to be positive. 

One of the most significant results of this 
experiment is obtained by a comparison of the 
H, resonance minima with the source at liquid 
nitrogen and at dry ice temperatures, respec- 
tively, as shown in Fig. 1. At dry ice tempera- 
tures, in addition to the six resonance minima 
found for Hz at the lower temperature, where 
the only nonzero rotational state is the state 
one, a seventh minimum appears midway be- 
tween the two central minima of rotational 
state one. This minimum may be assumed to be 
due to the second rotational state of H» which, 
though unexcited at liquid nitrogen temperature, 
is partially excited at dry ice temperature. 

Since only para-Hz may be in the second rota- 
tional state and since para-H: has a zero resultant 
nuclear spin, there is no spin-spin or spin-rota- 
tional interaction so only a single unperturbed 
resonance minimum is obtained from which the 
rotational magnetic moment may be directly 
deduced from Larmor’s formula. The magnetic 
moment as thus obtained is 1.757+0.014 
= 2(0.8786+0.0070) nuclear magnetons or just 
twice the rotational moment of rotational state 
one. This result may also be obtained directly by 
observing that the new peak due to the second 
rotational state falls midway between the two 
central peaks of the first rotational state, which 
by Table I are unaffected by the second-order 
perturbations. Therefore, ur/J must be the same 
for J=1 and J=2 to within the experimental 
error of 0.2 percent in the determination of ratios 


of magnetic moments; in other words, as far as 
the first two rotational states of Hz are concerned 
the rotational magnetic moment is proportional 
to the rotational angular momentum. 

As a check upon the consistency of the results, 
observations were made upon Hg at half the 
frequency and consequently at half the value of 
the magnetic field used in the published curves. 
The results are completely consistent with those 
at the higher frequency. Also the rotational 
magnetic moments of Hz and Dz molecules were 
measured with an earlier form of the present 
apparatus in which the path of the beam in the 
radiofrequency field was only 2.7 cm long in 
contrast to the present 13.5 cm length. Because 
of the shortness of the field and the smallness of 
the beam intensity, it was necessary to use a 
much stronger oscillating magnetic field whence 
the resolving power of the apparatus was much 
less (the half-widths of the resonance minima are 
proportional to the strength of the oscillating 
magnetic field®). For this reason it was impossible 
to resolve the individual peaks and. there was 
consequently a much greater uncertainty in the 
evaluation of the magnetic moments than now. 
However, the moments could be determined with 
an estimated error of two percent and these 
values were reported.* The earlier measurements 
agree with the present ones to within the former’s 
estimated error. 


DISCUSSION OF RESULTS 


A summary of the results obtained in this 
experiment is given in Table VII. Those quanti- 
ties which have also been measured by KRRZ 
are (up/r*)y, (up/r*)», H’, and 5e’gQ/4up. The 


‘1. I. Rabi, Phys. Rev. 51, 652 (1937). 
*N. F. Ramsey, Jr., Phys. Rev. 55, 595(A) (1939). 
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portion of Table VII pertaining to these quanti- 
ties may be directly compared with Table XII 
of KRRZ.? It will be seen that the agreement is 
excellent, being well within the estimated experi- 
mental error in all cases. This agreement and 
the close agreement between theory and experi- 
ment in Tables IV, V, and VI provide very 
strong evidence for the validity of the theory on 
which this work is based and in particular for the 
quadrupole moment of the deuteron. A full 
discussion of this agreement, a final averaging of 
all the experimental results including those of 
this paper, and a computation of the deuteron 
quadrupole moment from the final average of 
the quantity 5e’gQ/4yp is contained in the paper 
of KRRZ? and will consequently not be discussed 
further here. 

The observed value of the rotational magnetic 
moment of Hz: falls within the limits 0.8 to 0.9 
nuclear magneton per rotational quantum num- 
ber set experimentally by Estermann and 
Stern.’ It also agrees with the much wider 
theoretical limits set by Wick.*? The uncertainty 
of Wick’s result arises from the sensitiveness of 
the computed rotational magnetic moment to the 
nature of the electron wave function of the mole- 
cule that is used, the predicted moment being 
+0.36 nuclear magneton with Wang wave 
functions and +0.93 nuclear magneton with 
another wave function of Wick’s invention. It is 
perhaps worth emphasizing that the observed 
rotational magnetic moment for He: is much 
different from that which would be obtained if 
the electrons were a rigid charge distribution 
which rotated with the same angular velocity as 
the nuclei. Such a naive theory would predict’® a 
rotational magnetic moment in the state J=1 
of —3 nuclear magnetons which is much different 
in magnitude and opposite in sign from the ob- 
served result, the electron contribution to the 
magnetic moment from this theory being over 
thirty times greater than its experimental value. 
The necessity for a theory such as Wick’s in 
which the electron contribution to the rotational 
magnetic moment is much less than on the rigid 
body model is thus evident. 

The result that the rotational magnetic mo- 
ment of Hz: is proportional to its angular 


1 R, Frisch and O, Stern, Zeits. f. Physik 85, 25 (1933). 
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momentum would be completely obvious if the 
electrons were assumed to behave as a rigid 
charge distribution, but, as just pointed out, this 
assumption is necessarily false. However, the 
proportionality of the rotational magnetic mo- 
ment to J is also predicted by Wick’s theory. As 
this prediction is the most important application’? 
of the theory, its direct experimental proof in 
these experiments is particularly useful. 

The experimental values of the rotational 
magnetic moments of Hz, HD, and D; as listed 
in the preceding section are in the ratio 
3.994 : 3.000 : 2.003; that is, to within the ex- 
perimental error of 0.2 percent they are in the 
ratio 4: 3:2 or in the inverse ratio of the re- 
duced masses of the molecules. Since the rota- 


tional angular velocities in the first rotational 


state are inversely proportional to the reduced 
masses the experiment shows that the rotational 
magnetic moments are proportional to the rota- 
tional angular velocities in a comparison of 
He, HD, and Dz molecules. That such should be 
the case would certainly be expected in a com- 
parison of Hz and Dz molecules; however, that 
it should be so in a comparison of Hz with HD is 
not immediately apparent, for not only are the 
angular velocities of the two molecules different 
but so also are the positions of the axes about 
which the molecules rotate, since the center of 
mass of the HD molecule is not midway between 
the two nuclei. Still even for this case if the elec- 
trons behaved as a rigid charge distribution 
attached to the nuclei the rotational magnetic 
moment would be proportional to the rotational 
angular velocity, since for a rigid molecule with 
zero total charge and with central symmetry of 
the charge distribution the rotational magnetic 
moment is invariant to the position of the axis of 
rotation. That the magnetic moment should be 
inversely proportional to the reduced mass for 
unsymmetrical molecules like HD is on the basis 
of Wick’s theory not immediately obvious. In 
fact, to make it so the electron contribution to 
the moment of HD must be more than twice as 
great as its contribution to Hz even though the 
rotational angular velocity is less. This, however, 
may be proved as is shown in the appendix to 
the present paper where the general expression 
for the rotational magnetic moments of He, Ds, 
and HD molecules is derived. 


| 

| 
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Although the 4:3: 2 ratio of the rotational 
magnetic moments holds to within the estimated 
experimental error of 0.2 percent in the determi- 
nation of the ratios of moments, the discrepancy 
is very near the maximum experimental error 
and seems fairly persistent. It is in the direction 
that the moments of the lighter molecules are 
slightly less than would be expected. Such a 
discrepancy could probably be explained as the 
effect of the zero state vibrational motion of the 
molecule upon the mean separation of the nuclei. 
A quantitative calculation of the magnitude of 
the effect was not attempted. 

From the theory of Wick or from the more 
general theory given in the appendix it is possible 
from the experimental measurement of the rota- 
tional magnetic moments to evaluate the high 
frequency term™ of the ordinary diamagnetic 
susceptibility of H2, Dz, and HD molecules. This 
is due to the fact that the summation in Eq. (6) 
of the appendix is just exactly the same as the 
summation which occurs in the theory of the 
diamagnetic susceptibility of molecules.® In fact, 
a comparison of the two shows that (in the nota- 
tion of Van Vleck and the appendix) 


e’Lp* (: 2MiM2 3) 
(Mi+M2)M 

So with the observed values of the rotational 
moments, the high frequency term in the dia- 


magnetic susceptibility of H:, D2, and HD are 
as listed in Table VII, the average value being 


Xmole-HF = 0.093+0.007 X 10-* per mole. 


Xmole—HF 


This result may be subject to a somewhat larger 
error than 0.007 due to possible inaccuracies in 
the perturbation theories on which the calcula- 
tion is based. The result calculated from Wang 
wave functions by Van Vleck and Frank" of 
Xmole-HF = 0.51X10-* thus appears to be about 
five times too big. However, if instead of the 
Wang wave function another H, wave function, 
used first by Wick,*’ is used in the calculation, 
the theoretical value for the high frequency part 
of the susceptibility may be made to drop even 
below the above experimental value. 

Blectré 
(Onford Press, 1933) 


” J. H. Van Vieck and A. Frank, Proc. Nat. Acad. Sci. 
15, 539 (1929), 
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APPENDIX 


Theoretical expression for the rotational magnetic 
moments of He, De, and HD molecules.—Al- 
though Wick’s treatment®’ applies only to a 
symmetrical molecule, it may readily be extended 
to the case of a molecule one of whose nuclei is 
of mass M, and the other of M; with the result 
that the rotational magnetic moment of the 
molecule in rotational state J is, in nuclear 


magnetons, 
My)M 
,|(O| Meo| m) |? 
E,—Eo 


(4) 


where M is the mass of the proton, m the mass 
of the electron, and p the distance between the 
nuclei. (0| M¢o|) is a matrix element for the non- 
rotating molecule of the electronic angular mo- 
mentum about an axis through the center of 
mass of the molecule and perpendicular to the 
internuclear line. Since the angular momentum 
is about the center of mass of the molecule and 
not about a point midway between the two 
nuclei, the quantity in brackets is different for 
HD than H; or Dz. The first term in the above 
expression is the nuclear contribution to the 
rotational magnetic moment and the second is 
the electron contribution, neither of which alone 
is inversely proportional to the reduced mass. 
To make possible a comparison of the H, and HD 
results the electronic angular momenta must be 
expressed about an axis which is the per- 
pendicular bisector of the line joining the two 
nuclei. This transformation of the quantity in 
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the brackets can be performed with the aid of 
the invariance properties of the diamagnetic 
susceptibility shown by Van Vleck." (It should 
be noted that the symmetry of the molecule 
is such that where 
P(0, m) are the total angular momentum matrix 
elements used by Van Vleck.) The result is 


(M.?+ My?) M 
MMM + M2) 

mp? E,— Eo 
1(M\+M2)M 

CX JJ, (5) 


where (0| /;,.|) is the matrix element of electron 
angular momentum about an axis which is the 


perpendicular bisector of the line joining the two 
nuclei and is consequently the same for He, Dz, 
and HD except for second-order effects and where 
(ro?)» and (r.*)y are the mean square distances of 
the electron distribution from the center of mass 
and the midpoint of the molecule, respectively. 
However, if d is the separation of the midpoint 
from the center of mass (ro?) —(r2)4=d?. From 
this and the value of d=}(M2—M,)/(Mi+M.) 
the above equation reduces to 


Ms 
|(O| M¢.| n) |? 
i-—— 


whence ur should be inversely proportional to 
the reduced mass as found experimentally. 
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The spinor equations for arbitrary spin and rest mass zero are examined in some detail. Fierz 
has shown that for all values of the spin f(>4) there exists only two ‘“‘really”’ independent 
plane wave solutions instead of (2f+1) when the rest mass is not zero. Fierz later showed, in 
rather a complicated way, that these two plane waves correspond to components of spin +f 
along the momentum vector of the wave. We will arrive here at the same result but in a much 


simpler and more direct -way. 


INTRODUCTION 


FIERZ! has given a detailed discussion of 
e the spinor wave equations for particles 
with arbitrary spin proposed by Dirac.’ Fierz 
showed that when the rest mass is zero there is 
a certain degeneracy of the equations, and that 
for a given spin f there are only two really 
independent plane wave solutions instead of 
(2f+1). The interpretation of this was not given 
in his first paper. In a subsequent paper,’ how- 
* 
1M. Fierz, Helv. Phys. Acta 12, 3 (1939). 


?P. A. M. Dirac, Proc. Roy. Soc. A155, 47 (1936). 
3M. Fierz, Helv. Phys. Acta 13, 45 (1 940). 


ever, he discusses the rest mass zero case in more 


detail, and arrives at the result that the two 


plane waves correspond to axial spins +f about 
the momentum vector of the plane waves. His 
work is rather complicated. Here we will arrive 
at the same result in a much simpler way. For 
coherence it is necessary to give here some of the 
work appearing in Fierz’s paper. The first section 
will be largely a restatement of Fierz’s results, 
although the notation I will use will be slightly 
different. The second section will contain a 
detailed study of the rest mass zero case. In the 
present discussion no distinction will be made 
between the cases of integral and half-odd spin. 


as 
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Such a distinction is quite unnecessary, since the 
work is exactly the same for both. The only 
essential difference between the two cases is the 
form of the energy momentum tensor and charge 
current vector, but these do not enter in the 
present work. 

The spinor notation which I shall use is tha 
of Veblen. A very good concise account of this 
formalism has been given by A. H. Taub,‘ to 
which the reader is referred. 

The present work was developed during the 
course of a number of reports on Fierz’s paper 
that I gave to a seminar on spinor theory, con- 
ducted by Professors O. Veblen and J. von 
Neumann. I am much indebted to Professors 
Veblen and von Neumann for helpful discussion. 


1. THE SPINOR EQUATIONS FOR PARTICLES WITH 
ARBITRARY SPIN AND REstT Mass NONZERO 


The equations are stated most simply for our 
purposes in the following form: 


(o, r=1, 2,3, 4), (A, B=1, 2), 


Bi B, 
where ¥(0) 


is a spinor field of the type indicated by its 
indices and is symmetrical in all its dotted and 
undotted indices, g’* is the metric tensor which 
has the components 


(1) 


\|g7"||= 


0 0 -1 


in all Galilean frames, and x= mc/h where m is 
the rest mass of the particle. Furthermore 


By 
¥(0) means 


In addition to (1) the spinor ¥(0) (spinor 
indices are omitted when not necessary) is to 


‘A. H. Taub, Ann. of Math. 40, 937 (1939). For typo- 
— reasons the usual dotted indices will here be 
replaced by indices with dots put in front of the indices. 


satisfy the following condition: 


Bagi 


¥(1) 
defined by 
¥(1) AwW(0);¢ 


is symmetric in all its dotted and undotted indices, 
where the g’“z are the components of the four 
basic spinor matrices. (See Taub, reference 4.) 

For the case of integral spin f=r, we must 
have s=r in the above equations. : 

For the case of half-odd spin f=r+4, we must 
have s=r-+1 in the above equations. 

The above two statements will become clearer 
in the later work. : 

The most important property of the matrices 


is 


cA +B tA B or A 
BO Bg c)=g (3) 


and the complex conjugate of (3). From this it 
follows that 


pAi B B, 


B 
= «*y(0) (4) 


The second-order equation (1) may thus be 
replaced by the first-order equations (2) and (4) 
where the ¥(1) and ¥(0) are spinors symmetric 
in all their dotted and undotted indices. This is 
the usual form in which the equations are given. 

We can now define in general a y(m) 
(m=1, 2, -+-r) by 


Any (m — 1) 50m 


Bayi om Baym 


=g 


AmW(0) ; (5) 


It is easily seen that ¥(m) is symmetric in all 
its dotted and undotted indices in virtue of the 
like symmetry of ¥(1) and ¥(0). Furthermore the 
¥(m) all satisfy equations like (4), viz., 


pAmsi Biss 


g p 


We note that in particular 


Baar 


¥(r) 


] 
| 
| 
| 
- 
0 0 | 
| 
—— | 
Beam 
| 
| 
| 
i 


has only dotted indices and is symmetric in all 


these indices. 
In a similar way we can define spinors ¥(—m) 


(m=1, 2, -++s) by 


Bs omArym Ai *Aramei 


¥(—m) =9 —[m—1]):00 
omAram 

= 9 By: ‘9 Bu (0) ; 
1 ai Arai omAram 

=—9 By ‘9 Bud Bayi 


Bagi 
XV(1) ; (7) 


using (4) to get the last part of the equation. 
Again from the symmetry of ¥(1) and (0) in all 
their dotted and undotted indices, it follows from 
(7) that ¥(—m) has like symmetry. 

In particular 


AiA2**Arye 


¥(—s) 


is a spinor with undotted indices only, and is 
symmetric in all of them. The ¥(—™m) satisfy the 
following equations (analogous to (6)) in virtue 
of (1): 

p Bn Bs 


ArsmW(—™M) ; 
**Aram--1 Bm** 


“Bs 
1) (8) 


We can write the above equations (dropping 
the spinor indices) as follows: 


and so on up to 
and (9) 
up to 


The notation will be clear if we note that the 
first two equations of (9) are Eqs. (2) and (4). 
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From any of the above pairs of equations, 
together with the symmetry in dotted and un- 
dotted indices of the two y’s involved, can be 
deduced all the other equations as well as the 
symmetry of the rest of the y’s in dotted and 
undotted indices. From this it appears that we 
could have started with Eq. (1) involving a y 
with any number of dotted and undotted indices, 
rather than the ¥(0) we started with, having 
either r=s or r=s+1. One reason why the ¥(0) 
has a preferred position is that the physical 
quantities like the energy momentum tensor 
and charge current vector are most simply ex- 
pressed in terms of ¥(0). As we will see later, 
when the rest mass is zero the y that we start 
with definitely has a preferred position. If we 
start with any y other than ¥(0), the energy 
momentum tensor and charge current vector 


formed from it is identically zero. This follows 


from the fact that these tensors defined in terms 
of y’s other than ¥(0) involve second and higher 
derivatives, i.e., terms like 


Pi orp 


say, and these are zero in virtue of (1) with «=0. 
Let us now consider plane wave solutions of 


the equations. 


¥(0) = ¥(0) e (10) 


will be a solution if 


— 2° (11) 


and the components of ¥(0) are constants. 
In addition (1) given by 


o Bayt Be 


Bagi 


W(1) = ikeg 


must be symmetric in all its dotted and undotted 
indices. (0) has (r+1)-(s+1) -independent 
components, considering only its symmetry in 
its spinor indices. The symmetry of the right- 
hand side of (12) imposes r-s conditions on the 
components of (0), however, so that it has 
(r+1)(s+1)—rs=r+s+1 independent com- 
ponents. 

The same reasoning applies to all the Y(m) 
(m=r, r—1, -+-1, 0, —1, —2, —s). Each 
has r+s+1 independent components, and any 


Qe 


i ( 
| 
| 
| 
| 
| 
| (12) b 
| 0 
n 
re 
tl 


v(m) is expressible in terms of 


Bags 
V(r) 


by the following relation: 


¥(r—m) 


Am 
XJ Besr-mar V(r) 


(m=0, 1, 2, aioli r+s), (13) 
A _ @A 
where B=ikeg RB. (14) 
By (13) the components of any W(m) are 
expressed in terms of the (r+s+1) independent 
components of V(r). 
Let us now choose a special coordinate system 
in such a way that 


Ky =Ke=K3=0, tx. 
Then 
44 
g B=ing B. 

If we now apply a space rotation to the system 
(rotation of the x', x*, x*, leaving x‘ invariant) 
the corresponding spinor transformation is a 
unimodular linear transformation leaving g*4z 
numerically invariant, i.e., a unitary transforma- 
tion; since in a special spin coordinate system 
the matrix (obtained from by lowering 
the dotted index) has the components 


Thus under a space rotation for a system with 
zero momentum (k;=Kk2:=x3;=0) the connection 
between the (r+s+1) independent components 
of V(r) and the components of ¥(r—m) is 
numerically invariant. Now the components of 
V(r) transform according to the irreducible 
representation D,,4.) of the rotation group under 
the unitary spin transformation associated with 
a space rotation. Further since the gg, are 
numerically invariant under the unitary spin 
transformation, the (r+s+1) independent com- 
ponents of ¥(r—m)—and in particular ¥(0)— 
also transform according to Dy:-+.).° 

“See B. L. van der Waerden, Gruppentheoretische 


Methode in der Quantenmechanik (Springer, Berlin, 1932), 
Chap. III. 
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Hence the equations correspond to a spin 
3(r+s). From this follows the earlier statement 
that when s=r, the corresponding spin is f=r, 
and when r=s+1, the corresponding spin is 


2. Rest Mass ZERO CASE 


The work of the preceding section has been 
largely a restatement of Fierz’s results. We now 
consider the rest mass zero case. Fierz showed 
that when the rest mass is zero, the ¥(0) had 
only two “really independent’’ components, 
since it was possible to add to ¥(0) a set of plane 
wave solutions of the equations depending on 
(r+s—1) independent constants, without chang- 
ing the energy momentum tensor or charge 
current vector for the system. The transforma- 
tion of the y’s by the addition of these terms he 
calls ‘‘gauge transformations.’’ The supposition 
was that by a suitable choice of the gauge trans- 
formation one could reduce all but two of the 
components of ¥(0) to zero. Here we will inves- 
tigate the nature of the two components of ¥(0) 
which are not zero. Working in a special coor- 
dinate system for which x2=«x;=0 we will show 
that there are two components of ¥(0) which are 
“gauge invariant” and nonzero, i.e., by no 
choice of the gauge transformation can these 
two components be reduced to zero. These two 
components then are the two “really’’ inde- 
pendent components of ¥(0). We will then show 
that under axial rotations about the x' axis 
through an angle @, these two components are 
multiplied by 


so that they correspond to components of spin 
+f=+}(r+s) along the x' axis. But the mo- 
mentum vector lies along the x! axis (since we 
have chosen coordinates in such a way that 
ke=«x3=0). Thus the two gauge invariant states 
correspond to components of spin +f in the 
direction of the momentum vector. Thus we here 
arrive at the result of Fierz’s second paper in a 
simple and direct manner. 

The equation for the rest mass zero case is (1) 
with x=0, together with the condition that ¥(1) 
defined by (2) is symmetric in all its dotted and 
undotted indices. Then (1) satisfies (4) with 
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x=0. The ¥(m) (m=1, 2, ---, r) are defined as 
before by (5) and satisfy (6) with x=0. 

The ¥(—m) (m=1, 2, ---, s) are defined as 
before in terms of ¥(0) by the first part of Eqs. 
(7), but the symmetry of ¥(—m) in dotted and 
undotted indices no longer follows by our 
previous argument, since there the last part of 
(7) (which no longer applies since x = 0) was used. 
It is still true, however, as we will see later, that 
W(—m) is symmetric in dotted and undotted 
indices in virtue of the like symmetry of (0) 
and (1). So that for the rest mass zero case the 
equations (9) are replaced by similar equations 
with 

Let us now consider plane wave solutions of 
the equations as before 


¥(0) = W(o) e (15) 


where now 
2° Kok, =0. (16) 


We then have the following state of affairs: 
W(0) still has (r+s+1) independent components 
—as before a result of the symmetry of ¥(1) and 
W(0) in their dotted and undotted indices. We 
will show that for a special coordinate system all 
the remaining components of (0) are zero. 
Further, in this coordinate system W(m) has 
(r+1—m) nonzero components and in particular 
W(r) has only one nonzero component. Similarly 
W(—m) has (s+1—m) nonzero components and 
in particular ¥(—s) has only one nonzero com- 
ponent. 

In what follows we wish to be quite explicit 
and give the actual matrices ||g*“s|| in a par- 
ticular spin coordinate system. Let the 


ol 


le 


a2 
be given by 
VA 0 1 2:4 1 
= ’ B|| = ’ 
0 0 -1 
(17) 
3A 0 4A 0 
0 7 1 


Then let us choose a special space time frame of 
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reference in such a way that 
x3=0, x= (18) 


Then we have 


Bi) = TK 
0 
(19) 
0 
= (x, = — 
0 


For definiteness let us take x,;=x, and denote 
the matrix 


oA A 
Bi = by lg 


0 


Now ¥(1) is by definition given by 


=1keg Ai ¥(0) (20) 


In the special coordinate system, we have from 
(19) and (20) 


Ar Ar Bie 
=0. (21) 


Now (1) is symmetric in all its dotted indices. 
Hence the only nonzero components of ¥(1) are: 


W(1) 


and there are r such nonzero components.* It 
follows further from (20) and (19) and the 
dotted and undotted index symmetry of (1) 
that the following components of (0) are zero 
in our special coordinate system: 


Ba 
all terms v(0) 

excepting 


(A, B=1,2). (22) 


The number of terms of ¥(0) that are zero is 
thus r-(s+1)—r=rs. Since the total number of 
terms of (0) is (r+1)(s+1) the number of 
nonzero terms is (r+s+1), i.e., just the number 
of independent components of ¥(0). 


* In the following we shall frequently use component and 
term for independent component and independent term. 


q 
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Similarly we have (in the special coordinate 
system) 


Baym 
v(m) 
‘Basi 
=g . 


“Baim 


Bi 
AnW (0) (23) 


and it follows from the form of ||g“ || and the 
symmetry of ¥(m) in all dotted and undotted 
indices, that the only nonzero terms of ¥(m) are 


v(m) (r+1—m) in number. 


In particular the only nonvanishing term of 
V(r) is 


V(r) ¥(0) (24) 


Similar results will hold for the ¥(—m) pro- 
vided ¥(—1) is symmetric in all dotted and un- 
dotted indices. That this is so is easily verified 
in the special coordinate system. For 


¥(—1) =§ 


Now from (25) 


Arsi Bs 


B, (25) 


Ayt**Ar-1 12 Bs 
v(-1) 


Bs 


=0 from (22) and 


21 Bs 


v(-—1) =0 from (25) and (19). 


So the symmetry of ¥(—1) in dotted and 
undotted indices follows in a special coordinate 
system, and hence in any spin coordinate system. 

The rest of the argument then proceeds as 
before. ¥(—m) has (s+1—m) nonzero compon- 
ents and in particular ¥(—s) has only one non- 
zero component given by 


W(—s) ¥(0) 
Let us now consider the effect of the ‘‘gauge 
transformations”’ of Fierz on the ¥(0), ¥(m) and 
v(—m). 

Fierz showed that the energy momentum 
tensor and charge current vector for the system 
was unchanged if one replaced the plane waves 


o 


by 


where the ¥’(0) is given by 


v'(0) 


(0) 


@BiAs 
ting 


+ikeg 


AiA3***Ar Bi 


+ikeg 


+1keg (0) 


where 
A B, 
(0) 


is symmetric in all its dotted and undotted 
indices and in addition (1) defined by 


Bayi Ay B, 


(1) 


Bayt 
=tkeg (28) 
is symmetric in all its dotted and undotted 
indices. 

We must verify that ¥’(0) and W’(1) are sym- 
metric in all their dotted and undotted indices. 
For ¥’(0) this is apparent from the definition 


(27). W’(1) is given by 


w'(1) A, W'(0) 


=ikeg (1) 


Bayi 


= 


+ikeg (1) 

@B:Az By Bo Bagi 
(1) 


@BrAs Bi Bayt 


In deriving the above we use (3) and (16). 
Thus W’(1) is given in terms of (1) in the same 
way as ¥’(0) in terms of (0), and is symmetric 
in all its dotted and undotted indices. The same 
applies for all ©’(m) for m<r—1. All the ¥’(m) 
are given by means of a relation like (29) in 
terms of (m), where 6(m) is defined in terms of 
(0) in the same way as ¥(m) in terms of ¥(0). 
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For m=r, however, V’(r)=0. This is perhaps 
most easily seen from (27), considering the 
definition of W’(r) in terms of ¥’(0) and noting 
(3) and (16). 

The position then is as follows: (0) has 

(r+s+1) independent nonzero components. To 
this may be added W’(0) depending on (0) 
which has (r+s—1) independent nonzero com- 
ponents, so that by a suitable choice of (0) all 
but two of the components of (0) can be made 
zero. 
Further ¥(m) (m=1, 2, ---, r) has (r+1—m) 
nonzero independent components and ®(m) has 
(r—m) nonzero independent components, so that 
each W(m) has only one “‘really”’ independent 
component. In particular V(r) has only one non- 
zero component, and W’(r)=0 as was shown 
earlier. A similar state of affairs holds for the 
V(—™m). 

It is now clear that the single nonzero com- 
ponents of W(r) and W(—s) cannot be trans- 
formed away by a gauge transformation. But 
these nonzero components are given in terms of 
the components 


| 


2 
v(0) and wW(0) 


of (0) by (24) and (26). 

Hence it is clear that these two components of 
W(0) are gauge invariant and are the two really 
independent components of (0). 

All that remains is to interpret these states. 
With the choise of the matrices ||g’4|| that we 


DE WET 


have made, an axial rotation about the x! axis 
through an angle @ corresponds to the spin 
transformation I'4g given by 


ee 0 


= 


Under this spin transformation 


v(0) 


is multiplied by a factor e}*"+*)*, while 


deed 
is multiplied by a factor e~!*("+")*, 

From this follows that the two components 
of (0) correspond to components of spin 
+3(r+s) along the x! axis. But we have chosen 
the frame of reference in such a way that the x! 
axis coincides with the momentum vector. So 
that the two states correspond to components of 
spin +f along the momentum vector (where 
f=3(r+s)). This is the result given by Fierz. 
A further result of Fierz is that the (total angular 
momentum)? for the system is always 2 f(f+1). 
The reason for this is easily seen. The orbital 
angular momentum for the system always has 
component zero along the momentum vector. 
Thus it follows that the total angular momentum 
always has the components +f along the 
momentum vector, and thus the (magnitude)? 
of the total angular momentum must always be 


S f(F+1). 
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Apparatus is described for scattering electrons with thin 
foils and measuring the intensity of scattering in a narrow 
pencil making any angle from 30° to 120° with the direction 
of incidence. Electrons with energies below 95 percent of 

the energy on incidence are excluded by a retarding voltage. 
The theory of nearly single scattering is discussed, and a 
general expression is found for the error introduced by 
slight plural scattering. In the case of Rutherford scattering 
this expression gives for the fractional error the value 
é{csc?(@/2)— }} where @ is the angle at which scattering is 
measured and é is the mean square deflection of particles 
traversing in the foil the same length of path as those 
scattered at the angle @. The relation of this result to 
Wentzel’s criterion is shown. Electrons accelerated by 50 kv 
were scattered by thin composite foils of Al on collodion. 


A difference was found, though not explained, between the 
intensities of scattering at 90° on the two sides of the foil, 
which was inclined at 45° to the direction of incidence. 
With allowance made for this asymmetry in comparing the 
scattering at angles above and below 90°, relative intensi- 
ties at angles from 30° to'120° were within 5 percent of 
those predicted by Mott. Comparison of the observed and 
predicted absolute intensities of the scattering from Al was 
subject to an uncertainty of about 20 percent, principally 
in the determination of the thickness of the Al film. Within 
this uncertainty, the observation agreed with the predic- 
tion. If the actual scattering differs from the predicted, it 
appears from these observations to be more probably less 


than greater. 


INTRODUCTION 


F a beam of electrons of speed v is incident 
normally on a thin foil of thickness ¢, having, 
per unit volume, m atoms of atomic number Z, 
then the fraction scattered into a small solid 


angle 2 around a direction making an angle @ 
with the direction of incidence may conveniently 
be written 


ntQ(Ze?/2mv*)*?(1 — B*)csc*(8/2)R. 


(1) 


For the factor R, Mott! has calculated from 
the Dirac equations, assuming a Coulombian 
field, the expression 


sin® (6/2) 
+ (ZBxr/137) cos? (6/2) sin (6/2) 


+terms of the second and higher powers in 
Z/137. (2) 

The results of several experiments have indi- 
cated exceptions to this equation and to the 
closely related prediction of an asymmetry in 
double scattering, and it has been suggested that 
they show an actual invalidity of the Dirac 
equations in a range of electron speeds where it is 
hardly to be expected on other evidence.? On the 
other hand, the experiments offer difficulties in 


1N. F. Mott, Proc. Roy. Soc. Al24, 425 (1929). 
* A general discussion of theory and experiment is 
in an article by M. E. Rose, Phys. Rev. 57, 280 (1 


technique, and not all the experimental results 
are in mutual agreement. It seems desirable 
therefore to have more observation on the 
subject. 

In the present experiment tests were made of 
Eq. (2) with 50-kv electrons scattered at angles 
between 30° and 120° by thin composite foils of 
Al evaporated on collodion. 


APPARATUS AND EXPERIMENTAL PROCEDURE 


The chamber in which the electrons were 
scattered is shown in Fig. 1(a) in vertical section 
and in Fig. 1(b) in horizontal section. Its vertical 
wall is a short section of 1-ft. steel pipe, supported 
for insulation on a glass cylinder of the same 
diameter and 2 feet high. The beam of electrons 
enters the chamber, as shown by the arrow, 
through a sylphon bellows and ground joint from 
the accelerating tube (not shown). This beam, 
already rendered parallel by the distribution of 
voltage among the electrodes of the accelerating 
tube, is defined by the diaphragm D,. This dia- 
phragm is at the end of a steel tube set in a 
bracket fastened to the top of the scattering 
chamber. The bracket is flanked by blocks of lead 
to absorb x-rays. The electrons passing the 
diaphragm strike the foil, which is mounted on 
one of two diaphragms carried by a cylindrical 
rod projecting through the top of the scattering 
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Fic. 1 (a) vertical section and (b) horizontal section of 
scattering chamber. The wall of the chamber is 1-ft. steel 
pipe. The support is a glass cylinder. Arrow shows direction 
of incident om. T, turning floor of chamber. D,, dia- 
phragm defining incident beam. D2, ae defining 
scattered beam. D;, latger diaphragm behind D:. Ci, brass 
cylinder enclosing region of retarding voltage. C2, collecting 
chamber for scattered beam. C;, collecting chamber for 
incident beam. B, collecting voltage for primary beam. M, 
microammeter for measuring incident current. F, scattering 
foil. P, collecting electrode for scattered beam. 


chamber. The rod is attached to a mechanism 
(not shown) housed in sylphon bellows above the 
top of the scattering chamber and evacuated with 
it. This mechanism permits placing either of the 
two diaphragms in the electron beam and also 
setting their plane for incidence at 45° to either 
side of normal, all by manipulation outside the 
vacuum. 

The floor of the scattering chamber is a steel 
disk with a circular groove on the under side near 
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the rim. This groove rides on three small wheels 
set on axles 120° apart on the inner wall of the 
scattering chamber. Two of the wheels are idlers; 
the third can be driven from outside the vacuum 
by a crank housed in a flexible metal tube. In this 
way the floor can be rotated about a vertical 
axis through the center of the scattering foil. 

A brass pipe C; projects through the floor from 
below almost to the top of the chamber. In the 
wall of this pipe is a copper cylinder carrying 
two steel diaphragms D, and D3. Inside C, but 
insulated from it by three supporting Pyrex glass 
rods is the copper chamber C; for the collection of 
the beam of scattered electrons. 

The potentials of the parts of the apparatus 
are shown schematically in Fig. 2. The chamber 
C2 was at ground potential. The pipe C; and all 
the parts of the scattering chamber outside it 
were at a potential V which was the highest 
positive potential of the apparatus. It was main- 
tained by a transformer and electron-tube recti- 
fier with a large condenser to eliminate ripple. 
The voltage was found by measuring with a 
microammeter the current through a calibrated 
high resistance. 

The filament was maintained below ground 
potential by means of a power pack producing a 
voltage Vr which was a small fraction of the 
total voltage used to accelerate the electrons. 
The electrons were thus incident on the foil 
with an energy corresponding to the total voltage 
V+ Vr. The scattered beam was defined in field- 
free space at high potential by the diaphragm D, 
(which was considerably smaller than D3). Pass- 
ing through D; the electrons were retarded by 
the voltage V, applied between D; and a grid of 
fine tungsten wire in the opening of C2. Hence 
only the energy corresponding to the voltage V, 
could be retained by any electron passing through 
the grid. For this reason Vr may be called the 
residual voltage. Any electron which lost in being 
scattered an energy as great as that correspond- 
ing to this residual voltage was repelled before 
reaching the grid. 

Not only these electrons were excluded by the 
retarding voltage V but also secondary electrons 
ejected from the foil by the primary beam and 
any that may have been ejected by x-rays from 
any part of the apparatus except the interior 
of C2. 
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SCATTERING OF ELECTRONS 


Electrons passing the grid of C2 were collected 
on the plate P. This was maintained at a positive 
potential of 270 volts with respect to the grid to 
prevent the escape of electrons by recoil or the 
escape of secondary electrons ejected from P. 
The current to P passed thence to an amplifier 
of the Barth type, described by Penick* and the 
amplified current was measured by a shunted 
Leeds and Northrup Type R galvanometer. 

The current incident on the foil was commonly 
about 10 wa. With the precautions taken against 
slow electrons and x-rays it was possible to 
measure scattered currents as small as 10~® ya. 
This allowed the diaphragm D, to be made small 
enough to give a rather high angular resolution. 
It subtended a plane angle of about 3° at the foil, 
its solid parallax, 2 in Eq. (1), being 0.00223 rad.? 

In measuring the fraction of the incident elec- 
trons scattered at any angle through the dia- 
phragm Dz, readings of the incident and scattered 
currents were alternated. To measure the inci- 
dent current, the foil was removed from the 
incident beam by the mechanism above the 
scattering chamber, which was operated from a 
distance by gas pressure, without impairing the 
insulation of the scattering chamber or jarring 
the apparatus. The incident current then passed 
directly to the chamber C; and was measured by 
the microammeter M, which could be read from 
a distance. When the diaphragm carrying the foil 
was removed from the incident beam it was 
replaced by an open diaphragm carefully made 
to be as nearly as possible identical with it. This 
diaphragm was used for ‘“‘blank” readings to 
detect any current to C, other than that of elec- 
trons scattered by the foil. Whatever such cur- 
rent there was, it was too small to measure. 


Fic. 2. Schematic diagram of apparatus showing po- 
tentials of different parts. S, source of electrons, filament of 
accelerating tube. 


, scattering foil. C,, chamber 
n of retarding voltage. C:, collecting chamber. P, 


ecting electrode. Vr+V, total accelerating voltage. 
V, retarding voltage. Vr, residual voltage. 


*D. B. Penick, Rev. Sci. Inst. 6, 115 (1935). 
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The zero of the galvanometer was subject to 
some drift. To eliminate the effect of this each 
reading made with the beam of electrons incident 
on the foil was bracketed between two readings 
with the beam cut off. 


RETARDATION OF ELECTRONS IN THE FOIL 


Any electron traversing the foil must lose some 
energy in remote collisions with other electrons. 
It would obviously be a mistake to set the re- 


Fic. 3. Variation with residual voltage of current to 
collecting electrode (in units of 1 cm galvanometer deflec- 
tion per microampere incident on foil). Scattering at 90° on 
side of incidence. Foil II. 


sidual voltage so low as to exclude from the 
measured scattered current any electrons which 
had lost only such energy as this. To find a proper 
value for the residual voltage, the current 
scattered to the collecting electrode was measured 
with various values of the residual voltage. 
Figure 3 shows the result of observations made 
with one foil (foil II of the later discussion). The 
scattered current was measured at 90°, with the 
electrons emerging from the same face of the foil 
as that on which they were incident. Incidence 
was at 45°, so that an electron traversing almost 
the whole oblique thickness of the foil and being 
reflected through 90° at the far side would have 
a length of path in the foil 2v2 times the thickness 
of the foil. No other conditions of observation 
would provide so great a possible length of path 
and hence so great a possible retardation. 

It is evident from the figure that the energy 
lost by ordinary retardation was less than 1 kev, 
since the current was not increased by the in- 
crease of residual voltage from 1 to 2.5 kev. The 
thickest foil used in the observations on scatter- 
ing was about twice as thick as the one used here. 
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The residual voltage employed was 2.5 kv, which 
is judged more than enough to compensate for 
the ordinary retardation of the electrons. 


TESTS OF SINGLE SCATTERING 


Wentzel* has given as a criterion for single 
scattering of the Rutherford type that the angle 
6 at which scattering is observed shall be several 
times greater than 4w, where 


2V/2\3 
w=2 cot! —(—) 
Ze\rnl 


eV being the kinetic energy of the electron. It has 
been customary to consider scattering as single 
at angles above 12w. In the present experiment 
w, as reckoned from the oblique thickness of the 
foil, was about 2° for the thinner foils used and 
something less than 3° for the thickest. Conse- 
quently the criterion would be satisfied at angles 
around 30° and higher. Instead of depending 
wholly on Wentzel’s criterion, however, we have 
applied a more empirical test of single scattering 
derived from the following considerations. 

When the foil is thin enough for most of the 
electrons to pass through with only small de- 
flections, it is very improbable that one electron 
should undergo two large deflections. The devia- 
tions from strictly single scattering at any large 
angle will rather be due almost entirely to the 
combination of a single deflection through an 
angle nearly equal to this one with one or more 
much smaller deflections. 

In Fig. 4 let OP be the direction of the beam 
of electrons incident on the foil. Let OQ be the 
direction of the observed scattered ray making 
an angle @ with the direction of incidence, OP. 
Let it be assumed that an electron observed in 
the ray OQ entered the ray not by a single de- 
flection from OP to OQ through the angle @ but 
by the combination of two deflections, one from 
OP to OP’ through a small angle ¢« and another 
from OP’ to OQ through the angle @’ nearly equal 
to 6. The angle @ between the planes of « and 6 
may be taken at random. Also the sequence of 
the two deflections has no significance in the 
result ; « may even be considered as the resultant 
of several small deflections of which some may 
precede and some follow the large deflection. 


*G. Wentzel, Ann. d. Physik 69, 333 (1922). 
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Fic. 4. Combination of small deflection ¢ with large 
deflection 6’ to ore deflection at angle of observation 8, 
OP, direction of incidence. OP’, direction after deflection ¢, 
OQ, direction after deflections « and 6’. @, angle, assumed 
random, between planes of ¢ and @. 


If powers of « above the second are neglected, 
6’ =0—ecos cot 6 sin? 


Let p(@) be the probability of single deflection 
through the angle @ from the ray OP to the ray 
OQ, and let p(@’) be the probability of single 
deflection through the angle 6’ from the ray OP’ 
to the ray OQ. Then to the same degree of ap- 
proximation as before 


p(0’) = p(@) cos }¢? cot sin? &) 


Let p’(@) denote the average probability of 
deflection into the ray OQ of an electron initially 
incident in the direction OP but deflected through 
the angle ¢«, the average to be taken with all 
values of & equiprobable. Then 


1 
=— 0’ )d®. 
- 


With the expression given above for (6), 
we obtain 


To find the error made when the actual, ap- 
proximately single scattering is taken as strictly 
single, it is necessary to average this expression 


Fic. § 
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over all values of e. The error is thus the right- 
hand member of the equation with & denoting 
the mean square deflection of electrons having a 
mean length of path in the foil equal to that of 
the electrons deflected at the angle @ from the 
direction of incidence. In our experiment, inci- 
dence being at 45°, the mean length of path was 
not very different among the electrons observed 
and was not greater at any angle than that of an 
electron passing undeflected through the foil. 
Consequently & was not greater than the mean 
square deflection of all the electrons incident on 
the foil, which could be estimated experimentally. 
Equation (3) is valid for any law of scattering 
provided the observed large deflections may be 
ascribed to the combination of a single large 
deflection only with much smaller deflections. 
The derivatives of p may in this case be replaced 
without serious error by those of the experimen- 
tally determined function p’. Thus the experi- 
mental data furnish all the means of the test. 

The observations in this experiment show p’ 
varying very nearly as csc‘ (@/2). In this case 
Eq. (3) gives for the fractional error caused by 
plural scattering : 


(8) — p(8)} /p(0) = csc*(0/2)—3}. (4) 


To estimate ¢’, the scattering at small angles 
was observed with the movable collecting cham- 
ber. The observations made with one foil (foil II 
of the later discussion) are shown in Fig. 5, the 
current collected in the chamber being plotted 
against the angle. Since the diaphragm in the 


2 83 ¢ 6 


Fic. 5. Scattered current at small angles, in microamperes. 
Foil II. The incident current was 10 ya. 


front of the chamber subtended approximately 
3°, the resolution was low at these small angles, 
and é could not be determined accurately. But 
since the collected current at 0° was 0.3 of the 
incident current, it is clear that 0.3 of the incident 
electrons were not scattered through more than 
1.5°, half the angle subtended by the diaphragm. 
From this and the shape of the curve, it seems 
reasonable to take 2° or 0.035 rad. as the value 
of «. With this value, the error caused by plural 
scattering would be 1.8 percent at 30° and 0.18 
percent at 90°. With the thickest foil used, the 
errors would be not quite twice as large as with 
this one. 

From the preceding calculation it may be seen 
why the simple form of Wentzel’s criterion, 
which states that @ shall be at least equal to 
about 12w, serves as well as it does. From Eq. (4) 
it follows that, if the fractional error is not to 
exceed a certain value 6, then it is necessary that 


6 1 
(6/2) 


Now @[csc?(6/2) —4]! varies only slowly with 
6. From the value 2 when @ is zero it goes to a 
minimum of 1.9 when @ is about 0.67 and rises 
to 2.2 when @ is x. To set a maximum for the 
fractional error is thus approximately the same 
as to set a minimum for 6/e (about 14 for a 
maximum error of 2 percent). 

If the coherence of the scattering is not impor- 
tant, € is proportional to (nt)!. So also is w if it is 
small. Hence setting a maximum value for the 
fractional error is also approximately setting a 
minimum value for 6/w, as is done in applying 
Wentzel’s criterion. 

If Wentzel’s criterion is to be valid, however, 
w must be reckoned with ¢ taken as the mean 
length of path in the foil of the electrons scattered 
at the angle @, not as that of an undeflected 
electron. A similar point has already been made 
in regard to «. The exceptions to Wentzel’s 
criterion reported by Neher’ may probably be 
ascribed to the fact that the average path of the 
scattered electrons in his experiment was much 
greater than the foil thickness, from which w 
was computed. 


5H. V. Neher, Phys. Rev. 38, 1321 (1931). 
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Fic. 6. Variation with angle of relative scattered current. 
Plotted points show averages of measurements with three 
foils. Curve shows Mott's theory. 


RELATIVE SCATTERING AT VARIOUS ANGLES 


Three foils were used in the final observations, 
all made in the same way. Al wire was wrapped 
around a tungsten filament and evaporated in 
vacuum by heating the filament. Some of the Al 
was allowed to condense on the very thin col- 
lodion film used as the support for the scattering 
foil. 

The collodion film was too thin to allow a 
measurement of its thickness, but it was not so 
thin that its contribution to the scattering could 
be ignored. Consequently in order to compare the 
absolute intensity of the scattered current with 
that predicted it was necessary to determine the 
difference in the scattering by a film before and 
after a known layer of Al was deposited on it. 
Moreover the scattering by the collodion alone 
could not be measured directly, for uncoated 
films of collodion were quickly destroyed by the 
beam of electrons. Probably they became so 
highly charged that they were torn apart by 
electrostatic forces. 

Foils I and II were collodion films thinly coated 
with Al. Foil III was the same foil as II but with 
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an additional coating of Al of about the same 
thickness as the first. The difference in the 
scattering by foils II and III was ascribed to the 
second layer of Al. The determination of this 
difference was thus subject to the inaccuracies 
of the observations with both foils. The measure- 
ment of all absolute intensities is also subject of 
course to inaccuracies of calibration and to the 
considerable uncertainty in the determination of 
the thickness of the Al layer. 

The comparison of the relative scattering at 
different angles with the prediction of the theory 
was free from several of the errors of calibration 
and from the error in measuring the thickness of 
the foil. Also it did not require taking the differ. 
ence between the intensities observed with foils 
II and III. For elements of atomic number as 
low as those of Al and the constituents of col- 
lodion, the third term in the expression for R, 
Eq. (2), is small, and the difference in its extreme 
values between 30° and 120° does not amount to 
more than 2 percent of R. Since the probable 
errors of even the relative intensities at different 
angles were greater than this, the third term in 
R was ignored in the comparison. The other two 
terms do not involve Z and are the same for 
collodion as for Al. Hence the data obtained with 
the three foils could be averaged together. 

With these three foils and others used in pre- 
liminary observations, an effect was observed for 
which we have as yet no explanation. It appeared 
in a comparison of the intensities of the beams 
scattered at 90° on the two sides of the foil, the 
side on which the beam was incident and the 
opposite side. Single nuclear scattering should be 
equally intense on both sides. But the observed 
scattering was consistently more intense on the 
side on which the beam was incident. The ratio 
of the two intensities was 1.1 for the thinner foils 
I and II and 1.2 for the thicker foil III. 

It is hard to see what sort of scattering could 
combine with the nuclear scattering to cause 
such a difference. Some kind of boundary reflec- 
tion, if it could be imagined great enough to be 


TABLE I. 


6 30° 40° SoO° 60° 70° 80° 100° 110° 120° 


tg/ig obs. 60 19.5 8.2 4.1 2.31 1.50 0.69 0.53 0.44 
ie/iw pred. 60 19.6 8.3 4.2 2.38 1.49 0.71 0.54 0.42 
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SCATTERING OF ELECTRONS 


detectable, might produce a difference in this 
sense, but the absolute value of this difference 
should remain constant once the boundary was 
established, and the relative value should di- 
minish with increasing thickness of the foil. The 
observed difference varied in the opposite way. 
Any effect of the coherence of the scattering at 
different scattering centers might increase in its 
relative value with increasing thickness of the 
foil, but it is hard to imagine either that co- 
herence would be important at 90° with 50-kv 
electrons or that coherent scattering would be 
asymmetrical. 

On the other hand we have not been able to 
ascribe the asymmetry to any defect of the 
experiment. Tests already described seem to 
make it unlikely that it was caused by scattering 
from parts of the apparatus other than the foil, 
or by a selective action of the retarding voltage 
on electrons traversing different lengths of path 
in the foil, or by plural scattering. To test for any 
accidental asymmetry in the apparatus, the foil 
was turned through 90°, so that its sides of inci- 
dence and emergence were reversed with respect 
to the two sides of the scattering chamber, but 
the asymmetry persisted. 

Measurement of the scattering on both sides 
of the foil could only be made around 90°. Over 
the higher range of angles, scattering could be 
measured only on the side of incidence, and over 
the lower range only on the other side. Conse- 
quently in comparing the scattering above 90° 
with that below allowance for the observed 
asymmetry could only be made in some rather 
arbitrary way. 

We have followed the simplest procedure 
possible by comparing the scattering measured 
at every angle with that observed at 90° and 
using the scattering measured at 90° on the side 
of incidence for the comparison at angles above 
90° and the scattering measured at 90° on the 
side of emergence for the comparison at angles 
below 90°. 

Table I shows the observed values of i¢/i99 with 
those predicted by Eqs. (1) and (2). The observed 
values are averages for all the readings with the 
three foils in which every reading was given the 
same weight except that 7 out of 260 were dis- 
carded, either because they were very far from 
the mean or because of some suspicious circum- 


e° 
30 90 


Fic. 7. Variation with angle of Re/Rw, Eq. (1). Dots, 
foil I. xxx, foil II. Crosses, foil III. Circles, average. Curve, 
Mott's theory. 


stance in the behavior of the apparatus noted 
when they were made. 

The same data are shown in Fig. 6, in which 
the points show the observations and the curve 
shows the prediction. To the scale of this figure 
only the average values obtained with the three 
foils can be shown, as the points representing the 
separate observations fall too near together to 
be plotted. Although the deviations of the plotted 
points from the curve are scarcely discernible, 
it should be noted that the graph of i¢/i99 against 
6 does not give a very critical comparison of 
observation with prediction, since the factor 
csc‘(6/2) must provide the dominant variation 
in any theory of scattering in a Coulombian field. 
A stricter comparison is made in Fig. 7, in which 
Re/Roo is plotted against 6. The curve is plotted 
from Eq. (2) with neglect of terms after the 
second. The plotted points show the observations 
made with the three foils and their average. The 
curve passes about as near to the open circles, 
which show the average values of the measure- 
ments, as any smooth curve could. The agree- 
ment is, indeed, perhaps closer than can be con- 
sidered significant in view of the spread of the 
measurements with the foils separately. 


ABSOLUTE INTENSITY OF SCATTERING 


To obtain the scattered intensity as a fraction 
of the incident current, it was necessary to find 
what part of the current passing through the 


249 
ne + 
he 1.10 
he 
lis 
re- 10} + + 
of 
he 1.00 
of R, 

Re ° 

ry + 
on 
of 090 
ils 
ol- 
R, 
me 
to 
ble 
ent 
in 
wo 
for 
ith 
re- 
for 
red 
ms 
the 
the 
be 
ved 
the 
tio 
oils 
uld 
use 
jec- 
be 


250 


diaphragms Dz and D3; was intercepted by the 
grid in front of the collecting plate P. This was 
done by swinging the diaphragm D, into the inci- 
dent beam with the scattering foil removed and 
measuring the current to P with a microammeter, 
and alternately swinging the chamber C, out of 
the beam and measuring the current to the fixed 
collecting chamber C;. The current to P was 
found to be 90 percent of that to C3. 

To compare the absolute values of i, observed 
and predicted, it was necessary to determine the 
thickness or, more directly, the mass per unit 
area, of the second Al layer deposited on foil II 
to produce foil III. Two estimates were made. 
The first was by assuming that all the Al wire 
wrapped around the tungsten filament was 
evaporated, as it appeared to be, and that the 
fraction deposited per unit area on the foil was 
1/(4xr?), where r is the distance from the filament 
to the foil. The surface density of the foil so 
estimated was 16.2 micrograms per cm?*, which 
would make the thickness 6.0X10-* cm. Two 
factors may have made this estimate too high. 
The Al may not have been entirely evaporated 
from the filament and some of the atoms striking 
the foil may not have adhered. Of course some 
atoms may have reached the foil after rebound- 
ing from other surfaces, but this possibility was 
limited by baffle-plates above the foil. The 
values of i» predicted from this estimate of the 
surface density are for these reasons more likely 
to be high than low, and the ratio of the observed 
ig to that predicted more likely to be low than 
high. 

The other estimate was kindly made for us by 
Mr. Hermann Yagoda by a colorimetric micro- 
chemical method applied to a layer of Al de- 
posited on a blank collodion film by the evapora- 
tion of the same quantity of Al wire and under 
the same conditions as in the preparation of 
foil IIT. 

A disk of this foil 4.8 mm in diameter was im- 
mersed in 0.5 cm’ of 10-percent hydrochloric 
acid, which was then warmed over steam until 
the film lost its metallic luster. The solution, to- 
gether with water used in two subsequent wash- 
ings of the collodion film, was transferred to a 
micro-Nessler tube. The contents of the tube 
were treated with 0.4 cm® of 50-percent ammo- 
nium acetate and 0.1 cm* of freshly prepared 
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0.5-percent solution of the ammonium salt of 
aurin tricarboxylic acid. The solution was diluted 
to 2 cm* and, 15 minutes later, when the pink of 
the Al lake became fully developed, the sample 
was compared in color with a set of standards 
made by treating standardized Al-ion solutions 
in the same way as the solution of the Al of the 
foil had been treated. 

The surface density of the foil found by this 
method was 13.6 micrograms per cm®*, which 
would make the thickness 5.0 10-* cm. 

Other samples of the Al layer meant to be 
measured by another chemical method were un- 
fortunately destroyed in the process without 
yielding a result. 

In reckoning the difference in the currents 
scattered by foil III and foil II, which is to be 
ascribed to the additional layer of Al, only the 
readings at angles from 30° to 60° have been used, 
since the random errors of the measurements 
were relatively larger at the higher angles. The 
ratio of the observed scattering to the predicted 
is 0.80 if the surface density of the foil is reckoned 
by the inverse-square calculation and 0.96 if it is 
taken as given by the colorimetric method. Both 
values would be lower if the readings at larger 
angles were included. The first and smaller of 
these values might be expected to be low, as was 
explained before. Either one is uncertain by as 
much as the difference between them. It would 
appear, however, from these observations that 
if the actual scattering is different from that 
predicted by Mott it is more likely lower than 


higher. 


COMPARISON WITH NEHER’S OBSERVATIONS 


Among the many experiments on the scatter- 
ing of electrons, a part of the work of Neher’ 
seems most nearly comparable to ours in respect 
to the experimental conditions. He measured the 
scattering by Al of electrons at voltages around 
that of our experiment and at higher angles above 
90°. He found, as we do, that the relative scatter- 
ing at different angles is in agreement with 
Mott’s theory, though because of the small differ- 
ence between Mott's distribution and Ruther- 
ford’s he could not discriminate between them. 
He found the ratio of the observed to the pre- 
dicted scattering to be 1.32, an apparent differ- 
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PARTIALLY INVERTED MULTIPLETS 


ence from our result given above. Since all his 
observations were made above 90°, where ours 
are least accurate, the two results may not be 
strictly comparable. If they may be compared, 
then our observation of the asymmetry between 
the scattering on the side of incidence and on the 
other side may perhaps be relevant. All of his 
observations were on the side of incidence, and 
if the excess scattering to this side is a general 
characteristic of scattering by foils and not a 
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peculiarity of our apparatus, the difference be- 
tween his observations and ours might be 
explained. 

We wish to thank Mr. A. C. Weid for help in 
the construction of the apparatus and in making 
preliminary observations, Mr. Hermann Yagoda 
for measurements on the foil thickness, our 
colleagues Professors F. E. Myers and O. Halpern 
and Dr. R. D. Huntoon for helpful discussions, 
and Mr. William Werker for technical assistance. 
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With the use of Fermi-Thomas statistical eigenfunctions the separations of the 3s3d *D term 
of Mgl are evaluated, taking into account all the magnetic interactions in the atom and the 
electrostatic interactions with other configurations. The smallness of the observed separations 
is found to be due to the large value of the exchange integrals of the magnetic actions. 


1. INTRODUCTION 


HIS note discusses some remarkable fine 
structures found by Meissner’ in the 
magnesium spectrum (MglI), namely partially 
inverted multiplets showing very small separa- 
tions. It will be shown that they arise from several 
causes whose importance can be evaluated; the 
exact calculation is, however, impossible because 
of the approximations which one is compelled to 
introduce and because of the unavailability of 
accurate atomic eigenfunctions. In any case the 
study of these phenomena in a particular ex- 
ample may be interesting in view of conclusions 
that are generally valid in similar cases. 
We shall discuss the 3s3d *D term, placed at 
13,712 cm~', whose levels, relative to the center 
of gravity of the term, are, respectively : 


3D; at +0.0032 cm-'; at —0.0140 
8D, at +0.0157 


These positions are those which would be ex- 
pected according to the theory if only the spin- 


'K. W. Meissner, Ann. d. Physik 31, 518 (1938). 


spin interaction of the two outer electrons, and 
no other magnetic energy, existed. In fact, ac- 
cording to this hypothesis, the positions of the 
levels would be: 


at *Deat —}La®; *D, at +}Le?’, 


(a= fine structure constant; 


2 d 2 
f 


Rss, Rsa=radial eigenfunctions of the 3s and 3d 
electron, respectively). Thus the ratios between 
the experimental values are very near to the 
theoretical ones, and also the value of the 
integral ZL, evaluated with Fermi statistical 
eigenfunctions, is 


Le? =0.055 


leading to the correct order of magnitude of the 
separations. 

Now the interaction considered is only one, 
and not the most important, of the causes which 
contribute to the splitting of the term. We have 
now to discuss all these causes. 
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First of all the splitting is determined by the 
magnetic interactions in the atom. As is known, 
usually a sufficient approximation is obtained by 
considering only the “‘spin-orbit interaction”’ for 
the outer electrons 


a? 


neglecting all the other magnetic interactions 
and the exchange, and assuming for the potential 
U a suitable expression that takes into account 
the screening of the nuclear field by the inner 
electrons. This approximation would give in our 
case, again using statistical eigenfunctions, the 
value +0.41 for 2(*D) (the “parameter of the 
spin-orbit interaction”), namely a normal triplet 
with separations larger by a factor 30 or 50 than 
the observed ones. 

That means that the approximation is not 
sufficient, as could of course be expected, and we 
have to calculate directly the energies of all the 
magnetic interactions of the electrons outside 
closed shells with the other electrons (and with 
the nucleus) and between themselves. The re- 
sults of these calculations, whose particulars are 
given in Section 2, have no great significance 
because of the inadequate approximation (the 
use of statistical eigenfunctions), but it appears 
that the value of 2(*D) is very much lowered by 
the exchange integrals of the magnetic actions, 
whose contribution to the parameter is negative, 
and which has a relatively large value (—0.27 
cm-'), practically equal to the contribution of 
the classical integrals. (We can introduce a 
parameter analogous to the ordinary parameter 
of the spin-orbit interaction because all the non- 
vanishing magnetic energies (excepting the 
above-mentioned spin-spin interaction of the two 
outer electrons) cause the triplet to be a normal 
one.) 

Among the other causes that can contribute 
to the splitting of the multiplet, hyperfine struc- 
tures do not occur for the atom in question and, 
as the magnetic interactions are very small 
compared with the Coulomb interactions, the 
Russell-Saunders approximation is strictly valid. 
On the other hand, it is found in Section 3 that 
an electrostatic configuration interaction causes 
separations comparable with those produced by 
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the magnetic actions. Also this cause produces a 
normal separation of the triplet and we must 
therefore conclude that, in order to obtain as qa 
resultant the observed practically null separa. 
tion, the sum of the magnetic actions must ip 
reality produce a completely inverted multiplet 
and then that the value of the exchange integrals 
of these actions should be larger than the above 
given value. Then the calculation does not lead 
to definite conclusions about the observed posi- 
tions of the levels. It is, however, likely that the 
sum of all the causes, which would make the 
term a normal one, happens to be zero, and so 
the observed features are due mainly to the spin- 
spin interaction of the two outer electrons, and 
possibly to minor causes which produce, for the 
other interactions, small deviations from Landé’s 
interval rule. 


2. EVALUATION OF THE MAGNETIC INTEGRALS 


As an expression for the interaction of two 
electrons the best approximation is Breit's 
formula,? namely 


a 


+ terms not dependent on the spin, 


where fi, fz are the coordinates of the two elec- 
trons, pi, p2 their momenta, sj, S2 their spins and 
r12 their separation. 

We have to consider the interaction between 
the nuclear field and the spin of the 3d-electron 


=a°Z([r/r*, p/2]-s), 
which gives the first-order perturbation energy 


E,= f Hp 


and the interaction between all the inner elec- 
trons and the 3d-electron; besides the classical 


2G. Breit, Phys, Rev. 34, 553 (1929). 
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integrals E. of these interactions 


we have to consider also the exchange in- 
tegrals E, 


f 


In the preceding sums we include the spin-orbit 
interaction of the two outer electrons which, like 
the other interactions, produces a normal triplet ; 
on the contrary their spin-spin interaction, 
treated already, which produces another kind of 
separation is not included. The sum of the inter- 
actions of the 3s-electron with the inner elec- 
trons vanishes. 

The calculations are rather long: they require 
mainly the use of well-known formulas on 
spherical harmonics ; for the details we may refer 
to a note of David? on the doublets of lithium; 
they are simplified by sum theorems. For in- 
stance, the integral with respect to the spin co- 
ordinates, when we sum over the pair of inner 
electrons with equal space eigenfunctions and 
opposite spins, is in every case simply m, multi- 
plied by a numerical factor (m, refers to the 
3d-electron) ; the result of the integration with 
respect to the angular coordinates, when we sum 
over all the inner electrons with the same m and 
the same /, is likewise m; multiplied by a numeri- 
cal factor; and so for instance E, becomes 


E.= > (21'+1) 


n'l’ 


x f f (n)riedn 


(n’, /'=quantum numbers of the inner electrons). 
The exchange integrals are more complicated, 
especially for /’=1; besides the R’s they contain 
also the dR/dr’s. 

All the integrals have a factor m,m:; thus the 
diagonalization of that portion of the matrix E 
which corresponds to a given configuration n, | is 
obtained exactly as in the ordinary theory: in 
particular the parameter 234 (corresponding to 
the ordinary parameter of the spin-orbit inter- 


*E. David, Zeits. f. Physik 91, 289 (1934). 
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action for the 3d electron) is obtained by sub- 
stituting formally in all the integrals +4 (=5/2 
in our case) in place of m,m;. The required sepa- 
ration parameter of the *D term is 2(*D) = 423. 
Table I summarizes the results, giving the con- 
TABLE I. Contributions to the separation parameter of the 
3s3d*D term of Mgl. 


Interaction E, with the nuclear field 
Classical interaction E, with the inner 
electrons 
Exchange interaction E, with the inner 
electrons 
Total 2(*D) 


=+1.62 cm 
=—1.30 cm 


=—0.27 cm 
= +0.05 cm 


tribution to 2(*D) of the various integrals. It is 
difficult to estimate the size of the error produced 
by the use of the statistical eigenfunctions. It 
can even be rather large, but it probably affects 
the classical and the exchange integrals in the 
same way. In any case, even if the results ob- 
tained, especially the value of 2(*D), are uncer- 
tain, the large value of the exchange integrals is 
ascertained. If we neglect these integrals, we 
obtain for 2(*D) the value 


2(*D) = +0.032 cm 


in fair agreement with the value +0.41 given by 
the ordinary calculation of the spin-orbit inter- 
action. 


3. CONFIGURATION INTERACTION 


The term we are studying is 3s3d *D. It can 
interact only with other even *D terms; the 
nearest is 3p4p *D. It is a negative term to which 
we can ascribe a position‘ of about — 22,000 cm-"', 
and it is almost the only one responsible for the 
perturbation. As it is a normal triplet, it produces 
a normal separation also for the term we are 
studying. 

It may seem strange that there can be an 
appreciable perturbation by such a distant term, 
but it depends essentially on the very small value 
of the required separations; moreover Bacher 
has shown that in the same Mgl spectrum the 
3s3d 'D term is greatly displaced from its normal 


‘ This is the position observed by F. Paschen (Ann. d. 
Physik 12, 509 (1932)), for 3p3d*D which should have 
about the same energy. 

5 R. F. Bacher, Phys. Rev. 43, 264 (1933); 56, 385 (1939). 
As I had the necessary eigenfunctions at my disposal, I 
have calculated again the integral R. of Bacher’s note, and 
I have found the value 23,060 cm™ instead of 21,353 as 
found by Bacher with another kind of eigenfunction. 
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position by the perturbation of a notably distant 
term. Another case of configuration interaction 
has been treated by Hartree and Swirles.*® 

As we have a complete Russell-Saunders 
coupling, each of the levels of the terms interacts 
only with the level of the other configuration 
with the same J, and the interaction is inde- 
pendent of J, as follows from the fact that the 
electrostatic interaction commutes with all com- 
ponents of LZ and S. 

The matrix element connecting the two con- 
figurations is’ 


T =(2/15)[R'(3s3d, 3p4p) — R'(3s3d, 4p3p) 
Let A be the difference between the wave num- 
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Again using statistical eigenfunctions: 


R'(3s3d, 3p4p) = — 3290 
R'(3s3d, 4p3p) = +4620 
= — 2890 


The experimental value of A is ~ 36,000 and that 
of z is 10.8 cm—.8 We have then a displacement of 
the center of gravity of the term 


T?/A =232 cm, 


and the separation parameter of the perturbed 
3s3d *D term is 


2T?/A*?=0.069 


It is of the same order of magnitude as the 


| bers of the centers of gravity of 3s3d*D and_ separations due to the magnetic actions. ing 
i 3p4p *D, and let z be the separation parameter of Terms such as 3p5p *D and higher ones pro- = 
| the 3p4p term; then, neglecting the very small duce only negligible perturbations as compared reg 
i separation of the 3s3d term, we obtain three with that now studied. The (2) 5(3s)*3p *p ba 
if similar secular equations for the three levels, term has a very large separation and is inverted, = 
q e.g., for the *D; term: so that its action is in the opposite sense to, that lov 
| of 3p4p*D; but, owing to its great distance bat 
, =(. (>3 Ry), it has no appreciable effect. cor 
T -—A—2s—x pot 
q From these equations we obtain as displacements 4. CONCLUSION tro 
| of the 3s3d *D levels: The smallness of the separations of the 3s3d *D shi 
J mT? 2 term of Mgl is due to the large value of the —- 
C 
3 exchange integrals of the magnetic interactions 
] "A +22 A between the 3d electron and the inner electrons. shai 
in They probably compensate also the separation ont 
3 produced by the electrostatic interaction of tran 
f 3s3d *D with 3p4p *D which is found to be com- reat 
| Tm OT? 3s parable with that due to the magnetic actions. cart 
| 3 ws—f 14+— The particular position of the levels seem to be enti 
"ee 3z A A due only to the spin-spin interaction of the two true 
outer electrons. It 
6D. R. Hartree and B. Swirles, Proc. Camb. Phil. Soc.  —————_— Paul 
! 33, 240 (1937). 8 We deduce this value from the experimental separations 
\ 7 Notations as in Condon and Shortley, The Theory of of the 2s3p*P and 2s4p?P triplets, remembering that for Prt 
Atomic Spectra (Cambridge, 1935). npn'p *D terms 2(°D) = the < 
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Following the methods outlined in Condon and Shortley’s Theory of Atomic Spectra, the 
energy states of praseodymium IV were calculated. It was assumed in the solid and liquid 
phases that the sharp absorption spectrum would arise from transitions within the incomplete 
4f shell, in other words from states arising from the f? configuration. Excellent agreement was 
found with the experimentally observed values of Pr IV. 


HE absorption spectrum of praseodymium 
ion in solutions and in the solid state has 
been studied by a great many investigators dur- 
ing the past forty years. The spectrum always 
consists of three narrow bands in the blue-green 
region and an additional band in the red. These 
bands show some structure at room temperatures 
and resolve into many extremely sharp lines at 
low temperatures. While the positions of these 
bands in the spectrum remain approximately 
constant regardless of the praseodymium com- 
pound and are therefore to be attributed to elec- 
tronic transitions involving the Pr ion, they do 
shift slightly, and the fine structure of the bands 
may vary considerably. 

On account of the weak intensity and extreme 
sharpness of these bands Van Vleck! concluded 
on theoretical grounds that they must be due to 
transitions between electronic states arising from 
rearrangement of the 4f electrons of the rare 
earth ion. In other words they are transitions 
entirely within the incomplete 4f shell. If this is 
true and if Russell-Saunders coupling is assumed, 
it can readily be shown by the application of 
Pauli’s principle that the only states allowed for 
Pr+++ are *F, *P, ‘J, 'G, \D and 'S and that 
the absorption bands must arise from transitions 
between these states. Hund’s? work on magnetic 
susceptibilities shows that at least for the lowest 
state Russell-Saunders coupling is closely ap- 
proximated. 

From analogy with the La II spectrum Ellis* 
concluded that the blue bands of praseodymium 
arose from a *J7,4—*Po,;,2 and that the red band 


x H. Van Vleck, J. Phys. Chem. 41, 67 (1937). 
*F. Hund, Zeits. f. Physik 33, 855 (1925). 
°C. B. Ellis, Phys. Rev. 49, 875 (1936). 


arose from the *H,-'J, transition. Bethe and 
Spedding* next showed that the spectrum of 
Tm IV could be completely accounted for on 
theoretical grounds. This is the companion rare 
earth to Pr IV and has the same energy states 
except that they are inverted due to the incom- 
plete 4f shell being more than half completed. 
They assumed that the bands were due to transi- 
tions between states arising from the 4f?=4f? 
configuration. Following the methods outlined in 
the book by Condon and Shortley,® they calcu- 
lated the separation of the levels due to electro- 
static interaction and then took into account the 
spin-orbital interactions including perturbations. 
It was found for the upper states that Russell- 
Saunders coupling was far from being obeyed, 
since in some cases the perturbations amounted 
to as much as 10,000 cm~. The theoretically 
predicted bands, however, are in excellent agree- 
ment with the experimentally observed ones, the 
discrepancies in all cases being much smaller than 
might reasonably be expected when one con- 
siders the very simple assumptions which were 
made in the theoretical calculations. This holds 
also for the four new bands® which now have been 
observed, but which were not reported at the 
time the predictions were made. 


* H. Bethe and F. H. Spedding, Phy. Rev. 52, 454 (1937). 

5 E. U. Condon and G. H. Shortley, The Theory of Atomic 
Spectra (Cambridge 1935), p. 207. 

Selection rules state that the *H,—*P, *H,—*P» should 
be absent or at least very faint. Gobrecht, Ann. d. Physik 
31, 600 (1935), finds bands between 34,000 cm™ and 38,000 
cm™ which he attributes to Tm IV but states his salts 
contained Gd as an impurity, Meehan and Nutting, J. 
Chem. Phys. 7, 1001 (1939), fail to report any bands in 
this region. Spedding and Gladrow in work now in progress 
have photographed Tm,.(SO,);-8H:O crystals 2.5 mm in 
thickness and find no trace of these Sendo. They are 
neni to photograph much thicker crystals in the near 
uture to see if any trace of them exists. 
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Lueg,’ Freyman and S. Takavorian,* Go- 
brecht® and Rosenthal,” have investigated the 
infra-red spectrum of praseodymium and have 
reported bands at about 5200, 6500 and 9800 
cm~, Lange" by arbitrarily evaluating the F 
constants in the formulae for the f? configuration 
given by Condon and Shortley” in their electro- 
static-interaction calculations, forces a rough fit 
with experimentally obtained levels. It is tacitly 
assumed that good Russell-Saunders coupling is 


TABLE I. Electrostatic interaction. 


Lall Priv. TmIV 

3H 0 0 0 
oF 1,718 4,120 5,670 
1G 2,750 6,600 9,080 
1D 3,260 7,820 10,750 
Wy 5, 14,330 19,700 
3p 8,015 19,230 26,450 
1S 13,756 33,000 45,400 


obeyed by all levels, and that the perturbations 
can be neglected. Gobrecht!*® extending Lange’s 
work does the same thing for Tm and also obtains 
a forced agreement. He objects to the spectro- 
scopic designations given the bands by Bethe and 
Spedding‘ on the grounds that they used the 
wrong screening constant. This objection is 
groundless, since Bethe and Spedding calculated 
the spacing of the levels for various values of 
Condon and Shortley’s’ ¢ which is directly 
related to the screening constant o and chose the 
best value of ¢ to fit with experiment. This value 
is in excellent agreement with the o determined 
from x-ray measurements. Further the designa- 
tions given by Gobrecht must obviously be 
wrong since he has levels with the same J quan- 
tum number crossing each other which arose 
from the same electronic configuration. There- 
fore, the perturbations which he neglected must 
be extremely large and in fact do amount to as 
much as 10,000 wave numbers. It also follows 
that most of the conclusions in these papers 
regarding screening constants of the rare earths 


7P. Lueg, Zeits. f. Ng 39, 391 (1926). 

§ R. Freyman and S. Takavorian, Comptes rendus 194, 
963 (1932). 

*H. Gobrecht, Ann. d. Physik 28, 673 (1937). 

10 Gottfried Rosenthal, Physik. Zeits. 40, 508 (1939). 

1H. Lange, Ann. d. Physik 31, 609 (1938). 

2 E. U. Condon and G. H. Shortley, reference 5, p. 207. 


13H. Gobrecht, Ann. d. Physik 31, 600 (1938). 


F. H. SPEDDING 


are invalid on account of the neglected per. 
turbations. 

Mukherji' also gives a designation to the 
transitions of Pr IV and has done this by ex. 
trapolating the spacings of the observed levels 
of La II to the larger spacings which they would 
occupy in the stronger fields of Pr IV. He has 
therefore taken into account the perturbations 
as they exist for La II, but has neglected the 
greatly increased perturbations which would 
occur for some levels in Pr IV, and therefore gives 
the wrong designation for them. 

In this paper the levels for the free ion of 
Pr IV were calculated as follows. For the electro- 
static interaction, between the electrons, the 
calculated value of La II which is analogous to 
Pr IV was taken as a basis. It was assumed that 
the spacing between the energy levels without 
spin orbital interaction is simply multiplied by a 
constant factor on going from: La II to Pr IV, 
This arises from (a) the higher degree of ioniza- 
tion (b) the greater internal nuclear charge. As- 
suming Ellis’s designation to be correct this turns 
out to be about 2.4 (Table 1). 

The constant {' characteristic of the spin- 
orbit interaction was calculated for La II from 
the splitting of the H level which has most nearly 
Russell-Saunders coupling and was found to be 
about 170 cm™. It is assumed that ¢ increased 
again by a factor of 2.4 from La II to Pr IV and 
should therefore be about 400 cm for Pr IV. In 
column 2 of Table II the positions of the Pr 
levels are given for this value of ¢ when the 
perturbations are omitted, in other words, pure 
Russell-Saunders coupling is assumed to exist. In 
column 3, the positions of the levels are given 


TABLE II. Spin orbital interaction with and without 


perturbations. 
=400 cm™! = —1,400 
Pri Priv TmIV 

WITHOUT WITH 4 WITHOUT WITH 4 
3H, —2,400 —2,639 +239 | *He -7,000 —7,490 +490 
—400 —400 0 1,400 1,400 0 
He 2,000 41,923 +77 | 8400 1,925 +10,325 
Fs 2,520 1,857 +663 | 1,470 5,780 — 4,310 
5F; 3,720 3,720 0 5F; 7,070 7,070 0 
5,320 4,391 +929 | *F2 ~=:11,270 3,250 +8,020 
6,600 7,767 —1,167 | 9, 15,090  -—6,010 
7,820 8,239 | 10,750) 16,240 — 5,490 
Ue 14,330 14,407 —77 | 19,700 201909 
5Po 18,430 +509 | 25,000 27,530 —2,530 
Py 18,830 8,830 0 3P; 27,800 27,800 0 
19,630 19,873  —243 | 200 24,665 +-4,535 
1So 33,000 33,509 —509 | 'So 45,400 49,935 —4,535 


“ P, C. Mukherji, Ind. J. Phys., January, 1938, p. 399. 
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after perturbations are taken into account. This 
was done according to the methods outlined by 
Condon and Shortley'® or M. H. Johnson, Jr.’® 
The secular equations for the f? configuration are 
given below. 


6 
J=6 sy, | H+5¢ 
J=5 | H-¢ 
°F, 
—2(10/3)'¢ 0 
J=4 | —2(10/3)'% G 2(11/3)'¢ 
0 2(11/3)'t F+3¢ 
J=3 *F; | F-$ 
1D» 
3F, | F—4¢ —2(6)'¢ 0 
J=2 | —2(6)% D 3(2)8¢ 
0 3(2)i P+ 
J=1 
3Po 1So 
~~ | S 


These equations were solved for ¢ equals 350 
cm~!, 400 cm and 450 cm-', and the results 
compared with experiment. The best agreement 
was found for a value of ¢=400 in good agree- 
ment with our preliminary assumption. The 
agreement with experiment was exceptionally 


% E, U. Condon and G. H. Shortley, reference 5, p. 267. 
1M. H. Johnson, Jr., Phys. Rev. 38, 1628 (1931). 


TABLE III. Comparison of calculated and observed terms, 


Priv Priv 

CaLc. Os. TmIV TmIV 

cm cm" CaLc. 
1H, 0 0 Hy 0 0 
—«2,239 5,565 
4496 4,750\* | *H, 8.490 8,248t 
4562 5,650 10,740 —-12.750¢ 
6359  5.950\** | 13.270 14,600 
1G, 10,406 9,469 4 22'580 
1D; 10,878 10,200 1D, 23,730 
«17,046 16,630 27,680 28,000 
20,560 20,660 «32,155 35,000t 
21,469 21,150 35,020 36, 300 
«22,512 sP, 35,290 38,000t 

24.050 

1S, 36,144 1S» 57,425 


*Maximum at 5180 
_ **Maximum at 6557 cm~ shows a distinct indentation on short wave 


~ om bands were not observed when levels were first calculated. 
See reference 6. 

good (see Table III), the deviations in every case 
being much smaller than one might expect con- 
sidering the simple assumptions used. Further in 
every case where a band was predicted one was 
found except for the *H,—*H; transition which 
occurs so far in the infra-red that that region of 
the spectrum has not yet been investigated. In no 
case was a band found where none was predicted. 
Also the interpretation of the spectrum of 
praseodymium is completely consistent with the 
interpretations of the spectrum of LalII and 
Tm IV. By shifting the arbitrary constants 
slightly (2.4 for the electrostatic interaction and 
400 cm-! for {) an apparently better fit with 
experiment might have been obtained. It would 
not, however, have had any real meaning unless 
the fine structure of the multiplets due to crystal 
splitting had been investigated and other possible 
perturbations taken into account." 

17 The multiplets due to crystal splitting are always 
associated with superimposed crystal and molecular fre- 
quencies. This causes the multiplets to be composed of 
many lines. As the sharp absorption spectra are obtained 
at low temperatures, the observed associated lines are 
usually due to vibrations superimposed on the excited or 


upper state and therefore the lines usually occur on the 
violet side of the multiplets. 
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Effect of Hyperfine Structure upon the Magnetic Rotation of the Plane of Maximum 
Polarization of Resonance Radiation 


BERTRAND J. MILLER 
State University of Iowa, Iowa City, Iowa 
(Received April 23, 1940) 


An equation has been developed for the polarization, and the angle of maximum polarization 
of resonance radiation as observed along a magnetic field perpendicular to the electric vector 
and the direction of propagation of the incident beam. Explicit application has been made to 
the polarization and the angle for the line 'So—'P, for various values of the ratio of the hyper- 
fine separation constant to the natural breadth of the levels, for a nuclear moment of }. It was 
noted here that the form of the curves might differ widely from that to be expected from classi- 
cal considerations. The polarization, and the angle of maximum polarization for the sodium D 
lines, for several values of the same ratio, were computed and compared with experimental 
observations, Thus a range of possible values of A was determined, which is not in disagreement 
with previous figures, and a value for r established which agrees well with that calculated from 
optical dispersion data. 


INTRODUCTION 


HE polarization of resonance radiation in magnetic fields has been extensively studied, and 
since the realization of the role played by hyperfine structure, most of the early discrepancies 
between theoretical predictions and experimental results have been removed. However, one experi- 
ment of some interest has not yet been compared with modern theory ; i.e., the rotation of the plane 
of polarization by a magnetic field perpendicular to both the direction of propagation and the electric 
vector of the exciting beam, the emitted radiation being observed along H. 

Experimental data for the polarization and the azimuth of the plane of maximum polarization 
are available for the D lines of sodium ;! a theoretical equation for the effect has also been published,’ 
but needs some extending, as it is based on the assumption of a separation constant large in com- 
parison with the natural breadth of the levels, whereas, in this experiment, they are of the same order. 

It seemed worthwhile to undertake the comparison of this experiment with theory, as it should 
lead to values of the hyperfine separation constant A, and the mean life, t which could be compared 


with those available from other sources. 
THEORY 


The most convenient starting place is Breit’s? Eq. (159) 
= (C/ga) (C8) gm my (CE) —2xir 


Here J; is the intensity of radiation polarized in the direction £, C is a constant depending on the 
experimental conditions, 1 and m represent sub-levels of the upper state (b) and the lower (a), 
respectively, and g, is the statistical weight of the lower state, and is also a constant throughout 
this paper. 7 is the mean life of 6, and »v,, is the energy difference between the states »’ and x 
expressed, of course, as a frequency. The incident radiation is supposed to be broad compared to the 
emitted ; and is polarized in the x direction. 

Let us equip the matrices with more descriptive subscripts ; e.g., if the states » and m are those which 
would, in a weak field, be characterized by the quantum numbers F=¢, Mr=yu and F=f, Mr=m, 
the matrix element ez,» is written ezf,,. Further, let the axis of z be that of the magnetic field, and 
y be the direction of propagation of the exciting beam. We choose the standard representation for 
the matrix elements, ez being diagonal, and e(#+iy) is zero unless n.=m-+1. Finally, iet @ be the 


angle bet) between x and é. 


~ 1A, Ellet Ellett, a: Opt. Soc. Am. and Rev. Sci. Inst. 10, 427 (1925). 
2G. Breit, Rev. Mod. Phys. 5, 91 (1933), especially Art. 4, p. 117 et seg. 
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am Substitution into (1) and simplification gives | 


(1) 


Differentiating J; with respect to 6, and equating to zero, one has, for the angle of maximum | 


polarization 
where 


(3) 


f , , 
Br” =(C/ge)| Ye - 


As is customary, the frequency differences and matrix elements required will be obtained from 


ind the Goudsmit-Bacher* modification of the Darwin‘ method. It is necessary to calculate the ampli- 
cies tudes of the electric moments themselves, however, rather than the transition probabilities, owing 
eri- to the appearance of quantities such as (ex)%, (ex)/%, which are not transition probabilities. 
ine A formal simplification is perhaps not effected, but calculation is facilitated if the sum in the 
tric denominator of P, is split into two, parts. Noting that a,**’ is symmetric in ¢ and yw, we may write 
ion ay /1+(2 as” /14+(2 
Introducing the definition 26= we have 
THE TRANSITION 'So— 
To illustrate the method and to lend support to a peculiar feature of the D line curves, calcula tions 

have been made on the simpler transition 'Sp—'P:, and the assumption of a nuclear moment of }. 
P Following Goudsmit and Bacher, we assume that the perturbed functions may be written 
JF JF J Mi 
ly where the ¥, are the characteristic functions of the atom with no nuclear spin. The energy dif- 
he ferences and the X’s are then determined by the chains of equations, 

JF JF JF 
ich — XM Ms+1)(1+ M1 — 
- together with the normalizing equation, 


(XMoMy) (I+ (I — M1) (J+ Mz) (J — My) !=1. (6) 


+S. Goudsmit and R. F. Bacher, Phys. Rev. 34, 1499 (1929). 
*C. G. Darwin, Proc. Roy. Soc. A115, 1 (1927). 
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Fic. 1. P, and angle of rotation of plane of maximum polarization of line 'S)>—'P,, as functions of 2. 
The matrix elements are now given in terms of the X’s by the equations 
Mr, Mpst =} XM Ms) (J— Mo) (1) 


In the above equations, w= g(J)eH/4rmc. 
By equating the determinants of the chains to zero one obtains the quantities vy,/A as functions 


of 2(=w/A). The v3/,/A will hereafter be donated by EX, the J superscript being unnecessary, as 
the only energies needed are those of the upper level, and these will all have the same J in the cases 


treated. 
The chains and the normalizing equation, together with the Z’s found from the determinants, 


determine the X’s, and from these last, the matrix elements may be calculated. 
Inserting the quantum numbers of the lower (1So) level in Eqs. (5), (6), and (7), one gets no splitting 
with H, and all the X%3/,4,=1 independently of H, because the coupling of M; and H has been 


neglected. 
The upper level chains contain two members at most, so that the E’s are the roots of quadratic 


equations, making a calculation in closed form possible. Inserting the values of the Z’s and the X's 
in (7), and the resulting matrix elements in (3), one has 


where , K=2rrA, 


It will be observed that the terms in the 8 sum, which is the numerator of P, and the denominator 
of tan 20 may have sharp maxima, while those in the 7 sum, which is the numerator of tan 2@ may 
have both maxima and minima. The maxima for the terms in the 8 sum occur approximately where 
the frequency difference in the denominator vanishes (an impossible occurrence in the terms of the 
@ sum, since the frequency differences there are between levels of the same Mr); these maxima are 
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Fic. 2. P, and angle of rotation of plane of maximum polarization of lines *S;—*P4,, as functions of H. 


very sharp for large K. The terms have maxima approximately at KE¥/;~,=1 and minima at —1, 
and pass through zero where the frequency differences vanish. Observing the roles of these sums in 
P, and tan 26, one might be led to expect “humps” in the curves in the neighborhood of values of 
Q at which two levels cross, especially for large values of K. 

This expectation is borne out, as is shown in Fig. 1. The term responsible for the behavior is that 
connected with b, which vanishes at Q=1. 

The scales for tan 2@ are not the same in the various curves; in the curve for 2r7A =0.01, the 
ordinates of tan 2@ have been multiplied by 50, and in the curve for 2x7A =2, these ordinates have 
been divided by two, to facilitate comparison of the curves. In Fig. 2, however, all curves have the 
same scales. 


Tue D LINEs 


Considerably more labor is involved in the calculations for these transitions ; fortunately, much of it 
had already been accomplished. Thus, it had already been shown’ that D; was unpolarized for all 
fields and any value of J that might be assumed, so that Dz alone needed to be considered. In the 
same paper, a closed expression was given for the X’s of the lower levels; this enabled their elimina- 
tion from the expressions for the probability amplitudes. The energies and the X’s weré available 
for 2=1 and Q=v2. The secular determinants of the chains give rise to a quartic, two cubics, two 
quadratics, and two linear equations; solutions of the cubics for 2=} were also available. 

By eliminating the lower level X’s and considering D; alone, it is possible to omit the superscript 
J from the X’s, and write 


The equation for the polarization in zero field takes a very simple form, viz. ; 
Y+0.81 y 0.45 1.425 0.315 


«144K? 149K? 14+25K" 
* A. Ellett and N. Heydenburg, Phys. Rev. 46, 583 (1934). 
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TABLE I. Polarization in percent of Dz and of the doublet. 


2nrA 0 1 2 2.1 2.2 2.3 3 
P,(D2) 60 24.24 20.46 20.32 20.21 20.11 19.65 18,97 
P(D,+Dz2) 50 19.71 16.59 16.48 16.39 16.31 15.93 15.37 


The polarization in zero field is seen to be determined by the product 2x7A which may be con- 
sidered the ratio of the separation of the lines to their natural breadth. Table I shows the polarization 
of D2, excited by D2 and of the doublet, excited by the doublet. The theoretical intensity ratio of two 
to one in the incident beam is assumed. 

The energies, or rather the quantities and the quantities a,**’, which depend 
on H and A only through @ have been calculated* for 2=0.05, 0.15, 0.25, 0.50, 0.75, 1, V2. The 
E’s are known to one part in the fifth place; the others are probably not significant beyond the 
fourth figure. To calculate P, or @, it is necessary to assume a value for 277A; these calculations 
have been made for 2x7A =2, 2.1, 2.3, another significant figure probably being lost in the process, 
because of the number of terms summed. The polarizations for zero field corresponding to these 
values are, as seen from Table I, 16.6 percent, 16.5 percent, and 16.3 percent. The observations of 
Ellett! and Larrick’ show that the most probable value of the initial polarization lies within this 
range. These curves are graphed as functions of H and compared with the experimental points in 
Fig. 2. Here one has another parameter to adjust for the best fit; 2x7A having been chosen, one 
can still select a value for A, which governs the scale of abscissae of the theoretical curves, since 
Q=g(J)eH/4rmcA. The values of A used in the graphs shown range from 6.72 to 8.1610- cm-!. 
Since, in all cases but one, the polarization predicted for zero field was higher than that observed in 
this experiment (which seems likely in the light of later observations) ; and since the theoretical curve 
has a horizontal tangent at H=0; and finally, since the curve is insensitive to small variations in 
A in the region of small H, it is impossible to fit the theoretical curve to the experimental for any 
reasonable value of A. Thus, this parameter has been chosen in each case primarily to give a good fit 
with the rotation curve in the region of small H. 

The curves for the rotation of the plane of polarization may be said to fit within the experimental 
error, which is admittedly quite large for H >4 gausses. The same behavior is noted as in the simpler 
case, the levels which cross in this instance being those characterized by the quantum numbers 
F=1, Mr=0, and F=2, Mr=-—2. It is perhaps worthy of mention that the hump in the curve 
disappears for smaller values of 2x7A. The form of the curve is not well determined for large H, 
due to the paucity of calculated points. This is of little importance, however, as the theoretical 
and experimental curves are not being compared in this region. 

Theory and experiment agree to approximately the same degree for all the values of 277A pre- 
sented here; thus these experimental data will not fix A with any precision. Qualitative agreement is 
indicated, however, with the value of 6X10-* cm~! obtained by Ellett and Heydenburg,*® from 
calculations for the case in which the magnetic field is parallel to the exciting beam. 

The curves are much more sensitive to variations in 7. The total range of +r is from 1.50 to 1.58 
10-8 sec., which value is in excellent agreement with that which can be calculated from optical 
dispersion data,’ namely 1.54 10-8. 

In conclusion the author wishes to thank Professor A. Ellett for helpful suggestions throughout 
the progress of the work. 


* The a's for 2= v2 have not been calculated. 


7 L. Larrick, Phys. Rev. 40, 1041A (1932). 
8S. A. Korff and G. Breit, Rev. Mod. Phys. 4, 471 (1932), especially Table XX, p. 499. 
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Optical Dispersion and Molar Refraction at Zero Frequency for Compressed Nitrogen, 
Argon, and Carbon Dioxide Measured as Functions of Density 


CLARENCE E. BENNETT 
University of Maine, Orono, Maine 
(Received May 25, 1940) 


Data obtained by displacement interferometry have been 
reduced to give the variation of refractive index at infinite 
wave-length with density for pure nitrogen up to 0.01656 
gram per cc. Precise data are also presented for dispersion 
in pure argon, and the Cauchy constants under NTP con- 
ditions are evaluated as follows for this gas: Ao—1 
=0.0002771+0.0000008, By= 1.477 X 10-. It is also noted 
that the Lorentz-Lorenz function, and hence the molar 
refraction at zero frequency, are independent of density for 
these two gases; the density range for argon extending up 
to 0.03342 gram per cc. The values of molar refraction are 
found to be 4.369+0.003 cc for nitrogen and 4.138+0.012 
ce for argon. Preliminary work on carbon dioxide indicates 


a slight negative deviation from linearity for (A —1) with 
density, and a marked negative variation of the Lorentz- 
Lorenz function with density. The latter means a decrease 
in the molar refraction at zero frequency, and suggests a 
decrease in the molar polarization with density. This is 
not in accord with the findings of Keyes and Oncley who 
report an increase of molar polarization with density for 
this gas. The NTP values of the Cauchy constants for 
CO, are also evaluated as follows: (A4o—1)=0.0004419 
+0.0000036, and Byo=2.791 X10“. The extrapolated value 
of the molar refraction at zero frequency is found to be 
approximately 6.65 cc. 


INTRODUCTION 


A PRECISE method for measuring the refrac- 
tive index of a gas by displacement inter- 
ferometry has been described by the writer in 
earlier papers." ? The present paper deals first 
with an analysis of data previously reported for 
purified nitrogen over a range of pressures from 
zero up to fourteen atmospheres.? The reduced 
data give the variation of refractive index at 
infinite wave-length with density, so that it is 
now possible to report with some degree of pre- 
cision the behavior of the Lorentz-Lorenz func- 
tion up to a density of 0.01656 gram per cc and 
to evaluate the molar refraction at zero fre- 
quency, or the molar polarization for this gas. 
The present paper also deals with a similar 
analysis for pure argon over a pressure range from 
zero to eighteen atmospheres, and presents some 
preliminary results for carbon dioxide observed 
over approximately the same pressure range. 
This work was interrupted by the transfer of 
the apparatus to a new location,’ which necessi- 
tated its complete reassembly with numerous 
alterations. The publication of earlier data was 
postponed until they could be confirmed in the 
new location. It is now gratifying to report that 
recent observations taken with the reassembled 
~+C. E. Bennett, Phys. Rev. 37, 263 (1931). 
*C, E. Bennett, Phys. Rev. 45, 200 (1934). 


‘Former work done at the Massachusetts Institute of 
Technology. 


equipment are in complete agreement with the 
earlier data, and hence lend added significance 
to them. The complete analysis of the new CO, 
data, however, in which the density range is 
extended appreciably, will be presented later. 


THEORY 


If the refractive index of a gas is represented 
by a Cauchy expansion 


(1) 


where A presumably represents the infinite wave- 
length value of the refractive index, the funda- 
mental equation of the displacement interferom- 
eter is 


(2) 


where Ax is the observable displacement of the 
fringe system, é is the length of the gas tube, and 
the other symbols have the same significance as 
in the Cauchy equation. It was shown in the 
earlier work that A —1 is linear with pressure for 
pure nitrogen and that only two constants A and 
B are necessary to describe the data up to the 
limiting pressures used, namely, 14.145 atmos- 
pheres. It was also shown that these two con- 
stants, reduced to NTP values and called Ao 
and Bo, respectively, are as follows: for nitrogen 
Ao—1=0.0002932 to better than one-tenth of 
one percent, and By= 1.637 X10-™. 
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Since, to a first approximation, pressure is 
proportional to density and A—1 is propor- 
tional to (A?—1)/(A?+2), it was concluded 
that the above results constitute a rough check 
on the constancy of the Lorentz-Lorenz relation 
(A?—1)/(A?+2)(1/p) in this pressure range, 
where p is the density. Also, since the refractive 
index measurements here reported are very 
accurately determined, the question of convert- 
ing pressures to densities by an accurate equation 
of state is suggested to test the validity of this 
relation more precisely. 

The density values reported in this paper are 
obtained by the use of the Keyes equation of 
state,‘ the constants for which are now available 
for several gases. This equation is usually ex- 
pressed in terms of the pressure as follows: 


RT A (3) 


but it can be rearranged to express the specific 
volume (the reciprocal of density), in terms of 
pressure and temperature, thus 


TABLE I. Refractive index at zero frequency and the Lorents- 
Lorenz relation as functions of density for pure nitrogen. 


1 At-1 
p(g/cc)  p(cemHg) A-1 p Att2 


0.0043061 309.75 0.0010127 0.15676 

0.0010781 0.15645 
0.0011541 0.15662 
0.0012384 0.15633 
0.0012333 0.15538 
0.0016812 0.15528 
0.0017282 0.15608 
0.0019032 0.15650 
0.0019664 0.15591 
0.0019692 0.15612 
0.0021146 0.15546 
0.0021770 0.15612 
0.0025535 0.15551 
0.0027770 0.15570 
0.0027915 0.15591 
0.0030521 0.15548 
0.0033389 0.15622 
0.0037327 0.15565 
0.0038584 0.15524 


Mean 0.15596+0.00008 
REFRACTION Ry=4.369+0.003 


0.014219 1022.20 
0.015977 1075.36 
0.016559 1000.57 


o 

= 
Ss 


*F. G. Keyes, Proc. Nat. Acad. Sci. 3, 323 (1917). 


CLARENCE E. BENNETT 


The values of A —1 at each density are readily 
calculated either by use of the precise value of 
Bo properly adjusted for density and treated as 
a correction factor to Ax/e by Eq. (2) rewritten 


thus 
A-—1=Ax/e—3B/™, (5) 


or by a direct least-squares analysis of the indj- 
vidual values of Ax/e at each of several wave- 
lengths. In the latter method (A —1) and B are 
determined simultaneously for each density ob- 
servation. In the present work data for Ax/e were 
available for each of the following wave-lengths: 
5780A, 5461A, and 4359A, and calculations were 
made by both methods. 


RESULTS 


The results for nitrogen are listed in Table I, 
which shows the constancy of the Lorentz- 
Lorenz function up to a density of 0.01656 gram 
per cc. Attention is called to the fact that up to 
this density, at least, there is no observable 
variation of this function. The average value of 
the molar refraction at zero frequency is found 
to be 4.369+0.003 cc. 

Table II contains hitherto unpublished data 
for pure argon. The refractive index and disper- 
sion values are given for two different tempera- 
tures and six different pressures, reduced to 
density values by a Keyes equation of state as 
indicated above for nitrogen. The purity of the 
argon used was better than 99.8 percent as certi- 
fied by the Air Reduction Company. It was also 
passed through the same purifying and drying 
system previously used for nitrogen. 

The linearity of (A —1) and B with density are 
clearly shown, as is also the constancy of the 
Lorentz-Lorenz function up to 0.033418 gram/ce. 
The values of Ao—1 and Bo are found to be 
0.0002771+0.0000008, and 1.477 X10-, respec- 
tively. The last column gives the values of molar 
refraction at zero frequency, or molar polariza- 
tion and shows it to be independent of density 
with an average value of 4.138+0.012 cc. 


PRELIMINARY WoRK ON CARBON DIOXIDE 


A limited amount of data on pure CO, at 50°C 
is presented in Table III. As in the case of argon, 
the values of (A—1) and B were obtained by a 
least-squares analysis of the observations of Ax/e 


- 


At 
| 
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TABLE II. Dispersion in pure argon as a function of density. 
Ax/e 
t p 

No. cm HG °c GRAMS/CM? 5780A 5461A 4359A 
1 305.63 00.00 0.0071834 | 0.0011684 0.0011750 0.0012057 
2 573.20 30.00 0.012143 0.0019766 0.0019880 0.0020507 
3 567.26 00.00 0.013372 0.0021815 0.0021932 0.0022529 
4 1071.72 30.00 0.02275 0.0037168 0.0037359 0.0038450 
5 1093.21 00.00 0.025916 0.0042235 0.0042449 0.0043693 
6 1405.05 00.00 0.033418 0.0054422 0.0054699 0.0056398 * 

A-1 B At—11 

No. A-1 Bx10" X10" At42 Re 
1 0.0011199 5.438 0.15590 756.43 0.10391 4.1506 
2 0.0018786 10.903 0.15470 897.84 0.10310 4.1182 
3 0.0020877 10.465 0.15612 782.57 0.10405 4.1562 
4 0.0035465 18.899 0.15572 829.83 0.10375 4.1442 
5 0.0040295 21.513 0.15548 830.10 0.10359 4.1374 
6 0.0051780 29.2319 0.15495 874.74 0.10321 4.1226 

MEAN 0.15548 828.58 0.10381 4.138+.012 


Ao—1=0.0002771 +0.0000008 


Bo=1.4777 X10™™ 


for each of three wave-lengths. In the case of CO, 
the deviations from the ideal gas law make it 
absolutely necessary for all pressures to be con- 
verted to densities before linearity considerations 
are even attempted. Here again, the Keyes equa- 
tion of state was used. The observations were 
taken with the original apparatus for which the 
tube length e was 100.30 cm and the correction 
factor for the half-millimeter micrometer screw 
was 1.0104. The temperature was maintained 
constant to 0.01°C and pressure measurements 
were made with a Keyes type dead-weight gauge 
to approximately one millimeter of mercury over 
the entire range. The zero pressure reading was 
determined with a McLeod gauge to better than 
10-* mm Hg. 

It is to be noted that although B is very nearly 
linear with density, A —1 deviates slightly. If one 
disregards for the moment this slight negative 
deviation, and considers the average values of 
A-—1, and B, the following values of Ay—1, and 
By, are obtained by the use of standard density 
po=0.0019769: A o—1=0.0004419+0.0000036, 
and By=2.791 X10-". 

Incidentally, values of ~—1 were also de- 
termined for several wave-lengths by Eq. (1), 
and the results were found to be in good agree- 
ment with those generally accepted. They are not 
listed here, however, since they are readily ac- 
cessible in the literature. 

Table III also shows the variation of the 


Lorentz-Lorenz function with density and the 
molar refraction at infinite wave-length for 
carbon dioxide. There appears to be a definite 
increase in molar refraction with density for this 
gas whereas this quantity seemed independent of 
density for nitrogen and argon. The molar re- 
fraction at zero frequency is found by extrapo- 
lation to be approximately 6.65 at zero density. 


CONCLUSIONS 


In view of these results it is concluded that the 
Lorentz-Lorenz function, and hence the molar 
refraction at zero frequency, are independent of 
density for pure nitrogen and argon over this 
density range, although Keyes and Oricley' find a 
slight increase in the molar polarization for 
nitrogen based upon dielectric constant measure- 
ments. The value of the molar refraction at zero 
frequency here obtained for nitrogen, however, 
is in good agreement with their value of molar 
polarization (4.37 cc), thus confirming the va- 
lidity of the Maxwell law over this density range. 

The results for argon represent a contribution 
to the rather limited amount of information 
available concerning the Lorentz-Lorenz function 
at high densities for the rare gases which are 
structurally quite simple. That this function 


5 F. G. Keyes and J. L. Oncley, Chem. Rev. 19, [3] 195 
(1936). 
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TABLE III. Dispersion in CO: as a function of density (data at 50°C) 


MICROMETER READING Ax/e 

No. asus 5780A 5461A 4359A 5780A 5461A 4359A 
0 0.0 00.5105 00.5790 00.9690 

1 354.96 04.2500 04.3390 04.8500 0.0018835 0.0018939 0.0019548 
2 426.61 05.0090 05.1050 05.6495 0.0022658 0.0022797 0.0023575 
3 514.08 05.9565 06.0585 06.6285 0.0027431 0.0027600 0.0028506 
4 612.54 07.0260 07.1345 07.7425 0.0032818 0.0033019 0.0034117 
5 724.10 08.2410 08.3575 09.0180 0.0038938 0.0039180 0.0040542 
6 846.31 09.5950 09.7260 10.4130 0.0045758 0.0046072 0.0047568 
7 982.21 11.1395 11.2700 12.0415 0.0053537 0.0053849 0.0055771 
8 1210.00 13.7780 13.9305 14.8160 0.0066827 0.0067250 0.0069746 
9 1575.90 18.1455 18.3165 —— 0.0088823 0.0089342 — 

B At—1)\1 
No. c/ec ans Bx10" (4053) 5 Ro 
1 0.0079023 0.0017886 10.529 0.22634 1332.4 0.15084 6.638 
2 0.0095386 0.0021442 13.504 0.22479 1415.7 0.14981 6.593 
3. 0.011557 0.0026011 15.806 0.22507 1367.7 0.14998 6.601 
4 0.013857 0.0031100 19.105 0.22443 1378.7 0.14955 6.582 
5 0.016503 0.0036813 23.635 0.22307 1432.2 0.14862 6.541 
6 0.019452 0.0043398 26.435 0.22310 1359.0 0.14863 6.541 
7 0.022798 0.0050419 34.177 0.22115 1499.1 0.14731 6.483 
8 0.028580 0.0062938 43.100 0.22021 1508.0 0.14665 6.454 
9 0.038381 —_ — — 
Bo=2.791 X10" 


(Ao— 1) =0.0004419 +0.0000036 


would be constant for a monatomic gas is prob- 
ably to be expected. 

The negative variation of the Lorentz-Lorenz 
function with density for carbon dioxide, on the 
other hand, is also contrary to the findings of 
Keyes and Oncley who find a positive variation 
of the molar polarization. This point, however, is 
to be investigated further. Phillips® has reported 
such a negative variation although not so marked, 
but it is doubtful if the Phillips data are as ac- 
curate as that reported here over this density 
range. It may be that molar refraction at zero 
frequency is not the same as the molar polariza- 
tion for carbon dioxide although Keyes and 
Oncley® show that the Phillips data, corrected for 
the effects of infra-red absorption bands, as sug- 


¢P. Phillips, Proc. Roy. Soc. London 97A, 225 (1920). 


gested by Fuch,’ yield a value for molar refrac- 
tion nearly in agreement with the molar polariza- 
tion. The agreement would be even better, how- 
ever, on the basis of the results reported here. 
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The rate of increase of the entropy of a simple viscous fluid, capable of conducting heat, is 
investigated in greater detail than has hitherto been customary. It is shown that, if Kelvin’s 
hypothesis concerning the absolute temperature be adopted, and the usual law of viscosity be 
assumed, the requirements of the second law are satisfied. 


INCE the publication of Gibbs’ ‘Equilibrium 

of heterogeneous substances,’"' it has been 
obvious that the greater part of the existing 
thermodynamic theory would more appropriately 
becalled thermostatics, the term thermodynamics 
being reserved for the rather unorganized lot of 
equations that govern irreversible changes in 
which heat plays a part. Examples of such 
equations are Fourier’s law of heat conduction, 
Ohm’s law, Fick’s law of diffusion, etc., etc. 
Although this usage is dictated by the close 


analogy to mechanics, there has been no tendency — 


to adopt it. This reluctance has been partly due 
to the chaotic state of thermodynamics proper, 
and partly to the belief that the laws mentioned 
above are much less fundamental than Newton's 
laws of motion or than the energy and entropy 
laws, and require derivation from kinetic theory. 
The purpose of the present investigation has 
been to find the general equation of which the 
foregoing are special cases. The problem of 
providing a kinetic theory derivation has not 
been considered. The major part of the discussion 
will be presented in such a manner as to bring 
out the connection with Gibbs’ work, but by 
way of introduction, the case of a simple viscous 
fluid, conducting heat, will be discussed from a 
pre-Gibbsian standpoint. 

Let m be the mass density of the fluid, V its 
velocity, so that the conservation of mass leads 
to the equation of continuity 


am/at+V- (mV) =0. (1) 
Defining the specific volume v=1/m and the 
operator 
D/Dt=0/dt+V-V 
4 f. W. Gibbs, Scientific Papers (Longmans, 1906), Vol. I., 
p. 55. 


267 


we may also write 
mDv/Dt=V-V. (2) 


Let S be any closed surface moving with the 
fluid, so that an element of its area, de, sweeps 
out the volume V-dedt in time dt. If f is any 
function of position and time 


ca/at) f mar 
f f J [a(mf) /at]dr-+ f mfV-de 


f f J fa(mf) /at-+-V-(mfV)\dr 


= m(Df/Dt)dr. (3) 


In the absence of external volume forces, the 
motion of the fluid will be governed by the 
hydrodynamic equation 


mDV/Dt= —Vp+V-», (4) 


where p is the hydrostatic pressure, and p is the 
tensor of the viscous stresses. The precise ex- 
pression for the latter will not be needed until 
later. 

The kinetic energy of the fluid enclosed by S is 


k= f ff amv 


and by Eq. (3) 


dK /dt= f f f mV -(DV/Dt)dr 
= J J 


| 


or, on using the divergence theorem 


aK /at=— (V—p-V)-de 


+f vier (5) 


The surface integral is the rate at which the 
fluid inside S does work on the fluid outside S. 
The energy law as applied to this simple 


fluid, states that there is a function w(p, v) (the 


internal energy per gram) and a vector q (the 
heat flow) such that 


(d/dt) f f (4m V24+-mw}dr 
f f f (pV—p-V)-de. (6) 


The first integral on the right is the rate at which 
heat enters S and the other has already been 
identified. Subtracting Eq. (5) from Eq. (6) 
gives 


f f f {m(Dw/Dt) +p¥-V—(p-¥)-V 
+V-q}dr=0, (7) 


and since this must be true for all closed surfaces, 
it follows that the integrand itself must vanish. 

Since w is a function of p and v only, and 
because of Eq. (2), the first two terms may be 
transformed : 


m(Dw/Dt)+pV-V 
=m{[(dw/dv) +p ](Dv/Dt) (8) 
+ (dw/dp)(Dp/Di)}. 


Now, there are always two functions @ and 7 
such that 


(dw/dv)+p=0(dn/dv), (dw/dp)=0(dn/dp), (9) 
so that finally 
m(Dw/Di)+pV- V=m@(Dn/Dt), 


whence 
(10) 


Before making use of this equation, it is worth 
while to consider this definition of @ and 7, 
despite the fact that it is discussed in all standard 
texts on thermostatics. Differentiating the first 
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of Eq. (9) with respect to p and the second with 
respect to v, and subtracting gives: 


(00/dp) (dn/dv) — (30/dv) 0p) =1. 
Elimination of » by means of Eq. (9) results in 


(8 log 6/dp)[(dw/dv) ] 
— (0 log @/dv)(dw/dp)=1, (11) 


which can be solved for log 6, since w is supposed 
known. If 6=40(p, v) is one solution of Eq. (11) 
and ¢(~, v) is any solution of 


(8¢/dp)[(dw/dv) +p ]—(d¢/dv)(dw/dp) =0, (12) 
then the most general solution of Eq. (11) is 

0= 6of(9), (13) 
where f is arbitrary. To complete the definition 


of @, this function must be determined in 


some way. 
Thereafter, the function 7 is determined except 


for an additive constant. From Eq. (9) 


n(P, = Vo) 
(1/6) {C(aw/a0) +p 


Povo 

and Eq. (11) is the condition that the integral 
be independent of the path. 

The definition of 6 is completed on the basis 
of what may be called the Kelvin hypothesis: 

The arbitrary function f may be chosen so that 
the heat flow q is related to 0 by Fourier’s equation 

q= —kV8, (14) 

where k is a function of p and v that does not 
change sign. 

When f is so chosen, @ is unique except for 
sign and unit. For, let 60 also satisfy the equation 


q=- 
then it follows that 
0= F(60) 
or 
= F(80)/60. 


But, because of Eqs. (11) and (12), ¢ cannot be 
a proper function of 6, so that both sides of 


the equality must be constants; hence 


6=ao. 


| 
4 
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The constant @ can always be chosen so that k the second law of thermodynamics. It remains to 


. is positive. When thus defined, @ is called the consider the case of a viscous fluid ; the commonly 
absolute temperature and 7 the entropy (per accepted expressions for the components of ) are? 
gram) of the substance. Dee = (2A/3)[2(8 —(8V,/ay) —(aV,./8 

Returning to Eq. (10), we may write it v/dy) — 
m(Dn/Dt)+V- (q/6@ . >0); 

o = —[q-¥0]/e-+[(p-¥)-V1/e where A( 2 0) is the coefficient of viscosity. Hence 

. (3) and the diverge ce theorem (p- V)- V=(2d/3){ V,/dx) V,/dy) P+ coe} 

) (d/dt) f f f mndt+ f f (1/0)q-do If, now @>0, it follows that the right side of 
s 8 Eq. (15) is never negative even for a viscous 
2) fluid. If this be accepted as an empirical fact, or 
= J f f {[k(V6)*]/6?+[(p-V)-V]/@}dr. (15) included in Kelvin’s hypothesis, the latter and 

” Eq. (15) have the inequality 

3) The first term on the left is the rate of increase 

oi of the entropy of the matter inside S; the surface ca/at) fmnar+ f (16) 

_ integral is the familiar ratio of heat loss to tem- 

perature. Thus, if VWW=0 and Vé=0, this equa- 
tion expresses the ordinary relation between 4S Consequence. As has been remarked, this 

4 entropy, temperature and heat in reversible imequality is a partial expression of the second 
processes. ’ law of thermodynamics. 

Furthermore, if p=0 (nonviscous fluid), the ?H. Lamb, Hydrod ics (Cambridge. 1924). fifth edi 
’ ’ ’ ed 
right side is never negative—as is asserted by tion, p. $44. 
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The Thermodynamics of Irreversible Processes 

II. Fluid Mixtures 
at Cart Eckart, 
on Ryerson Laboratory, The University of Chicago, Chicago, Illinois 
4) (Received April 19, 1940) 
ot The possibility of constructing a systematic theory of irreversible processes is surveyed in 
general terms, by utilizing some of the results established in later parts of the paper. Three 

or assumptions underlying Gibbs’ application of the second law to equilibrium problems are 
formulated in explicit but general mathematical form. It is shown that they restrict the equa- 

mn tions governing irreversible changes. The theory of a general fluid mixture is developed in some 
detail, and is then applied to mixtures of ideal gases. It is shown that the usual equations for 
the velocity of chemical reactions are consistent with the second law provided that the de- 
parture from equilibrium is not too great. Mathematical complexities make it difficult to de- 
cide whether this is the case for larger deviations also. A somewhat general theory of diffusion 
and heat flow is considered and the requirements of the second law are formulated as the positive 
definiteness of a certain matrix whose elements depend on the diffusion coefficients, thermal 
conductivity, etc. 

De SURVEY thermodynamics. Before proceeding to more 

of 


N the preceding paper, the rate of increase of complicated cases, it is well to consider what may 
the entropy of a simple fluid was calculated be accomplished. 
and found to be in accord with the second law of In the case of the simple fluid, it was possible 


| 

1 

| 
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to establish the existence of entropy and tem- 
perature functions. This depended on mathe- 
matical properties that are peculiar to Pfaffian 
forms in two variables (p, v), and it is not to be 
expected that the discussion can be extended to 
the case of more variables.’ It will therefore be 
advisable to follow Gibbs and suppose that the 
concepts of temperature and entropy are not in 
need of definition or existence proof. 

It was also possible to deduce the usual incom- 
plete form of the second law from Kelvin’s 
hypothesis, the positiveness of the absolute tem- 
perature, and the usual expression for the viscous 
stresses. This may seem of doubtful interest, 
since these hypotheses are all less general than 
the second law itself. Thus the object of the con- 
siderations will not be the derivation of the 
entropy law from other laws. 

It will rather be to derive equations of the 
form? 


=f (1) 


which, when combined with the second law, 
yields the inequality 


G20. (2) 


At this stage of the considerations, it will be 
shown that G is always the sum of a number of 
terms, each of which is the product of two 
factors: 


G= (3) 


This decomposition of G, and more particularly, 
the classification of the factors into X’s and Y’s, 
is to a great extent arbitrary. However, it will 
then be shown that it can be done so that 


1 In this connection, the work of C. Carathéodory should 
be mentioned: Math. Ann. 67, 355 (1909), Berl. Ber. p. 39 
(1925); also T. Ehrenfest Afanassjewa, Zeits. f. Physik 33, 
933 (1925), and S. Chandrasekhar, Introduction to the 
Study of Stellar Structure (Chicago, 1939), Chap. I. 

There is reason to believe that Carathéodory’s theory 
affords a better basis for a theory of irreversible processes 
than does that of Gibbs. The author has chosen the latter 
because of its greater familiarity. 

? The notation is the same as that of the previous paper; 
m, n, 8, q are the mass density, entropy per gram, tempera- 
ture and heat flow, respectively. S is any closed surface 
moving with the matter. 


270 CARL ECKART 


A. the definition of thermostatic equilibrium is 
X,=0 
together with the vector equation 
dV/at=0. 


For example, the X’s are the components of the 
heat flow q, of the velocity gradients, of the diffy- 
sion velocity, the rates of production of chemical 
compounds, etc., etc. 

The Y’s are determined by this choice of the 
X’s, and it will then be shown that 
B. the Gibbsian criteria of thermostatic equilibrium 
are equivalent to the equations 


Y,.=0 
together with the vector equation 
Vp=0. 


The Y’s corresponding to the X’s enumerated 
above are proportional to the components of the 
temperature gradient, of a stress tensor, of the 
gradients of chemical potential, and the differ- 
ences of chemical potential, in the same order as 
above. 

It is desirable to have an inclusive terminology 
for the X’s and Y's; because of the diverse na- 
tures of the quantities thus combined into 
classes, it is difficult to find a terminology that 
is free from serious objections. The most satis- 
factory that has occurred to the author is D- 
factor for the X’s, (D for definition) and C-factor 
for the Y’s (C for criterion). The components of 
the acceleration and pressure gradient should be 
considered D- and C-factors, respectively. 

It remains to give an explanation of these 
observations. This is to be found in tacit as- 
sumptions that enter into all discussions of 
thermodynamics and thermostatics, but are 
rarely or never formulated clearly : 

a. The irreversible changes of thermodynamic 
state are governed by equations expressing the 
D-factors as functions of the C-factors and certain 
other variables (say &): 


F,(Y, é). (4) 


These may be called the thermodynamic equa- 
tions of motion; it should be noted that the C- 
and D-factors have been defined without refer- 
ence to the equations of motion. The definition 
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of thermostatic equilibrium then leads immedi- 
ately to the equations 


F,( y, =0, (5) 


but it must be assumed that 

b. The Eq. (5) are equivalent to the equations 
under B, i.e., to the vanishing of the C-factors. 

c. The inequality 


(6) 


is a consequence of the thermodynamic equations of 
motion. 

Unless a, b, and ¢ are true, it is not true that 
“For the equilibrium of any isolated system it is 
necessary and sufficient that in all possible varia- 
tions of the state of the system which do not 
alter its energy, the variation of its entropy shall 
either vanish or be negative’’*® which is the prin- 
ciple on which Gibbs based his theory of thermo- 
statics. In order that this principle follow from 
a, b, c, it is further necessary that B be true. 

In the next section, the results that have been 
anticipated in this survey will be established for 
the case of a fluid mixture. 


FLuip MIXTURES 


For the sake of formal simplicity, it is essential 
to consider the possibility that the mixture may 
contain all of its chemical elements in the free 
monatomic form; let cz, a=1, 2, --- be the con- 
centrations of the free elements in moles/gram 
of mixture. Let the other possible components of 
the mixture be molecular compounds, the ith 
containing vj. atoms of element a in each of its 
molecules. For simplicity, the numbering may be 
such that 7 runs from n+1 to N, and let c; be the 
concentration of compound 7. The word substance 
will be used to designate either free elements or 
compounds. 

Let V; be the velocity of the substance k and 
m the density of the mixture in g/cm’; then 


O(mc,) =mT;, (7) 


defines I',, the rate of production of substance k 
in moles/g/sec. Since the atoms of the elements 
are not created or destroyed by the chemical 
reactions, the I’, are not independent, but re- 


+ J. W. Gibbs, Scientific Papers (Longmans Green & Co., 
New York), Vol. I, p. 56. 
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lated by 
Vials =(, a=1---n. (8) 


If M;, is the molecular (atomic) weight of 
substance k, the concentrations are related by 


Mice =1, (9) 
and the M, by 
M;= VieMa, t=n+1, (10) 


Equations (8) and (10) combine to give the 
equation of conservation of mass 


=0. (11) 


The momentum density of the mixture is mV, 
where 


V= Mick Vi, (12) 


and the Eqs. (7), (9), (11), (12) combine to give 
the equation of continuity : 


dm/dt+V-(mV) =0 
or (13) 
mDv/Dt=V-V. 


The diffusion velocities of the various substances 
are 


U.=V.-V, (14) 
and Eggs. (7), (13), (14) combine to give 
(15) 


The hydrodynamic equation will be assumed 
in the form 
mDV/Dt=V-», (16) 


the stress tensor p now including the hydrostatic 
pressure as well as the viscous stresses. There is 
now the possibility that p may also depend on the 
U;,,, but this complication may be ignored. In any 
case, the kinetic energy equation in the form 


famvrar= ff (17) 


will still be valid. Following Gibbs, it will be 
assumed that the internal energy (per gram of 
mixture), €, is a function only of v, of the c,, and 
of the entropy 7. The possibility that « may de- 
pend on the U; will be ignored, but this consti- 
tutes an approximation. 

Because of the arbitrary zero of energy, and 
of Eq. (9), all equations must be invariant when 
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« is replaced by e+f(2:Micx), where f is arbi- 
trary. This will be called an e-substitution.‘ 

An application of the energy law to the volume 
whose boundary S moves with the mass of the 
fluid yields 


cavat) f f V24meldr 
f 


ff (18) 


The first integral on the right is the rate at which 
heat enters the volume; the second is the rate 
at which the fluid outside S does work on the 
fluid inside S; the remainder is the rate at which 
energy enters the volume due to the diffusion 
processes. It is readily verified that this equation 
is invariant under an e-substitution. 
Subtracting Eq. (17) from (18), transforming 
the surface integrals into volume integrals, and 
recalling that S is arbitrary, one obtains: 


(19) 


Following Gibbs, the temperature is defined by 
6=de/dn, (20) 
the chemical potentials® by 
(21) 
and the hydrostatic pressure by 
p= —de/dv. (22) 


(Under an e-substitution, uw, is replaced by 
M;,f’(1), so that Vux is invariant.) Then 


De/Dt=0Dn/Dt— Dt, 


‘ It will be noted that the theory of ¢-substitutions is not 
used to prove any equation, but only as a check and guide. 
Because of the equivalence of mass and energy, the status 
of the e-substitutions in a fundamental theory is not clear, 
and it has seemed preferable to restrict the use of these 
methods, Other types of invariance are also associated with 
Eq. (9) and can “eager be used in a similar manner. 

Gibbs calls the uw; simply “the potential of the sub- 
stance 4.” Other writers call them “partial molal thermo- 


dynamic potentials.” 
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which combines with Eqs. (13), (15), and (19) 
into 


m {xT e+ =0. (23) 

From Eq. (8), it follows that 
= MAT, (24) 


where 
(25) 


is the difference in chemical potential between 
the compound 7 and its equivalent in free ele- 
ments. Because of Eq. (10) the A; and therefore 
also the Eq. (23) are invariant under an ¢- 
substitution. 

Dividing Eq. (23) through by @, integrating 
the result over the volume bounded by S, and 
comparing with Eq. (1), show clearly that 


The derivation of Eq. (26) completes the first 
stage of the program outlined in the survey (cf. 
Eq. (3)). It remains to verify the propositions 
A and B. 

It will be obvious that when the mixture 
reaches thermodynamic equilibrium, the follow- 
ing quantities will vanish: (1) the heat flow, q; 
(2) the velocity gradient, VV (and therefore also 
V-V); (3) the rates of production of the various 
substances, mI';; (4) the diffusion currents, 
mc,U,.; (5) the acceleration, V/dt. The first four 
of these quantities appear as factors in the vari- 
ous terms of G, thus verifying A. The correspond- 
ing C-factors are 

(1) (1/8); (3) —(1/0)A;; 

(2)§ (1/6)(pi+p); (4) —(1/6)Vun; 
while the Gibbs criteria for thermostatic equi- 
librium are 
(1) Va=0; (2) p=—pi; 

(4) Vue=0; 
thus verifying B. 
However several remarks are called for in this 


connection: the second criterion is not found 
explicitly in Gibbs’ paper, but is implicitly as- 


(3) A:=0; 
(S) Vp=0; 


* Here i is the unit tensor. 
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sumed at an early stage of the calculations. The 
criterion (5) follows from criterion (2), the hydro- 
dynamic equation (Eq. (16)), and aV/at=0. The 
exceptional character of the fifth criterion is thus 
seen to be related to the exceptional position 
accorded the hydrodynamic equation as com- 
pared to the other equations of motion. 


ELIMINATION OF THE ENTROPY FROM 
THE EQUATIONS 


The equations just derived involve functions 
of v, the c;, and the entropy, 7. This proves to 
be inconvenient, for it is an empirical fact that 
these functions are very complicated. The equa- 
tions are very much simplified by eliminating the 
entropy and replacing it by the temperature, 
using Eq. (20) to accomplish this. This is an 
algebraic process, but can be very laborious, so 
that it deserves some consideration. 

It appears that the process can be much 
simplified by introducing the function 


Ci, 6) = e(v, Ci, n) 6n, (27) 


the entropy being eliminated from the right side. 
The function y is one of those introduced by 
Massieu under the name characteristic function, 
and is usually called the free energy, or work 


_ function (per g of the mixture). Once y has been 


determined, the quantities p, uw; and 7 can be 
expressed as functions of v, c; and @ by means of 
the equations 


p= —dy/dv, (28) 

(29) 

n= —dy/d8, (30) 

which follow from Eqs. (20), (21), (22) and (27). 


GENERAL CONSEQUENCES OF a, ¢ 


The assumptions a, b, c have important conse- 
quences for the theory of irreversible processes, 
inasmuch as b, and c impose marked restrictions 
on the equations of motion. From 6 it follows 
that 

(0F./dt)o=0, (31) 


the subscript zero indicating that the Y,=0 
while the ¢’s have general values.’ Equation (6) 
‘If the quantities dX. o/dt are the é’s, this 


ality must be restricted, since these 
equilibrium. 


erivatives at 
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imposes further restrictions on the F,; it can be 
given another form if the F, are linear functions 
of the Y,, which is probably the case whenever 
the departure from equilibrium is not very 
great. In this case, Eq. (4) becomes 


Y,, (32) 
where 
Aw= (0F,/d Y5)o, (33) 


is in general a function of the £’s. Equation (6) is 
then 
0, (34) 


which must be true for all values of the Y, and 
the £’s. Therefore the matrix 


ba) (35) 


must be positive definite. This implies such 
inequalities as 
Bua? 0, (36) 


BaaBos— (Bar)? 2 0. (37) 


MIXTURES OF IDEAL GASES 


In order to work out special cases of the general 
principles, it is convenient to consider a mixture 
of ideal gases, for which the work function is 


Cis 6) 
= ci {ui(0) +RO[log (ci/v)—1]}, (38) 


the u; being functions characteristic of the sub- 
stances and R the gas constant per mole.® 
Equation (29) then yields 


log (c:/v), (39) 


and Dalton’s law of partial pressures, etc., can 
also be derived. 


This ideal mixture may be the scene of’ 


chemical reactions, proceeding at rates given by 
the usual formulae. It is to be considered, whether 
these formulae are consistent with a, b, c. Before 
doing so, however, it will be shown that Eq. (39) 
leads to the usual form of the law of mass 
action. 

Let M;, stand for the chemical symbol of 
substance k, as well as for its molecular weight. 
Then the general form of a reaction may be 


See, F. Herzfeld, Handbuch der Physik, Vol. 9, 
p. 101, Eq. ( 
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written (the »’s are integers) 
ve! vy!’ My. (40) 
Since the reaction balances, 


Veq being the number of atoms of element a@ in 
one molecule of substance k. (If a{n, via=Sxa-) 
The criteria for chemical equilibrium (4;=0) 
yield 


” 
— 


which may also be written 
vi’) mi =0 (42) 


because of Eq. (41). Substituting from Eq. (39) 
gives 
— vi’) [u;/RO+log (c;/v)]=0, (43) 


which is essentially the Arrhenius form of the 
law of mass action. 

The commonly accepted (and empirically veri- 
fied) formulae for the rate of a chemical reaction 
in the gaseous phase were first deduced by 
Guldberg and Waage from rather vague kinetic 
theory considerations. These have since been 
made somewhat more specific and brought into 
agreement with Eq. (42).° According to these 
formulae, the net rate at which the reaction 
proceeds to the right is the difference p’—p”, 
where 

=oexp 


@ being a positive quantity. It is clear that 
Eq. (42) follows from this and 


p’—p”=0, (45) 


but that Eq. (45) does not follow from the 
definition of thermostatic equilibrium. It is 
usually deduced from the principle of detailed 
balancing, but also follows from the Gibbsian 
criteria of equilibria. 

If there are many possible reactions, these 
must be distinguished by an index (¢,, ps’, Mas, 
etc.). The rate of production of substance 7 is 
then 

= vie’ Vis ) ps"), (46) 


*R. H. Fowler, Statistical Mechanics eee Uni- 
versity Press, 1936) second edition, p. 700. 


and the definition of equilibrium is T;=0. Eq, 
(45) cannot follow directly from this and Eq. (46) 
unless there are at least as many substances as 
possible reactions. It may, however, follow in- 
directly because of other considerations leading 
to the Gibbsian criteria of equilibrium. 

From Eq. (25), 


Ps EXP res’ RO), 
fe" =a; exp ves Ré), 


EXP RO) > 0. 


(47) 


where 


Since p and A are invariant under an e-substitu- 
tion, a must also be; there is thus a presumption 
that the interpretation of a will be simpler than 
that of ¢. Equation (46) then becomes 


—exp Vee’ Ax/RO)}. (48) 


For sufficiently small values of Aj, 
i= vis” — vis’ — Ves’ )Ar/RO, 
so that 
Vieng —Ai/8) 
=(1/R)Zo, — ven’) (Ae 


These last two equations are in complete accord 
with the general considerations leading to Eq. 
(34). The more general Eq. (48) is in accord with 
Eq. (31), and the Gibbsian criteria A4;=0 are a 
sufficient condition for T;=0. It is not easy to 
see whether they are also the necessary condition, 
as Gibbs asserts, nor whether 


Ti(—A,/0) 20 


for all values of A;, as is required by the present 


assumption c. 
Turning now to the conduction of heat and 


the diffusion of matter in the mixture, it may be 
supposed that the equations governing these 
processes are 


q = kooV (1/0) — (Ko: /0) Vu, 
me ;= xjoV (1/0) — /0)V ui. 


(49) 


In order that these equations be consistent with 
the general principles, it is necessary and 
sufficient that the matrix }3(«;;+«,;:) be positive 
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definite. Insofar as the mixture is an idealized 
one, there is no need to justify the assumption 
of Eq. (49), but the question of its departure 
from actual mixtures does arise. 

The conduction of heat is usually supposed to 
be given by Fourier’s law that q is proportional 
to V6; similarly, Fick’s law that U; is proportional 
to V(mc;) is supposed to govern diffusion. How- 
ever, more general cases have been considered, 
in which U; is a linear function of all the concen- 
tration gradients. Such “‘diffusion-drag” forces 
have been invoked to explain certain biophysical 
phenomena, and it seems likely that they are 
important in the separation of isotopes. Except 
for a difference in notation, and the terms in 
koi, the Eq. (49) are identical with those con- 
sidered in the diffusion-drag theories. 

To see this, note that because of Eq. (39) and 
m=1/v 


(50) 


where 
Ga;=du;/d0+R log (mc;), 


so that Eq. (49) becomes 


—q koiV(mc;), 
($1) 
with 
jo = 0? + 


(52) 
kj; = Rxji/me;j. 


The coefficient kop is the thermal conductivity, 
the kj; (¢=1---N) are the diffusion coefficients, 
while the k;; ({4j=1---N) are the drag coeffi- 
cients. If the coefficients kj94#0, this amounts to 
assuming that a temperature gradient tends to 
produce a diffusion of matter; similarly, if ko;#0 
the concentration gradients tend to produce a 
flow of heat. Both of these phenomena occur; 
they are closely related to the various thermo- 
electric effects. 

Since the matrix x;; is positive definite, the 
thermal conductivity and diffusion coefficients 


must be positive quantities, and 
jh iik (Rises +k 


etc. However, these particular inequalities are 
based on Eq. (40) ; they may have different forms 
for liquid mixtures. The general conclusion that 
the numerical values of the drag coefficients are 
restricted by the numerical values of the diffusion 
coefficients is valid in all cases. 

Finally, it would be possible to consider 
generalizations of the usual formulae for the 
viscous stresses, but this seems of slight interest 
here. 


SUMMARY 


The restrictions imposed by the second law of 
thermodynamics on the equations governing 
irreversible changes have been investigated, first 
somewhat generally, and then for the particular 
case of a mixture of ideal gases. It was found 
that they are automatically fulfilled by the usual 
law of mass action for chemical reaction veloci- 
ties. In the case of general diffusion processes, 
the diffusion coefficients must be positive, and 
certain other inequalities must be satisfied. The 
general methods used can be applied to other 
irreversible processes, and will always yield such 
inequalities. The experience of the author sug- 
gests that many of these inequalities have already 
found their place in physical theories, having 
been introduced on the basis of “physical 
intuition,”’ but that the latter cannot deal with 
the more elaborate cases. 


Note added April 26, 1940. Since writing the above, my 
attention has been called to the work of E. Lohr on similar 
problems.* A first reading of his extensive papers shows 
that some of his conclusions practically coincide with mine, 
though there seem to be appreciable differences in our 
points of view. The entropy law has also been formulated 
as a differential equation by R. C. Tolman,f, independently 
of the work of Lohr. 


*E. Lohr, Wiener Denkschriften, 93, 339-421 (1917); 
99, 11-37, 59-91 (1924). Vektor- und Dyadenrechnung, 
p. 312 (Berlin, 1939). 

} Literature references in R. C. Tolman and H. P. 
Robertson, Phys. Rev. 43, 564 (1933). 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may be secured by 
addressing them to this department. Closing dates for this department are, for the first issue of the 
month, the eighteenth of the preceding month, for the second issue, the third of the month. Because of 
the late closing dates for the section no proof can be shown to authors. The Board of Editors does 
not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


Induced Radioactivity of Rhenium and Tungsten 


Rhenium, with two stable isotopes, mass numbers 185 
and 187,! has been made radioactive and one new activity, 
with a half-life of 5242 days, found in addition to the 
18-hour? and the 90-hour*® periods reported previously. 
The probable assignment of these activities to their 
isotopes is shown in the diagram below. 


184 185 186 187 188 
52days stable 90hr. stable 20hr. 


Experimental results supporting these conclusions 
follow. Slow neutron and deuteron bombardments, fol- 
lowed by chemical separations, produce two activities in 
rhenium, whose half-lives are 18 hours and 90 hours. Fast 
neutrons yield three activities, chemically identified as 
rhenium, with half-lives of 18 hours, 90 hours, and 50+5 
days. From the relative intensities of the 18-hour and 90- 
hour periods produced by slow and fast neutrons, Yamasaki 
and Sinma’s‘ assignment of the 18-hour period to Re'** 
and the 90-hour period to Re!**, has been confirmed. 

The rhenium separation from the bombardment of 
tungsten with deuterons has given three activities with 
half-lives of 20 min., 90 hours, and 52 days. Production 
of the 90-hour activity in this way, no 18-hour activity 
being present, is an additional confirmation of the assign- 
ment of the former to Re!**, and indicates a d,2n reaction 
from W'**, Deuteron capture would probably give both 
the 18-hour and 90-hour activities from the equally 
abundant W'* and W'® isotopes, or if the placement were 
reversed and the d,2n reaction alone occurred, only the 
18-hour activity would be produced. 

The 52-day activity, which was identified as a rhenium 
isotope through a series of chemical reactions, has been 
assigned to mass number 184 for the following reasons. It 
is produced from rhenium by fast neutrons (m,2n reaction 
from Re!®), and appears in rhenium, chemically separated 
from tungsten bombarded with deuterons (probably a d,n 
reaction from W!*, although a d,2n reaction from W!™ is 
an additional possibility). It is a negatron emitter, having 
a gamma-ray of 0.85+0.1 Mev associated with it. 

The 20-min. activity, a positron emitter, follows the 
chemistry of rhenium but is probably due to Mn® formed 
from iron, an impurity in the tungsten foil, since the in- 
tensity differed in different tungsten samples. 


Investigations, during the past year, upon the radio- 
activity produced in tungsten, by deuterons, fast and 
slow neutrons, have given two activities, 24.1+0.1 hours 
and 74.5+1.5 days. The 24-hour period is identical with 
that found by Fermi,? while the 74.5-day period cor- 
responds to that of 77+3 days, just reported by Minakawa.s 

Both activities are negatron emitters, the upper limit 
of the 8-ray energy for the 24-hour period being 1.40-1,44 
Mev by absorption in Al, and 1.40+0.05 Mev from pre- 
liminary cloud-chamber measurements. A y-ray of 0.87 
+0.03 Mev, as measured by absorption in Pb, is associated 
with the 24-hour activity. The upper limit of the 8-ray 
energy of the 74.5-day period is 0.55-0.65 Mev, from 
absorption measurements in Al, and 0.64-0.72 Mev from 
preliminary cloud-chamber measurements. 

The results of extensive investigations of the relative 
intensities of the two periods, produced by fast and slow 
neutrons, does not give as clear a picture as would be 
expected from the assignment, made by Minakawa, 
of the 24-hour to W'®’ and the 74.5-day to W'®. However, 
from the deuteron bombardment of rhenium, a weak 
activity of 70+6 days is found in the tungsten separation, 
no 24-hour activity being present. Such an activity could 
be produced by a Re!®? (d,a)W!® reaction. This additional 
evidence supports the conclusions of Minakawa, thus the 
following scheme is most probable. 


180 182 183 
0.2% 22.6% 17.3% 


184 185 186 187 
30.16% 74.5days 29.8% 24hr. 


In conclusion, we wish to thank Professor J. M. Cork 
for his active cooperation during all phases of the work 
and we are very grateful to others of the Physics and 
Chemistry Departments who gave assistance, especially 


to Dr. J. L. Lawson. 
KASIMIR FAJANS 


WILLIAM H. SuLtivan 


Chemistry Department, 
University of Michigan, 
Ann Arbor, Michigan, 
July 2, 1940. 


1 F. W. Aston, Proc. Roy. Soc. A132, 487 (1931). 

2 E. Amaldi, E. Fermi, et al., Proc. Roy. Soc. Al49, 522 (1935). 

31. V. Kurtchatow, et al., Physik. Zeits. Sowjetunion 8, 589 (1935). 

4 H. Yamasaki and K. Sinma, Sci. Pap. Inst. Phys. Research (Tokyo), 
37, 10 (1940). 

5 O. Minakawa, Phys. Rev. 57, 1189 (1940). 
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On the Efficiency of y-Ray Counters 


The efficiency curve for y-ray counters (efficiency 
defined as the ratio of the number of incident quanta to 
the number of registered quanta) is a valuable key for the 
interpretation of every counting experiment involving 
y-rays, but a careful determination of such curves meets 
with some difficulty. It implies an accurate knowledge of 
the number of quanta emitted in some transformations, 
but only a few (if any) energy-level schemes involving 
y-rays are quite definite at present. The coincidence 
method is one of the best ways of testing such schemes, 
but it, in turn, makes use of efficiency curves. Because of 
the low efficiency of y-ray counters, whereby counting 
experiments often become laborious, it is also of interest 
to see whether different materials show different efficiencies. 

Counters of brass, aluminum and lead, all of the same 
geometric size (effective length 6 cm, inner diameter 1.4 
cm, wall thickness 0.2 cm), have been compared (Fig. 1). 


% A 
oul + Fb 
Brea 
$40 
Nas ' 


40 25 30 AV 
5030 20 s 


Fic. 1. Comparison of efficiency of counters of different materials. 


The curve for brass was published by Droste,! who 
measured the end points and estimated the form semi- 
empirically; it was recently verified also by Dunworth.? 
Yukawa and Sakata® have calculated the shape of the Al 
curve. The point reported here for 2.7 Mev will thus fix 
this curve. The measurements at 2.7 Mev were made by 
a radio-thorium source of known strength (ThC” y-rays, 
filter 2.5 cm Pb). Points on the Pb and Al curves at about 
0.9 Mev were determined by 8—vy-coincidence measure- 
ments of *Mn** with the source surrounded by a 0.24 
g/cm? brass filter to eliminate the low energy 8-group. 
According to Curran, Dee and Strothers‘ the high energy 
8-group should be coupled chiefly to a 0.9-Mev line. At 
still lower energy *Au™, having y-components of 0.28 Mev 
and 0.44 Mev,® was used. The 8—y-coincidence rate cor- 
responds to an efficiency of about 1.1 percent for the 
0.44-Mev line, if the 8-spectrum is simple, which is 
plausible. In any case the efficiency of the lead counter 
for the 0.44-Mev line was about 5.6 times greater than the 
efficiency of the brass counter (after correcting for absorp- 
tion). This fact, which must be related to the large photo- 
effect in lead at these energies (see Fig. 1b where 7 is the 
photoelectric absorption coefficient, R the range of photo- 
electrons), makes the lead counter very effective in the 
region of 0.5 Mev, where the other counters have low 
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2 8 1/0 mm Pb 


Fic. 2. Comparison of absorption of y-rays and coincidence rate. 


efficiency. Therefore, a brass counter with thin lead foils 
on the inside probably should be most useful. For RaC 
y-rays, filtered through 2.5 cm Pb, (~2 Mev) the relative 
efficiencies were: Pb=100, brass=84+4, Al=109+4, 
and for unfiltered Ra y-rays correspondingly: Pb = 100, 
brass =67+3, Al=75+3. 

It may be pointed out, that a y-component in favorable 
cases (highly different absorption coefficients of the 
components) can be correlated to other particles without 
use of an efficiency curve by measuring the decrease of 
the coincidence rate when absorbers are inserted in front 
of the y-ray counter. It seemed worth while to try this 
method in the important case of *Mn** above. In Fig. 2 
the upper curve is a simple absorption curve (Pb-counter), 
whereby the bad geometry may be taken into account. 
The coincident rate curve below it (8-filter 0.24 g/cm? 
brass) gives some support to the assumption of Curran, 
Dee and Strothers, that the most energetic 8-particles are 
coupled to the soft y-component, but the curve is not quite 
conclusive, since a still softer y-component of low intensity 
(possibly “Bremsstrahlung”’) seems to be present. 

Because of the higher effectivity of obliquely incident 
rays the coincidence measurements have been reduced to 
a “large distance” between source and counter by a 
distance function: 


R 
F(x, k) (1) 
ab 


and not by the relative space angle, approximately given 
as: 


out Rb (2) 
arctan 
(2R internal diameter, a the distance from source to counter 
axis, F elliptic integral of first kind, / effective counter 
length, x=arctan (//2c), b=a+R, c=a—R, 
The difference between these formulas is appreciable only 
for small distances such as are common in coincidence ex- 
periments. For high y-energy and light elements the 
formula (2) should be best (Compton recoil forward). 
A more detailed account will be published later else- 
where. 
I wish to thank Professor M. Siegbahn and Professor 
L. Meitner for their interest in this work. 
FoLke NorLinG 
Forskningsinstitutet for Fysik, Vetenskapsakademien, 


Stockholm n, 
June 4th, 1940. 


1G. A Droste, Zeits. f. Physik 104, 474 (1937). 
: Dunworth, Rev. Sci. Inst. 11, 167 (1940). 
: i. ‘aoe and 8. Sakata, Sci. Pap. Inst. Phys. and Chem. Research 
(Japan) 31, 187 (1937). 
4S. C. Curran, P. I. Dee and J. E. Strothers, Proc. Roy. Soc. 174, 


553 (1940). 


‘J. R. Richardson, Phys. Rev. 55, 609 (1939). 


277 
50| rete 
30 
dio- 
and | 
urs | 
vith 
mit 
1.44 
pre- 
).87 
ited 
ray | 
om 
om 
tive 
low 
ver, 
eak 
uld | 
nal 
the 
| 
ork | 
ork 
and | 
ily | 
| 
| 


278 ' LETTERS TO THE EDITOR 


Pair Production of Mesotrons at 29,000 Feet 


In a second airplane flight! made on April 30, carrying 
a cloud chamber in a magnetic field of 700 gauss up to 
29,000 feet altitude, we found one particularly interesting 
photograph which shows the occurrence of a mesotron 
pair. A reproduction is given in Fig. 1. 

The radii of curvature of the two partner tracks are 45 
and 75 cm. Both particles are stopped in a copper plate of 
0.32 cm which was placed in the middle of the chamber. 
They enter this plate at an angle of 45°, which gives a 
maximum path length inside the copper of 0.45 cm. If the 
particles were electrons and were stopped in this amount 
of copper their maximum energy is 6 Mev. This corresponds 
to a radius of curvature in the magnetic field smaller than 
28 cm. The measured radii are larger by an amount far 
beyond the errors of measurement. Thus the tracks cannot 
be produced by electrons. For mesotrons of a mass 200, 
similar considerations give 25 Mev as the maximum 
kinetic energy of the particles and a radius of curvature 
smaller than 350 cm. This curvature limit covers well the 
measured values. The same holds true, of course, if one 
considers the particles to be protons. 

One can further estimate the masses from the ionization 
produced in the gas of the chamber. Fortunately the 
picture shows two slower particles and a fast one. Their 
ionization is markedly smaller and corresponds to the kind 
of tracks we get from 8-particles and penetrating cosmic 
rays in the laboratory. There is no doubt that the par- 
ticles of small radii of curvature are electrons. The ioniza- 
tion of the two particles in the pair is distinctly higher 
although the radii of curvature are greater. This proves 
again that the pair must consist of heavier particles. 

For the observed curvature of the pair, mesotrons should 
ionize from 6 to 15 times as strongly as electrons, and 
protons 10 times as strongly as mesotrons. Such an enor- 
mous ionization excludes positively the possibility of 
assuming that the pair is protons. The theoretical ioniza- 
tion of mesotrons, on the contrary, lies well within the 
limits one estimates from the pictures. 

We therefore must conclude that this pair consists of a 
slow positive and a slow negative mesotron. 

The two stereoscopic pictures taken show that the slow 
electron emerging upward from the copper plate is possibly 
associated with this pair and may be a decay electron of 


one of the mesotrons. It is difficult to state from which 
mesotron this electron might come. Assuming that jt 
moves upward from the plate its charge is negative. [t 
has a kinetic energy of 0.2 Mev. 

Other pictures taken on this flight likewise show meso- 
trons stopping in the copper plate. In several cases electrons 
appear which may originate from mesotron decay. A fyl] 
report on this material will shortly be published. 

G. HERzoG 


The Texas Company, 
Houston, Texas, 


W. H. Bosticx 


Ryerson Physical Laboratory, 
University of Chicago, 
Chicago, Illinois, 
July 15, 1940. 


1G. Herzog, Phys. Rev. 57, 337 (1940). 


Transport Cross Section of Helium 


Endeavors to fit the low temperature viscosity measure- 
ments of helium! have in the past not yielded very satis- 
factory results. Use of the Slater-Kirkwood field,? the 
Lennard-Jones potential,’ the elastic sphere‘ and the 
modified elastic sphere models® have all failed in predicting 
correctly the quantum effect at low energies. It has been 
pointed out by the author that the general behavior of the 
transport cross section which is required for a correct 
temperature dependence is: the existence of minima at 
low energies, a broad maximum at higher energies, and at 
energies corresponding to room temperatures an asymptotic 
approach to a value approximating 4.46 X 10~"* cm?, 

The physical characteristics of the interaction under- 
lying this behavior are that two elements of differing order 
of magnitude must exist: a strong repulsive potential with 
a rather sharply defined boundary, and a weak interaction 
composed of attractive and repulsive terms. It is essential 
for the existence of minima in the low energy range that 
the interaction contain attractive terms, as, although the 
elastic sphere and modified elastic sphere models show 
minima in the kinetic cross section, these disappear in the 
viscosity cross sections due to the presence of two com- 
pensating series. 

It has seemed of interest therefore to determine as 
exactly as possible the form of the cross section necessary 
to fit the experimental data. That is, one might attempt to 
solve the integral equation for @(y7) the cross section, using 


the relation given by Enskog® 


15 (mkT)} 
7 
where 
k* 82x2mE 
*)2— 
2x (mkT)¥ (k*) (2) 


The exact solution of the above equation presupposes 
the knowledge of (7) for all values of T. One has, however, 
only its values for a limited set of values of the temperature. 
The solution of @(7) is then subject to a certain degree of 
indefiniteness. It has, therefore, seemed most advantageous 
to assume for @ the expansion 


0(k*) A; exp [—j(k*)*o] (3) 
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and so determine A; so as to fit the experimental points. 
The arbitrary choice of the parameter ¢ (unless deter- 
mined in the same manner as the A) introduces the inde- 
terminacy of definition previously mentioned. Retaining 
the idea of a hard core, the term for 7=0 was included in 
the sum although the behavior of 9(7) does not necessarily 
justify this step. 

In carrying through the calculation it was found possible 
to fit the points with better than a one percent accuracy 
by retaining only the first six terms of the sum. 

The values of @(k*) are given in Table I for values of 
k*X X2. 

TABLE I. Values of 0(k*) for various values of k*. 


0.00 9.39 4.52 3.82 8.02 2.75 
2.05 4.41 4.89 3.72 9.60 3.04 
2.55 3.62 5.26 3.51 9.81 3.04 
2.99 3.41 6.05 2.96 10.05 3.03 
3.39 3.49 6.48 2.74 10.61 2.96 
3.78 3.67 6.94 2.63 12.51 2.58 
4.15 3.79 7.45 2.63 16.63 2.25 


The results reveal the existence of two minima. In 
general the trend is that of the elastic sphere cross section 
for S=2.11X10-* cm. The values of @ in the neighborhood 
of k*->0 permit large variations due to the form of (1) 
and hence for both k->0 and k=, must be regarded as 


indeterminate. 
E. J. HELLUND 
Department of Physics, 
University of Washington, 
Seattle, Washington, 
July 10, 1940. 


1A. v. Itterbeek and W. Keesom, Aye) - 257 (1938); A. v. 
Itterbeek and O. v. Paemel, Physica A 265 (1940 

2H. S. _ Massey and R.A . Buckingham, i. Roy. Soc. A168, 
378 (1938 

4H. S. W. . Proc. Roy. Soc. Al4l1, 434 
(1933); Alaa, 188 (1934). 


SE. J. Hellund (unpublished). 
*D. Enskog, * ‘Kinetische Theorie der Vorgange in Massig Verdiinnten 
Gasen” (Dissertation, Upsala, 1917). 


Electric Breakdown of Alkali Halides 


In discussions of the electric breakdown of alkali halide 
crystals an expression for the force between a free electron 
in the crystal and an ion in the lattice has proved to be of 
importance. Such a force is of the Coulomb type, of course, 
but it also contains a factor that depends on the dielectric 
constant and other properties of the lattice. Designating 
this factor by the symbol e«, we may express the force 
between a free electron of charge e and an ion of 2-fold 
electronic charge by eze*/r?, where r is their distance of 
separation. In the case of an electron that is moving suf- 
ficiently slowly, the following formula for ¢ has been given! 


in which yu is the “reduced mass” of the ions, a is the 
distance between neighboring ions, k’ is the dielectric 
constant extrapolated from high (visible) frequencies, and 
k is the actual dielectric constant. The frequency » which 
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appears here is one appropriately averaged over the crystal 
vibrations that are excited by the passage of the electron. 
In making such an average, vibrations of rather long wave- 
length have considerable weight because of the interaction 
between the electron and somewhat distant ions. Since 
transverse vibrations are not excited by an electric field 
with VX E=0, only longitudinal ones need be considered 
here. According to Mott and Littleton* the frequencies of 
the longitudinal vibrations differ considerably from those 
of the transverse ones, the observed values of which 
have previously been employed in the calculation of 
breakdown fields. For in a longitudinal vibration in which 
positive and negative ions are displaced by unit length 
relative to each other, a surface charge density accumulates 
near the nodes, viz., +e*/a*. (The formulae used in this 
letter include the effect of electronic polarizability in a 
rough approximation; they are not based on the customary 
assumption of a full electronic charge rigidly attached to 
an ion.) Here e* is the apparent charge attached to an 
ion; it may differ from the electronic charge on account 
of the distortion of the ions. This surface charge at the 
nodes gives rise to an additional restoring force, which is 
equal to 2re**/a*k’ for longitudinal waves and which is 
altogether absent for transverse ones, there being no 
accumulation of charge in the interior of the crystal in 
this case. We find, therefore, the relation 

for the frequency v, of the longitudinal vibrations in 
terms of the frequency vy of the transverse ones. (Our 
formulae become equal to those of Mott and Littleton for 
k’=1 and for e*=e.) It can be shown that the transverse 
frequency satisfies also the following equation of Born 
and Géppert-Mayer,’ which has only been modified by 
the replacement of e with e*: 


=e**/(k—k’). 
From these two expressions we obtain 
Finally, the factor « becomes 


Using this value for « and substituting vr for » throughout 


TABLE I. Comparison of calculated and Speed ontuns of the breakdown 
potentials in alkali hal 


or 


ALKALI 

HALIDES Ay E Bg-calc EB-obs 

LiF 18.9%u 58.0 (10)§ v/cm 31.0 (10)§ v/em 
aF 41.3* 22. 24 

NaCl 61.1 7.2 15 

NaBr 74.7 3.5 8.1 

KCl 70.7 5.3 10.0 

KBr 88.3 3.7 7.0 

KI 102 2.7 5.7 

RbCl 8 4.8 8.3 

RbBr 118 * 2.8 6.3 

RbI 130 * 2.3 4.9 
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the derivation for the formula for the breakdown field 
Esa we find 


1.82x2mev,/1 1\/k\* 


where m is the so-called ‘effective mass” of the electron 
in the lattice. 

In the second column of Table I we give the values of 
the wave-lengths \, corresponding to the frequencies »7, 
which we have used in the calculations of the breakdown 
fields of the alkali halides listed in the first column. In 
general, the values are those that have been obtained from 
absorption measurements in thin crystal layers; those 
marked with an asterisk, however, are the “reststrahl” 
wave-lengths multiplied by the number 1.10 (the factor 
1.10 being due to the fact that the wave-length measured 
by the method of maximum reflecting power does not 
coincide with that of the lattice vibrations). In the third 
and fourth columns, respectively, the values of the break- 
down fields calculated by means of the new formula are 
compared with the values originally observed. The com- 
putation has been carried out with the substitution of the 
mass of a free electron for the mass m in the formula. 


THE EDITOR 


The deviations between the calculated values and those 
observed cannot be ascribed merely to the coarseness of 
the theoretical approximations, but may be due even more 
to the omission of several physical factors, of which two 
important ones are the following: (1) The “effective mags” 
of the electron, which should be used in the calculation, 
if known, differs from the mass of a free electron; it may, 
indeed, be greater or less than this mass. (2) Breakdown 
fields calculated for absolute zero temperature, as in this 
case, are necessarily lower than those observed at room 
temperature (on account of a temperature dependence). 
For example, in the case of KBr the calculated value of 
the breakdown field at 0°K is 3.7 (10)§ v/cm as compared 
with the observed one of 7.0 (10)* v/cm at 300°K, whereas 
the measured value at 200°K is only 5.3 (10)* v/cm, 

R. J. SEEGER 
E. TELLER 
July 3, 1940. 


and E- Teller, Phys. Rev. 54, 519 
Mott and M. J. Littleton, Trans. Faraday Soc. 34, 485 (1938); 
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4R. C. Buehl and A. von Hippel, Phys. Rev. 56, 946 (1939). 


i 
{ 
1 
t 
ri 
h 
Vv 
ce 
0. 
of 
a] 
at 
pu 
| 


